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Supersonic Civil Airplane Study and Design: 
Performance and Sonic Boom 

Samson Cheung 


This final report summarizes the work performed from July 1989 to 
Jan. 1995. The work is supported by NASA Co-operative Agree- 
ment NCC2-617. This report consists of four parts. The first part is 
the introduction of the research effort. The second part describes 
the work and results from July 1989 to June 1993. The third part 
describes the work and results from July 1993 to January 1995. A 
summary is given at the end of this report. 


1 INTRODUCTION 

The present supersonic civil airplane, the Concorde, is a technological break-through in 
aviation history. However, it is an economical disaster for two main reasons. The first is 
her low aerodynamic performance, that allows only 100 passengers to be carried for a 
short-range flight with expansive airfare. Another reason is that the shock waves, gener- 
ated at supersonic cruise, coalesce and form a classical N-wave on the ground, forming a 
double bang noise termed sonic boom, which is environmentally unacceptable. To 
enhance the U.S. market share in supersonic civil transport, an airframer’s market risk for 
a low-boom airplane has to be reduced. 

Since aircraft configuration plays an important role on aerodynamic performance and 
sonic boom shape, the configuration of the next generation supersonic civil transport has 
to be tailored to meet high aerodynamic performance and low sonic boom requirements. 
Computational fluid dynamics (CFD) can used to design airplanes to meet these dual 
objectives. The work and results in this report are used to support NASA’s High Speed 
Research Program (HSRP). 

In this five years of study and research, CFD tools and techniques have been developed for 
general usages of sonic boom propagation study and aerodynamic design, hi the beginning 
of the 90’s, sonic boom extrapolation technique was still relied on the linear theory devel- 
oped in the 60’s for the nonlinear techniques were computationally expensive. A fast and 
accurate sonic boom extrapolation methodology (Section 3.2), solving the Euler equations 
for axisymmetric flow, has brought the sonic boom extrapolation technique up to the 90’s 
standard. 

Parallel to the research effort on sonic boom extrapolation, CFD flow solvers have been 
coupled with a numeric optimization tool to form a design package for aircraft configura- 
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tion. This CFD optimization package has been applied to configuration design on a low- 
boom concept (Section 2.3) and an Oblique All- Wing concept (Section 2.4) prior to the 
wind-tunnel models are built and tested at Ames. The tunnel test results have validated the 
CFD technique and design tools. 

Moving to the world of parallel computing, the aerospace industry needs a numeric opti- 
mization tool suitable for parallel computers. A nonlinear unconstrained optimizer for Par- 
allel Virtual Machine has been developed for aerodynamic design and study. Study in 
Section 3.3 demonstrates the capability of this optimizer on aerodynamic design. 


2 PREVIOUS WORK/RESULTS 

The work and results described in this section was begun in July 1989. The first project 
was to use CFD tools and existing linear theory to predict waveform signatures at some 
distances from flight vehicles. The aim of this study was to demonstrate and develop the 
technique of sonic boom prediction by CFD. The next step was to apply this developed 
technique to low-boom configurations. 

The second project, which was the continuation of the first one, was to develop a CFD 
optimization package for design process on meeting the dual objectives of high aerody- 
namic performance and low sonic boom loudness. This optimization package was applied 
to three different High Speed Civil Transport (HSCT) baseline configurations and a 
generic body of revolution. 

A wind-tunnel model (Ames Model 3) was built based on one of the modified HSCT base- 
line configuration. This model was tested in June 1993. The test results were used to vali- 
date the design method. Publication of the result was limited due to the sensitive nature of 
the project. 

A counterpart of the conventional HSCT concept was the Oblique All-Wing (OAW) con- 
cept. CFD computational supports, as well as optimization calculations, were provided to 
the OAW design team consisting personnels from NASA Ames Research Center, industry, 
and university. The aim of the project was to design a realistic configuration for wind-tun- 
nel test. The model was built and tested at Ames in June 1994. 


2.1 Sonic Boom Prediction Technique 

In the early stage of sonic boom prediction activity, two major things were involved. The 
very first thing was to identity the capability of CFD in sonic boom prediction. The second 
thing was to apply these CFD tools to predict sonic boom signals of varies configurations 
after necessary code modification, grid refinement study, and comparison with supersonic 
linear theory. 
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2.1.1 Method Validation 


A three-dimensional parabolized Navier-Stokes code, UPS3D,* developed at Ames was 
used as the flow-solver. It is a space-marching code with finite-volume approach. The near 
field solution of a simple wing/body configuration was calculated by UPS3D, and the 
overpressure signal at some desired distances were obtained either by the axisymmetric 
option of UPS3D or a quasi-linear extrapolation code, based on Whitham’s F-function the- 
ory 2 . Later I realized that using Lighthill integral 5 to calculate the F-function for non-axi- 
symmetric aircraft was more accurate, I wrote a Fortran code, OIF, for sonic boom 
prediction based on Lighthill integral. This code is available from Ames Software Library. 
A copy of LHF is attached in Appendix A. The figure below is a brief summary of the 
sonic boom extrapolation process. 


CFD Calculation 



A series of studies on grid refinement, including solution adaptive grid, and on sensitivity 
of initial distance of extrapolation were conducted. It was found that viscous calculation 
was unnecessary for sonic boom prediction. However, the grid must be sufficiently fine in 
the regions of shock and expansion waves. In order to capture all the nonlinear effects in a 
three-dimensional flow, the near-field overpressure should be captured at about one span 
length below the flight track before extrapolating to the far field. The detail results were 
published in AIAA Journal of Aircraft 4 and NASA Technical note 5 . 

In summary, the tools for sonic boom prediction had been identify and validated in the 
above study. The combination of CFD and Whitham’s method gave a relative efficient 
tools for sonic boom prediction. Nevertheless, the CFD codes was still computationally 
expensive for design optimization runs. 
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2.1.2 Boom Prediction for Low-Boom Configurations 

With the experience on grid refinement study and the extrapolation procedure, the predic- 
tion tools were being used to predict the sonic boom of two low-boom configurations 
designed by Boeing aircraft company and Langley research center. 

Each of the two configurations consisted of two separated parts, namely, the wing and the 
fuselage. The wing was defined by data in spanwise cuts, whereas the fuselage was 
defined by data in streamwise cuts. In order to create a single wing-fuselage surface gird 
for UPS3D code, a grid generator (SAMGRID) was written to defined the wing in stream- 
wise cuts and aggregated the wing to the fuselage. Computation results of the two config- 
urations are shown below. 


Langley’s Mach 2 Configuration 



The sonic boom signals calculated from the CFD prediction tools were compared to the 
wind-tunnel data of the Langley’s configuration. The computational results of the Boe- 
ing’s configuration was used to validate the linear design method used by Boeing. 


2.2 Supersonic Airplane Design 

The need for simultaneous sonic boom and aerodynamic optimization was highlighted 
when it became clear that designed to a strict sonic boom constraint suffered an unaccept- 
able performance penalty. Therefore, low-boom design studies must carefully balance the 
trade-off between sonic boom loudness and aerodynamic performance. A CFD optimiza- 
tion package was developed to demonstrate the methodology for the optimization of 
supersonic airplane designs to meet the dual objectives of low sonic boom and high aero- 
dynamic performance. 

In this project, an optimizer with linear and nonlinear constraints was first identified, and 
then an efficient CFD flow solver was chosen. This CFD code had to be sufficiently fast 
because more than 90% of the computational time were used in CFD calculations. Before 
this optimization was used to design low-boom wind-tunnel model (Section 2.3), it was 
tested and exercised by improving aerodynamic performance of a low-boom wing/body 
configuration and a body of revolution. 
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2.2.1 CFD Optimization Package 

Several computational tools interconnect in the optimization procedure are listed below: 

• UPS3D: 3-D parabolized Navier-Stokes code; inviscid calculation only (Ref. 1) 

• NPSOL: numerical optimization code 6 ; a sequential quadratic programming algorithm 
in which the search direction is the solution of a quadratic programming subproblem 

• HYPGEN: hyperbolic grid generator 7 ; a sufficiently fast and robust to operate within 
an automated optimization environment. 

■ LHF: sonic boom extrapolation code (Appendix A); a routine based on Whitham's F- 
function and the equal-area rule 8 

• SAMGRID: wing/body surface grid generator (Appendix B); a sufficiently fast and 
robust to operate within an automated optimization environment 

■ DB: sonic boom loudness calculation; a code gives perceived loudness (PLdB) of the 
sonic boom can be determined by Stevens’ Mark VII method 9 which involves Fast 
Fourier Transform on the energy spectrum of the sonic boom 

This CFD optimization package is robust and efficient on Cray-YMP. The application of 
this package will be described in the following sections. 

2.2.2 Aerodynamic and Sonic Boom Optimization 

The optimization design package was exercised using a recently-developed low-boom 
wing-body configuration, Boeing 1080-991 (also called Haglund model), designed by 
George Haglund. This optimization technique was applied separately to the two objectives 
of high aerodynamic performance and low sonic-boom loudness. 

For aerodynamic enhancement, control points are set on the cambers of the wing, with the 
thickness kept fixed. The left figure below shows the differences on a inboard airfoil sec- 
tion of the original and the modified. The polar plot shows the improvement of L/D of the 
modified configuration over the original by 3.8%. The right figure below shows that the 
modified wing had less wave drag than the original one at the leading edge. This means 

Aerodynamic Optimization ofHSCJ Wing Camber at Mach 1.7 
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that the leading thrust is improved by the optimization process. The whole process takes 
about 4 CPU hours on Cray-YMP. 

For sonic boom improvement, F-function was employed as an entity to define the equiva- 
lent area distribution and sonic boom shape. The original Haglund model was supposed to 
give a flat-top pressure waveform at the ground. However, calculations showed that the 
waveform had an intermediate shock followed right after the bow shock; whereas the flat- 
top waveform would have no intermediate shock. The design code redistributed the equiv- 
alent area of the fuselage (without changing the wings), and re-captured the flat-top char- 
acteristic of the pressure waveform. The figure below compares the sonic boom signatures 
among the original, optimized, and target flat-top. Due to the sensitive nature of the con- 



figuration, the change of the configuration will not be shown here. The details of this opti- 
mization methodology and results were considered as sensitive materials and were 
presented in the 2nd Annual Sonic Boom Workshop. 10 

2.23 Drag Minimization on Haack-Adams Body 

The purpose of this study was threefold: 

■ to search for a design method to minimize the drag of a supersonic projectile 

• to demonstrate the capability of the CFD optimization package described above 

• to search for computational grid density effect on optimization performance 

The baseline configuration chosen for this study was called Haack-Adams body 11 , a body 
of revolution with a pointed nose and a base of finite area. This body was thought to be the 
minimum-drag body under the slender body theory. Wind-tunnel data were available for 
CFD validation. The method of optimization made use of the Fourier Sine expansion, 
which had three main advantages over the traditional techniques based on shape functions 
and control points: 

• The volume of the body was fixed without putting external constraints. External con- 
straints cost more computational time. For some cases, fixed volume is not feasible. 

• Global minimum was search. 

• Number of design variables was substantially reduced. 
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The figure below summaries the result of this study. The nose of the body was trimmed to 


Variation of Wave Drag Coefficient 
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reduce the wave drag. Since the total volume was constrained, volume was added near the 
end of body. Total wave drag reduction was by 6%. The results were presented in a AIAA 
meeting 12 and published in Journal of Aircraft Vol. 32, No. 1, Jan/Feb. 1995. 


2.3 Low-Boom Wind-Tunnel Configuration (Ames Model 3) 

Efforts were made to design a new wing/body /nacelle configuration, which had a lower 
sonic boom relative to the baseline, 1080-911 from Boeing Company, of low boom HSCT 
concept. The CFD optimization package described in Section 2.2.1 were employed to 
modify this baseline configuration. The result of the optimization was used to build a 
wind-tunnel model, Ames Model 3, tested at Ames 9’x7’ wind tunnel in June 1993. Due to 
the sensitive nature of the configuration, no planform shapes will be shown here. How- 
ever, the left and right figures below show the computational grid and the optimization 
result, respectively. The plot at the lower right-hand comer of the right figure shows the 



sonic booms of the baseline and Model 3 respectively. The baseline configuration has a 
loudness level about 100 PLdB; whereas Model 3 has about 92 PLdB. The results of this 
research were presented in the 3rd Annual Sonic Boom Workshop. 13 
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2.4 Oblique-All Wing (OAW) Computation and Design 

Oblique flying-wing 14 is an alternative supersonic aircraft concept. Ames, Boeing, Dou- 
glas, and Stanford University joined and formed a design team in 1992 to investigate the 
feasibility of OAW for commercial use. The study included aerodynamic performance, 
stability, structure, landing gear, airplane exits, and airport regulations. The design team 
decided to build a wind-tunnel model for wind-tunnel testing in June 1994. My job was to 
provide Navier-Stokes CFD supports and, if possible, optimization results. The figure 
below shows some of the wings that were analyzed since the beginning of this study. 


OAW Design History 



Time line 


The flow solver being used was Overflow code, a 3-D Navier-Stokes code using the diag- 
onal with ARC3D algorithm 1 ^. One of the most challenging works of this project was to 
reduce the separation on the left wing (trailing wing). TTie separation on the upper surface 
of the wing and the corresponding vortices are shown in the left side of the figure below. It 


OAW-O OAW-O (beading 



was found that bending of the wing could abate the separation, as well as improve the lift- 
to-drag ratio. The right side of the figure shows a weaker separation pattern on the ended 
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wing. Due to the sensitive nature of this study, the results can only be presented in the 
weekly group meetings at Ames and a controlled distributed NASA Contractor Report. 


3 CURRENT WORK/RESULTS 

Currently, research effort was concentrated on one theme that is sharpening the tools for 
HSCT design. Three research topics are focused: near-field CFD calculation and sonic 
boom softening of Boeing Reference-H, improvement of sonic boom extrapolation, and 
aerodynamic design on parallel computer. 

In order to study and design a real complex aircraft, a relatively fast CFD technique has to 
be developed for optimization environment. Coupling a fast space-marching code and a 
time iterative code with overset gird concept can take the advantage of marching code at 
the fuselage/wing region and solve the complex flow field near the wing/nacelle region at 
the same time. 

A very efficient wave propagation code for mid-field sonic boom prediction has been 
developed based on the method of characteristics. This code solves the Euler equations for 
1.2 minutes on Cray-YMP; whereas, the axisymmetric CFD method described in Section 
2.1.1 takes 40 minutes on the same computer. 

Number crunching problems, like CFD calculations, on parallel machines can be effi- 
ciently done in today’s computing environment. This may lead to the future of aerody- 
namic research and design. In order to exercise HSCT design on parallel computers, a 
nonlinear optimization routine has been developed for a network based parallel computer 
system in which a cluster of engineering workstations serves as a virtual parallel machine. 


3.1 Sonic Boom and Performance Study of Reference-H 

Research effort on low-boom configuration concept has been invested for the past four 
years. A new proposed route structure for HSCT’s incorporating supersonic corridors over 
land and water has relaxed the sonic boom constraint somewhat. The objective of this 
study is twofold. First is to exercise the methodology of combining two different CFD 
codes to solve the near-field solution of a realistic HSCT configuration in an efficient and 
accurate manner. Second is to reduce the sonic boom loudness of a performance configu- 
ration concept, Reference-H, without jeopardizing the aerodynamic performance. The 
basic components of Reference-H are a fuselage, a pair of swept wings, and four nacelles. 

3.1.1 Reference-H Near-Field Study 

The CFD codes used in this study are the UPS3D code and the OVERFLOW code. Both 
CFD codes has been described in Section 2.1.1 and 2.4, respectively. The former is an effi- 
cient space-marching code. However, it fails in the region where subsonic pocket exists; 
especially in the region of the wing/nacelle integration. The latter is a time-iterative code 
with Chimera overset grid concept, which makes the code more viable in solving the 
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region of wing/nacelle integration. In this study, only inviscid flow is considered. Figure 
below summarizes the result of the CFD calculations. 



The near-field solution is studied for the case of Mach number 2.4 and angle of attack 4.5 
degrees. Wind-tunnel data of the Reference-H validate the CFD method. Study shows that 
flow particles turn significantly over the outer nacelle compared with the inner nacelle. It 
indicates that the effect of the nacelle orientation might improve the aerodynamic perfor- 
mance. 

3.1.2 Sonic Boom Softening 

The sonic boom of the Reference-H configuration is also obtained. The calculation shows 
that the boom is an N-wave of 104 PLdB with 2.5 psf. bow shock on the ground. Details of 
the sonic boom prediction technique can be found in Ref. 10. Boom modification for per- 
formance aircraft is very much different from the low-boom aircraft for cruise Mach num- 
ber and lift are higher. Therefore, the technique developed previously can not be strictly 
applied to Reference-H. However, changing the equivalent area can be helpful. The result 
of this study was presented in the 4th Sonic Boom Workshop. 17 Another approach to 



reduce the boom is by experimenting the sweep angle. The figure above show one of the 
exercises done on the Ref-H. This exercise successfully shows Boeing how much boom 
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reduction can be achieved by redistributing the lift. An closer on-going technology com- 
munication with airframe industry is needed in order to achieve the goal of sonic boom 
softening on performance aircraft. A team consisting myself and other personnels from 
Boeing and NASA Langley has been formed to achieve the goal. 


3.2 Sonic Boom Mid-Field Extrapolation (WPSYM) 

In the beginning of 90 ’s, sonic boom extrapolation technique was still relied on the linear 
theory developed in the 60’s for the nonlinear techniques were computationally expensive. 
Today, a fast and accurate sonic boom extrapolation methodology is needed to bring the 
sonic boom extrapolation technique up to the 90’s standard for HSCT design. The objec- 
tive of this study is to develop an efficient and accurate higher-order computational 
method, solving the Euler equations, for supersonic aero-acoustic wave propagation. 

An axisymmetric wave propagation code (WPSYM) has been developed for mid-field 
sonic boom extrapolation. This propagation code has been demonstrated as an efficient 
and accurate tool over the previous CFD method, described in Section 2.1.1 and Ref. 4, on 
a generic wing-body configuration. The figure below shows that a 3-D near-field solution 



AX/ L 

is obtained from UPS3D code; the result is then interfaced to two axisymmetric sonic 
boom extrapolation codes, namely, the axisymmetric version of UPS3D and the recent 
wave propagation code (WPSYM). The former takes 40 minutes on Cray-YMP, and the 
latter takes 1.2 minutes on the same machine. The x-y plot in the figure compares the 
numerical extrapolation results to wind-tunnel data. The result has been shown in NASA 
Technical Highlight and the methodology has been presented in the 4th Annual Sonic 
Boom Workshop at NASA Langley in June 1994. 
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3.3 Optimizer on PVM (7/OWA) 

Moving to the world of parallel computing, the aerospace industry needs a numeric opti- 
mization tool in the parallel environment. One of the promising parallel computing con- 
cept is the network-based distributed computing. The Parallel Virtual Machine (PVM) is a 
software package that allows a heterogeneous network of parallel and serial computers to 
appear as a single concurrent computational resource. PVM allows users to link up engi- 
neering workstations to work as a single distributed-memory (parallel) machine. Merritt 
Smith and I wrote a manual on PVM for beginning users. A copy of the manual is attached 
in Appendix C. 


A parallel optimizer based on nonlinear Quasi-Newton method has been developed and 
coupled with an efficient CFD code for basic aerodynamic design and study. This opti- 
mizer is called 7/OWA (parallel Optimizer With Aerodynamics). The figure below is a 
demonstration of 7/OWA. A Boeing arrow wing/body configuration is chosen in this 
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study. The fuselage radius is changed so that the wave drag is minimized. The parallel 
CFD optimization process takes 24 wall-clock hours on 4 SGI workstations to reduce the 
wave drag by 6.5%. The optimized result is a “coke bottle” shape fuselage, as expected by 
supersonic area rule. The convergence history of the optimization process is also shown in 
the figure. The optimizer is also coupled with a parallel CFD code, MEDUSA, to perform 
viscous 2-D multizone airfoil optimization supported by overset grid concept. The results 
will be presented at NASA CAS conference in March 1995. 


3.4 Oblique All- Wing (OAW): CFD support 

The OAW design team has asked for CFD support on the latest configuration OAW-3 from 
which a wind-tunnel model has been built and tested at Ames in June 1994. The figure 
below shows the chimera grid topology on the OAW-3 with fin. The design team want to 
compare the CFD result with the result from pressure sensitive paint (PSP). Therefore, 
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CFD calculations have to be done prior to the wind-tunnel test because color map from 
CFD result is need for PSP calibration. 


4 SUMMARY 

The computational tools for sonic boom prediction, aerodynamic calculation, and configu- 
ration design of the current HSCT concept have been validated and applied to build wind- 
tunnel model for further testing and validation. The techniques developed in this five-year 
research and their applications, such as sonic boom prediction technique (Section 2.1), 
design of Ames Model 3 (Section 2.3) by CFD optimization (Section 2.2), and sonic boom 
softening for performance configuration (Section 3.1), have clearly shown support to the 
HSRP as it moved to its phase two period. 

An accurate sonic boom extrapolation tool has always been an issue. It is because the flow 
phenomena in the atmosphere are nonlinear, but the common technique for extrapolation 
is linear acoustic theory developed in the 60 ’s. On the other hand, CFD technique is too 
computationally expensive. Recently, a fast and accurate sonic boom extrapolation meth- 
odology (Section 3.2), solving the Euler equations for axisymmetric flow, has brought the 
sonic boom extrapolation technique up to the 90 ’s standard. 

Parallel computing is a fast growing subject in the field of computer science because of the 
promising speed in number crunching computations. A new optimizer (Section 3.3) for 
parallel computing concept has been developed and tested for aerodynamic drag minimi- 
zation. This optimimizer is also coupled with a parallel CFD code so the whole optimiza- 
tion process is parallel. This is a promising method for CFD optimization making use of 
the computational resources of workstations, which unlike supercomputers spend most of 
their time idle. 

Finally, the OAW concept is so attractive because of its overall performance in theory. In 
order to fully understand the concept, a wind-tunnel model is built. CFD Navier-Stokes 
calculations helps to identify the problem of the flow separation (Section 2.4), and also 
help to design the wing deflection for roll trim and alleviating the flow separation. 
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FORTRAN Program 
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LHF.f 
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Dr. Samson Cheung {Tel : <4 

MCAT rnStitdtS 

NASA Junes Re search Center 
tyS 25#~$i 

Moffett Field, CA 94039 

bite if 93/3/4 :::: : version 3 ;i 


ITel : < 4 11) -S 04 -44 62 J 


PARAMETER ( KMAX- 220 , LMAX- 1 , JMAX -3 51 , NMAX- 9 00 ) 

DIMENSION X ( KMAX , WAX , 3 MAX ) , Y < KKAX , ti 4 AI , JMAX ) , Z ( DtAI, tMAX , JMAX } 
REAL R(NMAX), S(NMAX) , SPfNKAX), TAU(NNAX), FTAU < NMAX ) 

COMMON/P AF/ FKACH,PFAC 
LOGICAL TODY 
TODY - . TALSE . 


OPEN{UNIT«3 , FILE** Lhf.in’, STATUS- ' OLD’ ) 

Read tie ' 'input ' parameters ‘ 

NAMELIST /PARA/ FKACH , PFAC , RO , R1 , IN ARE A , LIFT, TORX 
READ{ 3 , PARA) 

WRITEf 6 /PARA) 


Input free- stream Mach number » FMACB 

If TORX > 0 , sonic boom verses time, else verses distsRca 


Is the surface grid contains the whole configuration, or 
only half -plane or only qua ter- plane ? 

PFAC * 1 . \ Whole plane ' 

PFAC 2.,. ; ? belf-piene 

' PTAC - 4 , ' | Quater-plaoe 


‘ R 1 >111 ' be. the diatance yherW the signaturefe c^tnrSd^ " 


If ^ read in area distribution, I MARIA *• 0 


read the grid " 
the wlng-nody case “ : £ ; :| 
read in function 
read the R-f unction 
read a Signiture at ' R0 


I MARIA * 1 
1 MARIA 
IN AREA - 2 

R 0 > 0.0 


PI - 4 . *ATAN( 1 . ) 
JDIM - JMAX 
JDIM * 301 


If BO > 6 , We read the pressure sigsuture at RO and extrapolate 
IF(R 0 .GT. 0 .) GOTO 7 90 ""' v 

Find the area distribution of body configuration ' (sample ' Case) -• 
IF ( IN ARE A . LT . 6 ) THEN 
CALL WBODY ( JDIK, S , TAU) 

mn - 

CALL CONE ( JDiK, $ , TAU) 

CALL WING{ JDIM, S t TAU) 

CALL SEARS ( JD 1 M/S ,TAU> 

; . CALL BULLET! JDIK, S,TAU) ' 

GOTO 270 

ENDTF 


Read in the given area distribution 
ir! INAREA. EQ. 0 ) THEN 
OPEN( UNIT- 2 , FILE- 'area. in' ) 

DO 50 J-l , NMAX 
READ( 2 ,*, END- 75 ) TAU<J),S(J) 

S(J) - S(J)*PFAC 
CONTINUE 
CONTINUE 
CLOSE( 2 ) 

JDIM - J-I 

OPEN( UNIT- 2 , riLE-'area. io‘ ) 

do io j-i , jdim 

WRITE!!,*) TAU(J),S(J) 

CONTINUE 
GOTO 270 
ENDIF 

Read in the F- function or define a F- function by calling FUNC 
and integrate out the equivalent area by calling EAREA 
IF< INAREA. EQ. 2 ) THEN 

CALL FUNC ! TAU , FTAU , JDtM) 

OPEN! UNIT- 4 , FILE-'f .det' ) 

DO 100 J-l ,NMAX 

READ! 4 , * ,END=110 ) TAU ( J ), FTAU ! J ) 

CONTINUE 

CONTINUE 





SOURCE TEXT 


JDIM - J-l : 

CALL DISTARC ( TAU , FTAU , J-l ,TAU , FTAU , JDIM, 10 . , 0 ) 
CALL EAREA ( S , FTAU , TAU , JDIM) 

GOTO 270 
ENDIF 

Read the PLOflD surface grid "file (Planar fOXlUt) 
and find tbe equivalent »r« distribution •; 

OPEN (UNIT- 11 , FILE* ’grid, in’ , FORM- 'UN FORMATTED* ) 
READ(ll) KDIM,LDIM, JDIM 
DO 200 J*1 , JDIM ' 

READ(ll) ( (X(K,L, J) , K«1,IDIH) ,L-1,LDIH), 

4 ((Y(X,L,J), K-1,KDIK) ,L-i,LDIH>, 

4 ( ( Z ( K , L, J ) , X-l ,XDIH) i L-l ,LDIM) 

CONTINUE 

CLOSE (II) 

CALL EQUAREA( EDI M,LDI M , JD I H , X , Y ,1 , KMAX , LKAX , JMAX f S| 


DO 220 J-l ; JDIM 
S(J) - PFAC*S ( J J 
TAU( J )*X ( 1 , 1 , J ) 
CONTINUE 


c Obtain the a*ti»«ti?* of "tb* iWntrUwtion 
IF(LIFT.EQ. 1) CALL BFUNC( JDIM, S , TAU) 

jDikMi i 

DO 300 J-2,JDIKM1 

Al- ( TAU( j )-TAU( J+l ) )/( (TAU( J )-TAU( J-l) ) * (TAU( J+l ) -TAU( J-l ) ) ) 
A2-(TAU( J)-TAU(J-l) )/( (TAU (J+l )-TAU( J) )* ( TAU ( J+l ) -TAU (J-l ) ) ) 
SP(J) - Al* ( S ( J-l )-S( J) )+A2* (S(J+1)-S(J)) 

300 CONTINUE 

A3 - ( TAU( 3 ) -TAU(1 ) )/( ( TAU (2 )-TAU( 1 ) ) * (TAU (3 )-TAU<2 ) ) ) 

A3 - (TAU( 2 )~TAU(1) )/( ( TAU(3)-TAU( 2 ) ) * (TAU (3 )-TAU(l) ) ) 

Al - A2-A3 

SP(1) - -A1*S ( 1 )+A2 *S(2 )-A3*S( 3 ) 

Al - (TAU (JDIM) -TAU (JDIM- 2) ) 

4 / ( ( TAU( JDIM) -TAU (JDIM-1 ) ) * (TAU( JDIM-1 )-TAU( JDIM-2 ) ) ) 

A3 • (TAU (JOIN) -TAU (JDIM-1)) 

4 / ((TAU( JDIM-1 )-TAU( JDIM-2 ))*(TAU(JDIM)-TAU( JDIM-2))) 

A0 - Al - A2 

SP(JDIM) - A2*S( JDIN-2 ) -Al*3( JDIM-1 )+A0«S( JDIM) 

0 1st order Sl>(l> - ( 5(3 ) -$( 1 ) )/( TAV( *>-***( 1 ) ) 

0 1st Ota** SP(JOIM) - ( S( JDIM)-S( JDIMM1) )/(m( JDTKJ-TA0<JDIIBC1 ) ) 


c Redistribute the * In equal spacing 

CALL Dt STARC ( TAU , S , JDIM , TAU , S , JDIM ,10. ,0) 

CALL DI STARC ( TAU , SP , JDIM, TAU , SP , JDIM, 1 0 . , O) 

DO 340 J-l, JDIM 

R(J) - SQRT ( S ( J J/PI ) 

340 CONTINUE 

OPEN ( UNIT- 1 2 ri LE- ‘ area . out' ) 

WRITE( 12,400) 

400 FORKAT( 37H|Thie ia equivalent area distribution) 

DO 450 J-l, JDIM 

WRITE( 12,580) TAU ( J ) , S( J ) 

450 CONTINUE 

c wsi*E( 13 , * & y-'-'- ■ : : : : : : * : 

o WRITE(ia,451) 

o 451 FORMAT (ifHlThi s U tie derivative of the itm distribution) 

C Do 4 55 J-i ; JDtH ' 

c WRIT*(12,5i«j fAtr<dj f &m 

c 455 CONTINUE 

CLOSE ( 12 ) 


The intarference betvaen ving and body, aeiple cane only. 
IF(WBDY) CALL WB ( JDIM, SP, TAU) 


Obtain tbe tighthili F- function 
CALL LIGHT1 < TAU , R,SP, JDIM , FMACH , FTAU ) 


C Nr lie out iS» Lighfcbill r~ function at tie body surf nee. 

OPEN (UNIT-1 3, FILE-' fin. out* ) 

WRITE (13,554) 

556 FORMAT (498 I This ia the Ligbthill r-function at body surface) 
DO 560 J-l, JDIM 

WRITE(13, 580) TAU ( J ) , FTAD( J) 

560 CONTINUE 

580 FORMAT ( 2X,E16.8,1X,E16. 8} 

CLOSE (13) 

790 CONTINUE 


Obtained the pressure signature at distance ftl fro* the body 
CALL FFN( FTAU , TAU, FWtCH , JDIM, R0 , Rl , TORX) 

CLOSE(3 j 

STOP 

END 






U N IX® 

/ORTRAN Progrc 


f 


SOURCE PROCRAM 

LHF.f 


LINE# 


95 


pi** 4 . •atXn( i.) 

BETA“SQRT( FMACH* *2-1 01 
TAU{1)=0. ' 1 

FTAU(l)-0. 

DO 95 N-i,N 
FTAU ( M ) * 0 . 


DO 100 J-1, N 
DO 102 1-2, N 

IF(ABS{R(I)) THEN 

* l.E+10 
F4 - 0 . 

GOTO 98 
ENDir 

*B-2.0/(BETA*R(I ) 3 
ABl-ABS(AB) 

Fl-SQRT(ABl) 

F2-SP( I ) -SP ( 1~1 ) 

F3-F1*F2 
F4-F3/( 2 . 0*PI ) 

Jwi-f ' V » ' TAU < 1 ) )/( bbta* R ( r ,, 
IF (Zl.LT.XLO) GO TO 96 
ll <JJ- LT - 4 -0) GO TO 9 7 

BZlio 1 GE 4 60 10 * 8 

FTAU( J)-FTA0( J)+HZ1*F4 
GO TO 99 

HZ1*-02937*Z1*Z1- .2175*714 75 H 
FTAU(J)-FTA0( J)+HZ1*F4 
GO TO 99 
BB» 1 . 0/ ( 2 . 0 • Z 1 1 
BB1-ABS ( BB ) 

HZl-SQRT(BB) 

» ccS>™ J) — 

102 CONTINUE 
100 CONTINUE 
RETURN 


96 


97 


98 


DATE 


7/14/94 pace# 


TIME 


5:00:1 1 pm 


SUBROUTINE EQUAFEA ( 1DIM J W* , X , " 

till* subXouti»« fl®«® 

tjopisodial rule* 

D x«E« S io^ i Vt^ J - JMXX ’ 

REAL S(JOIM) 

tx> 10 j-1 , jdIM 

darea - o. 

M)D - 

DAFEA * DAFEA * 0.5*M)D‘B 
CONTINUE .,-,,...,.,; 

Tfc* UUVMtlKi fc**® tnTM 71 

n « vi i LDIM » J 1 “ Y(KDIM,LDIM» J) 

M)D - < i(l > LDlM,J)+Z(* DlM ' IjDlH ' J> 

base - ABS<0.5*ADD*H) 

; Tb. .r.. s«o«^ 

$(J) - ABS( DAREA) -BASE 

continue 

return 

END 




| LINE » SOURCE TEXT 


DATE 

7/14/94 

TIME 

5:00:11 pm 





SUBROUTINE FFN ( F , X , FKACH ,NP , RO , R1 ,TORX) 

c - V .V: •• :v : - 

c THU program uses r-functioo theory to predict the pressure 

C signature at fax field when an initial pressure signature is given 

C • " =! ■' • M. " -7 

PARAMETER (NMAX-1400) 

DIMENSION F (NP ) , X(NP) , Y (NMAX ) , P (NMAX ) 

DIMENSION YS TP (NMAX) 

DIMENSION DBLVL( 3 ) 

OPEN(UNIT-14, FILE-' p. out' ) 

C " " ! 

C Input of initial parameters and pressure signal 

C 

c TKACH - Free~*tresi* ffaoh number 
C RO ' * initial distance from the body (altitude) 

C R1 » Final distance front the body {altitude)' 

C Nf>’ "'i Nuaiber of data (NP < NMAX) ' 

C NMAX - this should be large enough to resolve the signatur* 

C Tbux > 0/ sonic iwdi Pereep tine, else verses distance 

NAMELIST /XPSCALE/ XSCALE, PSCALE , Ag , Pg , Pa, IRISE, ALT 

C Read the input par am# tax* 

READ ( 3 , IP SCALE ) 

WRITE ( 6 , X PSCALE) 

C " Define the jpii'i(dBe^ ^function theo.ry 

GAMMA - 1.4 

B - 5QRT(FMACH**2 - 1.) 

CAP - (GANMA+1 . )*FMACH**4/(SQRT(2 . )*B**1 . 5) 

SRRI - SQRT < SI ) 

C If RO > 0, ex trap olate from RO to il. First calculate the F-fS, 
IF(RO'.GT.O. ) THEN 

OPEN ( UNIT- B , FILE- 1 p. rO ' ) 

DO 15 I"1 , NMAX 

READ( 8 , * , END— 30 )X( I ) ,P ( I ) 

15 CONTINUE 

30 CONTINUE 

CLOSE ( 8 ) 

NP - 1-1 

5RR0 - SQRT(RO) 

DO 35 1-1 , NP 

F( I ) - SQRT( 2 . ‘B^RO) *P ( I )/ (GAMMA* FKACH* FKACH) 

1(1) - X(I) - B*I0 4 CAP*SRF0*F(I ) 

3 5 CONTINUE 

ENDIF 

C " ' ; ■ : : : S? • : M M P: > - • ; M% V: . ••• : ; ' ' ■ ■ 

C Y is transposed COdrdioate 7 : 7! 

C" •• v- '• : • , . 

c 

t>0 50 I-1,NP 

Y{ I ) - X(I) - CAP*SRR1*F( I ) 

C N*ITI<U,«) X{X),f<I) . • 

c 49 FORMAT ( 4 3R fthle U • transposed f-fpn<?ti«i et 'surface) 

SO CONTINUE 

C Find the largest and smallest esTuas of Y 

C . : ; : : • , ' : -;.V ' : i : i. ■' V V/V.-Ui- ^V:‘ : ^ : f!--:- ■ * ' " L : ^ - 

YMAX - -I.fc+B 
YMIN - l.E+B 
DO 55 1-1 ,NP 

IF(Y<I) .GE. YMAX) YMAX-Y<I) 

IF(Y(I) . LE. YMIN) YKIN-Y(I) 

55 CONTINUE 

C : '• > : '• ' 'x ; V ' \ : : : f*, / ■ - : ' 

c Prist out the Integral curt* of the shifted r- function 
c '••■•■•■•qat iHTr(OT,r,y) - m&W&m* 

C Need to ■arch in T-Oifactiot, define the step 

C 

YSTP(l) - YMTN 

YDIS - YMAX- YMIN 

DY - YDI S/FLO AT ( NMAX- 1) 

DO 80 J-2 , NMAX 

YSTP(J) - YSTP(J-l) + DY 
80 CONTINUE 
C 

C Perch through the shifted p«functionr Chech area-balance end 
C place the shock < 

CALL MARCH ( NMAX , NP, Y, YSTP , F) 

C obtain the solution 

C NOTE: It TW»>0, the sonic boas is in the for* (P~Pi*f ) Va ti*a 
C or it is In the tor* iP-Pinf )/Pi»f vs diatanos . 

DO 150 1-1 ,NP 

P<I) - GAMKA*FMACH*FMACH*F(I)/SQRT(2. *B*Rl) 

X<I) - Y( I ) + -B*R 1 
150 CONTINUE 

C KaXe the data points in evenly distributed »atot>er and 
C scale the sonic if desired 

DO 180 I«1,NP 

X(I) - X(I)*XSCALE 
P(I) - P( I ) UPSCALE 
180 CONTINUE 

C , Atnospher ic aspect;; _ 

C ALT • Altitude 

C Ag « speed of sound at ground-in ft/sec 

C P0 - reference pressure ib/ft'*! * SQRT{Pa*Pg) 

c , pg * pressure et the ground 

C Pa * pressure at flight altitude 

P0 - SQRT(Pg*Pa) 

VEL - FMACH * Ag 
TREF * X(1)/VEL 
I F ( TORX . GT - O . ) THEN 
DO 260 I-1,NP 

X(I) - X ( I J/VEL - TREF 
P(I) * P(I)*P0 
260 CONTINUE 


The signal (DP vs Ti»e) is calculated, use a enpirical progra* to 
calculate the rise tise, and «abed the rise tiee into the signature. 


C 








SOURCE PROGRAM 

LHF.f 


7 / 14/94 PACE * 



LINE # 

...41 

4 

4 

"‘41 

1 

'41 

41 

7 

4! 

if 


SOURCE TEXT 


Notes Unit need is still the stupid English unifct 
CALL RISETIME ( FMACH , P , X ,NP , ALT, I RISE ) 

Obtain 'tie noise level 
CALL NOISE (DBLVL,X,P,NP> 

Write the dft{PL) value Out 

WRITE ( 14 , 500 )DBLVL( 1 ) , DBLVL( 2 ) , DBLVL( 3 ) 

FORMAT 

l ( ’INoise level * , FI 0 . 4 , 1 PLdB 1 , 3X,F10 . 4 , ’ dB< A) 1 , 3X , F10 . 4 , ■ dB{C ) * ) 
END IF 


Write the sonic boon 

WRITE! 14 , 555 jIr! 

FORMAT! 2 8H 4 The pressure signal at Rl- ,F10.4) 
DO 670 1-1 ,NP 

NRITE( 14 , 700 ) X(I),P{I) 

CONTINUE 

FORMAT! 3X , E20. 8 ,2X,E16 -6) 

CLOSE (14) 







SOURCE PROGRAM 



DATE 7/14/94 PACE# 
T,ME 5:00:11 pm 


LINE 


438 

439 

c 

440 

44 1 

442 

c 

443 


444 

44$ 

446 


447 

448 

449 

100 

456 

120 

451 

130 

■'"' 457 ' 

c 

453 

454 



SOURCE TEXT 


SUBROUTINE INTF(NP,Y,F) 

Thin progra* print out the Integra^ carve of the shifted E-function 
DIMENSION F(NP) ,Y(NP) 

OPEN(UNIT«34, FILE= ' icuxve_F' , FORM- ’ FORMATTED • ) 

SUMF » 0 
WRITE (34,120) 

DO 100 J * 2 , NP 

DY - Y( J )-Y ( J-l ) 

SUMF - SUMF + 0. 5*DY*<F< J)+F( J-l} ) 

WRITE <34, 130 )Y(J) ,SUMF 
CONTINUE 

FORMAT ( 4 2H# Integral curve of the shifted F- function) 

FORMAT (2E1 6.6) 

CLOSE ( 34 ) 

RETURN 

END 







SOURCE PROGRAM 


LHF.f 

time 5:00:1 1 pm 



SOURCE TEXT 


SUBROUTINE SBXPT(NHAX,NP , Y, YSTP , F, INDEX , FS ,YS , III ) 

DIMENSION F(NP),Y(NP) 

DIMENSION YSTP(NMAX) 

DIMENSION INDEX (40) , FS ( 4 0 ) 

COMMON/SHOCK/ INSCT 

YEND - Y (NP ) 

FIRST - 1. 

DO 500 J-2 ,NHAX 
YS - YSTP(J) 

'Ciit tie point* on tW <rtrt:Ve for integration, etart searching from 

is to ;.i* : 

CXLt POINT ( NR, Y,F, INDEX, FS,YS,INSCt) 

After obtain the integration point*, m can integrate and 
find the *r*a ' ; 

iri iNSCT GTVi ) THEN ' ' 

IF( III . EQ . 3 ) INSCT - 3 
IS-INDEX(l) 

IE- INDEX (INSCT) 

CALL AREA(NP ,Y,F,TS,FS,IS,IE, IFLAT2 ) 

ELSE 

The tail shoe*. la already f creed, leave program 
IF ( INSCT. LE .1 - AND . YS . Gt . YEND* 1 . 05) RETURN 
FIRST-1. 

GOTO 500 
ENDIF 


tt( FIRST .GT . 0 . ) irtXTl - IFLAT2 
IF( IFLAT2 . EQ. 0) RETURN 
FIRST - -1. 

IT I FLAT » 0, IS la the point tlUi liiipO area balanced. 

If IFLAT2 and I FLAT! are in different *ign , i . e . , 

the correct point should be betvega* I and 1-1. 

tjae bisection eethod to fin® th* <wrect point Y(Istart) 

I F ( i FtATi * t FLAT2 LT . OO THEN 
Y1 - TSTP(J-l) 

Y2 - YSTP(J) 

NC - 500 
DO 200 IC-1,NC 
YS - 0 . 5* ( Y2+Y1 ) 

CALL POINT(NP,Y,F f INDEX, FS,YS, INSCT) 

XF(II2.CQ.3) INSCT - 3 
IS-INDEX(l) 

IE-INDEX{ INSCT) 

CALL AREA(NP ,Y,F,YS,FS , IS,IE»IFIAT0) 

IF(IFLATO.EQ-O) RETURN 
IF( IFLAT0*IFLAT1 .LT. 0) THEN 
Y2 - YS 

IFLAT2 - IFLAT0 
ELSE 

Y1 - YS 

IFLAT1 » IFLAT0 
ENDIF 
CONTINUE 

WRITE ( * , * ) ‘After ' #NC, ‘ steps of bisection* 

RETURN 

ELSE 

IFLAT1 - IFLAT2 
GOTO 500 
ENDIF 


FIRST - 
CONTINUE 
RETURN 
END 






SOURCE PROCRAM 

LHF.f 


7 / 14/94 


sssnss'riswisi* • r ■ inde]i - ™ scr, — XWCT 

DIMENSION INDEX (40) ,FS(40) 

- « o f po lnts being *■*»*»,• it iWw point to <Jo integration 

INSCT - 6 "'"■■■■ • . 

IF(YS LT. Y ( I j ) THEN 
INSCT - INSCT + 1 
FS( INSCT) - 0. 

INDEX (INSCT) - I 
ENDIF 

DO 100 i-2,NP 

FACI - rs - T(I) 

FAC2 - TS - T(I-l) 

IF( FACI *FAC2 .LE. 0.) THEN 

IF(ABS(Y(I)-Y(I-1)).LE.1.E-I4) THEN 
"*£*•< * p * ) * 2EROOOOOO ,YS , Y( r ) ,t* 

RETURN 

ENDIF 

INSCT =» INSCT + I 

SL - (F(I)-F(I-l))/(Y{I)-Y ( l- ln 

FS(INSCT) - F(I~1)4SIj*{YS-Y(I-1) ] 

INDEX( INSCT) - I ' Mi 1)) 

ENDIF 

CONTINUE 

ir<YS .GT. Y(NP) ) THEN ‘ * * ' : * ‘ 

INSCT - INSCT + 1 
FS( INSCT) - F(NP ) 

Y(NP) - YS 
INDEX ( INSCT ) - NP 
ENDIF 





SOURCE PROGRAM 

DATE 7/14/94 

LHF.f 

time 5:00:1 1 pm 



SOURCE TEXT 


SUBROUTINE WBODY ( JDIM, S , TAU ) 

This subrouti ne . f lad the area d Isfcr ibufc Ion of 
the wing-body configuration . 

DIMENSION S( JDIM), TAU< JDIM) 


PI«4.*ATAN(1. ) 

ANG-21.*PI/180. 

ANG^-BO. *PI/180. 

DX-2 5 . 52/FLOAT ( JDIM-1 ) 

MD - 6 . ' ■■■■■■' ' -■■■ 

TAU<1) - 0. 

DO 2 J-2,JDIM 

TAU(J) - TAU( J-l )4DX 
TTT - TAU(J )-7 .01 
IF ( TTT . GT . 0 . ) TTT“0. 

RR-0. 54-0. 011*TTT**2 
S< J) - PI *FR*RR 

IF(TAU< J) .GT. 8,21 .AND. TAU( J ) . LT. 12 . 25) THEN 
AA * 4 . *0 . 5*0. 05*TAN{ ANG ) * ( TAU{ J )-8.21)**2 
S<J) - S(J) 4 AA 
ENDIF 

IF<TAU<J) .GT.12.25 .AND TAU ( J ). LT . 1 5 . 77688849 ) THEN 
B2 - 0 . 05* < 16 . 29-TAU ( J ) } 

B2 - 2. 91*<<TAU<J)-12.25)/(15. 77688849-12. 25)) 

B1 » (TAU{ J)-8 ■ 21 ) *TAN{ ANGJ-B2 
HI « 0. 05*B1/TAN( ANG) 

AA - 4. *<0.5*B1*H140.5*(H1+H2)*B2) 

S(J) - S<J) 4 AA 


B2 - 2. 91*<<TAU<J)-12.25)/(15. 77688849-12. 25)) 

B1 » (TAU{ J)-8 ■ 21 ) *TAN{ ANGJ-B2 
HI « 0. 05*B1/TAN( ANG) 

AA - 4. *<0.5*B1*H140.5*(H14H2)*B2) 

S(J) - S<J) 4 AA 
ENDIF 

IF(TAU(J) .GT. 15. 77688849 .AND. TAU< J ) . LT. 16 . 29 ) THEN 
AA - 4 . *(0.5*0. 05*TAN ( ANG1 ) *( 16 . 2 9-TAU( J ) ) * *2 ) 

S(J) - S(J) 4 AA 
ENDIF 

IF (TAU ( J ) . GT . 18 . 93 .AND. TAU( J ) . LT. 17 . 52 ) THEN 
SLOP*<0. 15-0. 54)/<17. 93-17. 52) 

RKR-0 . 54+SLOP* (TAU( J )-17 . 52 ) 

S(J) - PI *RRR* *2 
ENDIF 

IF(TAU(J) .GT.17, 93) S(J)«PI *0.15*0. 15 

CONTINUE 

RETURN 

END 




SOURCE TEXT 


SUBROUTINE CONE ( JDIM, S ,TAU ) 

This subroutine find the are* distribution of the oone- cylinder 
with half-angle 3. 34 d*gr*e and 8.6 units of length. 

DIMENSION S(JDIM) f TAU< JDIM) 

PI-4 . * ATAN { 1 . ) 

ANG - 3 . 24*PI/180 . 

DX - 1 6 . /FLOAT ( JDIM- 1 ) 

S(l) - 0. 

TAU ( 1 ) = 0. 

DO 2 J*2 , JDIM 

TAU<J)-TAU(J-1)+DX 
I r { TAU { J ) . LE . 8 . 6 } THEN 
R - TAU { J ) * TAN { ANG } 

ELSE 

R - 8.6* TAN { ANG ) 

END IF 

S(J) « PI*R*R 
CONTINUE 
RETURN 
END 






SOURCE PROGRAM 

LHF.f 


SUBROUTINE SEARS < J DIM , S, TAU) “ — - — 

VHi tlo^l "** ■«•««»*», of M* i«r S -H,,ok >odj 

piWW^ON S < JDIM) , TAU^DIK) 

JMV* ' ° £ P ° tnt 

JSNG - JDIM - JBDY 


FI-4. *ATAN{1 . ) - 

BL - 1. 

TOTL - 1.9* BL 
F - 23.5 

D THETA - PI/FLOAT ( JBDY-1 ) 

DX - BL/FLOAT( JBDY-l ) 

RKAX - BL/j 2 . *E) 

S<1) - 0. 

TAU(l) - 0. 

Owatant* ot $«ar*~A4ajaa body 

vr lte( * , * ) • I nput Abaae/A»ax - 

r«ad ( * , • ) AR 

wrlt«(* ,*)' input XMAX' 

raad(*,*) AA 

XMAX - AA*BL 

CONST - AR/PI 

\ * < ? : - _ 1 . ) ) 

’ DO i J-2 /JBDY 

. tfftA - P1-0TRBTA'* FXOAY* J~1 ) 

) " U ^OS(tHETA) ) *Bt/2 , 
tAB(J) - FL6AT(J-1)*DX 

THETA - ACOS(2. *TAU(J1/BL - 1 i 

:£eaxa-&a.aok body :::: . 

• JW* ” (5t#(*HEtA) ) * *3 
lUack- Ada** Body 

POS - CONST* ( PI-THETA+0 . 5*SIN{ 2 . *THETA) -f 
(4 ./3 . ) *C1* (SIN< THETA) )* *3 ) 

R - RMAX*5QRT( ABS/POS > ) 

S(J) - PI*r*r 
CONTINUE 

Add a ating- 

DX - (TOTL-BL)/FLOAT(JSNG) 

DO 5 J-JBDY4l,JDIM 

TAU(J) - TAU(J-l) + DX 
S(J) - S(J-l) 

CONTINUE 

RETURN 

END 


^onix™ r i: 

Fortran Program 


SOURCE PROGRAM 

LHF.f 


7 / 14/94 


SOURCE TEXT 


H - a . 

DIMENSION S(JDIM),TAU<JDIK> 

PI-4. * ATAN ( 1 • ) 

Bt - 4. 

GAMA - 0.65 * 

RBASE - 0.25 

A - RBASE/ (BL**G AHA) 

TOTLEN - BL 4 2.*BL 

DI - TOTLEN/FLOAT { J DI M- 1 ) 

S(l) “ 0. 

TAU ( I ) - o. 

DO 2 J-2 , JDIM 

TAU{J)-TAU(J-1)+DX 
IF(TAU(J) GE.BL) THEN 
R - RBASE 
ELSE 

R - A*TAU(J)**GAKA 
ENDIF 

S( J ) - FI*F*R 
CONTINUE 
RETURN 
END 


SOURCE PROGRAM 


DATE 



LHF.f 




gag 2 


7 / 14/94 
T,ME 5:00:1 1 pm 


PACE # 


15 


LINE # 


r 42 | 
'■431 c 
'44 | 
'4S c 
W\ 

Wt 

>w\ 

*ht 
'53 j 
'54 
'55 
'56 
'57 
'56 
'55 



SOURCE TEXT 


SUBROUTINE WING ( J DIM, 5,TAU) 

This subroutine find the Area distribution of the low^aspeot-ratio wing 
DIMENSION S(JDIM) ,TAU<JDIM) 

PI= 4 . * ATAN(1 . ) 

STING ** 0. 

DX - 3./FLOAT(JDIM-l) 

S<1) - 0. 

TAU(l) - 0. 

DO 2 J«2,JDIM 

TAU(J)-TAU(J-1>+DX 
I F ( TAU ( J ) . LT . 2 . ) THEN 

Z - (PI/12 . 5 ) • ( TAU( J ) -0 . 5*TAU( J ) *TAU( J ) } 
s ( J) - Z 

IF(TAU(J> .GT. 1. 70897) THEN 
5TING-PI*0. 062 5*0 .0625 
S(J) - Z + STING - Z*0 . 125 
ENDIF 
ELSE 

STING-PI*0. 0625*0.0625 
S(J) - STING 
ENDIF 
CONTINUE 
RETURN 
END 











U N IX™ _ . 

FORTRAN Program 


,L.''.«ife ; .v 


SOURCE PROGRAM 

LHF.f 


DATE 


TIME 


i MM 




LINE # 


SOURCE TEXT 


SUBROUTINE BFUNC ( JDIM , S , TAU ) 

This subroutine obtains the B-fun'otion fr-qn fort... VP. and add it 
''into the equalivaqt frtea. 

PARAMETER (NMAX-800) 

DIMENSION S(JDIM) ( TAU(JDIM) ,B<NMAX) ,X<NKAX) 

COMMON/PAR/ FMACH , PFAC 


OPEN{ UNIT-33 , FILE* ' coef . dat ' , FORM- 'FORMATTED ’ ) 

READ( 33,12) 

READ{ 33,12} 

READ( 33,12) 

READ{ 33,12) 

READ( 33,12) 

DO 10 I«1,NMAX 

READ (13,15, END- 17) XfX) ,CL,CD,SLOD,B(I) ,CM 
READ( 33,*, END-17) X(IJ,B(I> 

B { I ) - B ( I ) *PFAC 
CONTINUE 
CONTINUE 
FORMAT(IX) 

FORMAT { 6E13 . 5) 

CONTINUE 
CLOSE (33) 

NPOINT - I-I 

OPEN ( UNI T- 3 3 , FILE- ' bf o . dat ' ) 

DO 20 1-1, NPOINT 
WRITER 33 , * ) X ( I ) , B( I ) 

CONTINUE 

I START- 1 
DO 50 J-l ,JDIM 

DO 30 I- I START, NPOINT 

XF{ AB5 ( X( I )-TAU( J )) .LE.1.E-I0) THEN 
5(J)»S(J)+B(I> 

I START- 1 
GOTO 40 
ENDIT 

IF{X(I) .GT.TAU< J) ) THEN 
IF(I-EQ-l) THEN 
BF-0. 

XF-0. 

ELSE 

BF-B(I-l) 

XF-X(I-I) 

ENDIF 

SLOPE- <B( I )-BF)/ ( X(I )-XF) 

BT - B(I) + SLOPE* (TAU(J)-X(I ) > 

S(J) - S<J) + BT 
ISTART-I-1 
IF(I.EQ.l) I START- 1 
GOTO 40 
ELSE 

IF< I. LT. NPOINT} GOTO 30 
S(J) -S(J> + B{NPOINT) 

ISTAFT-I 
GOTO 40 
ENDIT 
CONTINUE 
CONTINUE 
CONTINUE 
RETURN 
END 


7 / 14/94 

PACE # 


16 









SOURCE PROGRAM 


DATE 


PACE# 




SUBROUTINE WB( JDIK, SP , TAU ) 

This subroutine obtains the wing-body Inter f erenow correction 
and add it Into the derivative of equalivant area 
This is a teat for wing-body ca*« 

DIMENSION SP(JDIM) ,TAU( JDIM) 

COMMON/PAR/ FMACH, PFAC 

do id 

IF(TAU(J) .GE.8.21 .AND. TAU< J ) . LE . 12 . 2 5 ) 
t SP<J)-SP<J>+4.*.05*. 54 

IF{TAU(J) .GT.12.25 .AND. TAU(J) . LE. 16 . 29 ) 
l SP(J)-SP(J)-4. V 05*. 54 
CONTINUE 
RETURN 
END 




LINE # 


SOURCE TEXT 


846 

847 

848 
84 $ 

850 

85 1 

852 

853 

854 
85$ 
858 
85? 
858 
85$ 
866 
861 
862 

863 

864 
86 $ 
866 

867 

868 
86$ 
876 

871 

872 

873 

874 

875 

876 

877 

878 
87 $ 
886 
88i 

6 

8 

8 


C 

c 6 


20 

75 

BO 


SUBROUTINE FUNC < TAU , FTAU , JDIM) 

DIMENSION TAU( JDIM) ,FTAU< JDIM) 

NAMELIST /FFUNC/ YF,ELAM f C , H , B, D, E , BL, YR , DEL 

■Read 'the Input parameters " 

RE AD (3, FFUNC) 

WRITER 6, FFUNC) 

TAU(1)*0. 

FTAU ( 1 ) * 0 . 

DY-YR/FLOAT( JDIM-1 ) 

DO 10 J«2,JDIM 

TAU<J)-TAU(J-1)+DY 
IF<YF.EQ.O. ) GOTO 6 

IF(TAU< J) .LE . YF/2 . ) FTAU ( J ) ■ 2 . *TAt?(J) *H/YF 
IF<TAU< J) .GE .YF/2 .0 .AND. TAD< J ) . LE . YF) 

4 FTAU< J)«C* (2*TAU(J)/YF-1 . } - H* < 2 V *TAU<J J/YF-2 . ) 

7r(TAU(7) CEYF .AND. TAD<J) . LEE LAN) 

4 rTAD(J)-8*(TAU(j)-YFJ+C • 

IF<TAU<J) GE .Yr .AND. TAU< J ) . LE DEL ) FTAU(J)-C 
IF<TAU( J ) . GE . DEL .AND. TAU< J) . LE .ELAM) 

4 FTAU{J)-BMTAU(J)-DEL)+C 

IF ( tAU ( J ) . GE . ELW(' * -"AiMb’.’' tMJ ( J ) . LE . BL ) 

4 FTAU<J)-B*<TAU( J)-ELAM)+(ILAM*B-D) 

IF(TAINJ) GE .BL) 

4 fTAU<J}» < ELAM*B“JHBM BL-ILAH) ) >/< BL-YR) 

4 FTAU{ J }--E/{ TAU( J )- (BL-ABS (BL-ELAM)/10- ) ) 

CONTINUE .v.v ■ .w\w . v .. ... .......... 

CALt D18TARC ( TAtJ FTAU , JDIM , TAU , PTXtf j 

WRITE{13,75) 

DO 20 J*1 r JDIM 

WRITE ( 13,80) TAU(J) iFTAU(J) 

CONTINUE 

FORMAT( 2 21 IF- function from input) 

FORMAT ( 2X,F8.4,1X,E16.8) 

RETURN 

END 


1 




SOURCE PROCRAM 

LHF.f 


SUBROUTINE EAREA( S , FTAU ,TAU , JDIM) 

DIMENSION S(JDIM) , TAU (JDIM) ,FTAU(JDIM) 

DIMENSION F( 900 ) 

are* fro* F- function via Abel Tranafon 


SOURCE TEXT 


A<XJ * < T ( t )/SQRT( y-fc ) dt ) dy 


s (1) - 0. 

TAU(I)-0. 

DO 10 J-2 , JDIM 
SS - 0. 

DO 7 1=1, J-l 
DY-TAU(l4l)-TAU(I) 

FINGRL » 0. 

DO 5 1=1 ,1-1 

DT - TAU ( K+l ) -TAU ( K ) 

FINGRL - FINGRL + DT*FTAU( I )/SQRT(TAU( I )-TAU( r ) ) 
CONTINUE 

SS - SS + 2.* FINGRL* DY 

+ 4 . * SORT ( t AU < a J -TAU ( I )j * F TAU <1 J*DY 

CONTINUE 
S(J) - SS 
CONTINUE 
TANN=S ( 3 )/TAU( 3 ) 

S < 2 ) -TANN * TAU ( 2 ) 

WRITE( 12 ,15) 

FORMAT ( 2 7H I Area froe given F- function) 

DO 20 J-l, JDIM 

WRITE( 12 , 80) TAU ( J J , S ( J ) 

CONTINUE 

FORMAT( 2X , F8 . 4 , IX , E16 . 8 ) 

RETURN 

END 



FORTRAN Program 


SOURCE PROGRAM 

LHF.f 


OATE 


7 / 14/94 


PACE# 


time 5:00:1 1 pm 


20 


•• SMi&i 








SOURCE TEXT 


SUBROUTINE RISETIM£( FMACH , P , T,NP , ALT, ^ * 
it. esDlxioal nethod to calculate the rise tine of a »o*ip boo» 
Rt.r?l»e derived from recession .MlyelJ of Air *£ * 

flight test data. Good for N-wa*e type of signal, nay b« «newbat 
' cbasefeatiee (shorter, rise tin®). 

' All unit' used, are English unit f J[ ! . . 

FMACH >• free^stxeae Mach n usher .;...■ 

P(T> » sonic boon ' 

*S» «* Shock strength* 

'■ At* - Altitude (ft) / 

>0 * Free- stxe a* pressure (lh/ft2) 

I rff - Rise tine (sec) 

"TEMP * Tenperature R*f4459.67«(#/5)F 
DIMENSION P(NP) ,T(NP) 

P0 - 2116. 2 
TEMP - 518-69 

I FOUNT - 0 
CONTINUE 

I FOUNT - I FOUNT 4 1 v --V-: ■-•-v; 

find cmt tie sfck strength 
psh - 0. 

ISH0 - 0. 

DO 30 1-1 ,NP 

I F ( T ( I } . EQ . T ( I 4 1 ) ) THEN 
IF(ISHO.EQ.O) XSB0-I 
PSH - ABS(P(t+i )-P(ISH0) ) 

ELSE 

ir(PSH.EQ. 0. ) THEN 
GOTO 30 
ELSE 

GOTO 40 
ENDIF 
ENDIF 
CONTINUE 
CONTINUE 
ISH - I 

I P ( PSH . EQ . 0 . ) RETURN . . - - ---- : : 

^^csiSii^S rise tine using Kir force data 
Y1 « 2. 92* FMACH - 7.3ft 

Y2 - YI 4 ((8.5* FMACH * * 2 - 45.9*FMACH 4 62.4))**(.5> 

AF1 - Y2 * 1000. 

AF - AF1 / (((ALT/1000.) - 5. >*2117.) 

VIS - 100. 4 0. 5* (TEMP-4 10. ) 

FT * (AF*VlS)*P0/< PSH* TEMP) 

FT - FT/1000. 

Nov calculate the rise tine assunning lp«t ha* 3»«*0 Else tine 
FT - 0 . 003 /PSH 

ENDIF . . . . . , . . ... . ... .^ . . . . . . Vt . ... . ...... .... .... .. . rw ^,.^.v.v.v,- - - 

WFITt(l4 / 80)RT . , _ . . ... ... 

FORMAT ( 3 8H| Rise tine (sec) of the bow shock la , rlO . 5) ....... .. 

originally, T(iSH0)-T<lSH) vith infinite •Abet a^hgtb, *6* uteat* a 
signal with the riae tine, between the index ISHO to Tfll.,.,,../-"----- 
Also extend the signal by the. anount of .ri*e....tM^f* . 

DRT - RT/rLOAT(ISH-ISH0) 

DO 200 I*ISH0 , ISH-1 
T(I41) - T< I ) 4 DRT 
0 CONTINUE 

DO 300 I*ISH+1 »NP 
T( I ) - T( I ) * RT 

0 CONTINUE x ; - 

i>( ikount.lt: 10 ) goto 12 

RETURN 

END 


SOURCE PROCRAM 


LHF.f 


-*1 -*1 fn n I >> V |i |jt' J. J. x t J, J. ^ w .. ' TEXT 

|iisTARC^x ,y ,h, xnew*y^e^ *NOTW^ piiLC ^ifut| * ^* * * ** * ^^'^ ^**^ *^* 



DIMENSION *(N).Y(N),»IEN(NNE«),YNEN(NNEH) 


bi^ P oH“ ”c £ !«lu l# n‘u« &*i 5 jr ( r b J outlJ,e «*«*’" 

*Jf •»* pot»t* of Ar^ttal^L 9 * 1 * 1 s ** cl “ 9 No ‘« «>»t 


«• « *?&* ssw&s^m *^ i “5 <** «.t 

m^rjsr^mim -v 2 ^** « i — «» .*,. 

o«put .„*? is &&& tj riiwjijj ; ii ?;. « m « 


PARAMETER < MAX* 2 000) 

DIMENSION S(MAX) , TOTARC (MAX) , XN(KAX) ,YN(MAX) 


' mirrv2^! b# L 0 * P 01 ® 141 aiw«a la MAX 

IF i2?iL LE * M OR - MAX LE NNEN) then 

WITEC,.) -sub DISTAFC : MAX Is t.. s thsn N or NNEN- 
END IF 


'. I*** lor toi*i arc length 

TOTARC (1) - 0 ^ :-::--::->>::-:::-:x:;.::.: : r : :; 

DO 10 K»2,N 

!s*sr‘ (i(f)-i(i-in*>2 , 

CONTI NOE 


obt * 10 ilP lilf I 

DE1T-FAC*(T0TA*C(N)/F10AT(NNEN-Il ) 

ELS£ LL DISTRI ^ DELT/NNE>, ' S ' ir,LA T) 

S(l) - 0. 

DO 25 K“2 ,NNEW 

ENDIF 


J-ai.txlbutloo, put u» « r ., i. « tu^r.ry .rroys w n 

YN(l)*y(i) 

XN(NNEW)-X(N) 

YN(NNEH)-Y(N) 


do «o j - 2,'nneW ' 

ARCNEW - S ( J ) * TOTARC ( N ) 

DO 55 K - 2, N 

IF( TOTARC ( X ) -EQ. ARCNEW) THEN 
XN(J) - X(X) 

TN{J) - Y(X) 

GOTO (0 
ENDIF 

IF(TOTARC(X) .GT. ARCNEW) THEN 
XI - X(K-1> 

12 - X(X) 

Y1 » Y(K-l) 

Y2 - Y(R) 

XX ■ XI + (X{K)-X(R-1) )« 

CALL ^ 2 ’ xi * J?T"' C < K > - -«> T A»C ( K - 1 ) , 

XN ( J ) - XX ' ' 

YN(J) - YY 
GOTO 60 
ENDIF 
CONTINUE 
CONTINUE 


S 1 ™ "**** lDt ° **• “tPUt XNBN, THEN 

XNEW(J) - XN( J ) 

YNEW(J) - YN( J) 

CONTINUE 

RETURN 

END 


I 


■UN 

FORTRAN Progrant: 


SOURCE PROCRAM 

LHF.f 


— DATE 

wm 3 ■% A tf\M 1 1 

/ / I j 

TIME 

5K>0:11 pm _[_ 




m - ^sa=iisa aaa — “ SOURCE TEXT 


SUBROUTINE blSTRI ( FANG , EFCS , S , I FINE ) 

PARAMETER (MAX-500) 

DIMENSION S(MAX) f DUH(HAX) 

. If iri5£l?*dUtrlbuttoo Is cuslterlog »t outer grid 

IF wllTl!rl*rSOB T DXSTRI = MAX U less then XrCS' 

STOP 

ENDIF 

IF(KFCS.EQ.l) THEN 
S<1) - 0. 

GOTO 40 

ENDir . • 


DZl - FANG 

KFH - FFCS-1 

DZETA - I * /FLOAT ( FFM ) 

CML™RKT(Dll , RRM, 0 . 0001 , 100 , RDBETA ) 
CALL FZ1< EFCS , HDBETA , DZETA , S ) ^ ^ 

IF (iFlNt.fiQl) TEEN 
DO 37 F-l t EFCS 

DUM(KFCS-K4l) - S<K) 

CONTINUE 
DO 38 K-1,*FCS 
S<K) - l.-DUN(I) 

CONTINUE 

ENDIF 

CONTINUE 

RETURN 

END 



SUBROUTINE FZ1 (LI ,TBETA, DET, Z) 


SOURCE TEXT 


DATE 

TIME 


7/14/94 
5-‘00:1 1 pm 


PACE# 


23 


• • ' ,U *. 

COMPUTES NORMAL! ZED NORMAL DISTANCE. z a , 

nTWrWCTftVI ' ™ , irx. 


DIMENSION Z( 250) 

IF(TBETA.EQ.l. ) then 
DO 10 L«I,L1 
Z(L)-0. 

10 CONTINUE 
ELSE 

DO 20 L-1,L1 
ETA-(L-1)*DET 
RR-(TBETA+1 . )/(TBETA-l . ) 

EEE-1 .-ETA 
EBETA-RR**EEE 

2o ^ii™ ETA_1 • 1‘iwwi/m*, . , 

END IF 
RETURN 
END 




SOURCE TEXT 


SUBROUTINE GRBET (DFM , NPTfFPCC , ICC ,BET ) 

N^TBpO . R » BET* ' 

MKICH GIVES DESIRED <*T AT THE "ALL 

DIMENSION Z ( 250) 

ICCL-ICC 
FPCCL-FPCC * DFM 
BETAl-BETA 
Zl-DFM 

det-i./npt 

BR-1. 

FF--Z1 

IICC-ICC/10 
DO 10 1-1,1 ICC 
BF-BETA1 „ . , 

BETA-0. 5* (BETA1+1 . ) 

CALI FZ1<2,BF,DET,Z) 
rr-Z(2)-zi 

ir<fr.GT.O. ) GO TO 15 
BETAl-2.*BETAl-l. 

10 CONTINUE 
15 CONTINUE 

DO 5 NIT-1, ICCL 
CALL FZ1(2,BETA,DET,Z) 
r-Z(2)'Zl 
ir<r.GT.O. ) THEN 

rr-r 

bf-beta 

ELSE 

FF-r 

bk-beta 

4 CONTINUE 

. ' CAlt TtU* # WWA/6BT # - • • ■ 

; ' V; ' 

BM1-BETA-1 - 
bfmi-bf-i . 

BRM1-BR-1 . 

RETURN 

END 



SOURCE PROGRAM 

DATE 7/14/94 

LHF.ff 

T,ME 5:00:11 pm 



SOURCE TEXT 


SUBROUTINE LININT( XI , X2 , Y1 , Y2 , XLOCAL, YLOCAL > 

Thi* subroutine linearly interpolate YLOCAL when given (ll r Yl > i <X2 , Y2 ) 
IF(X1.EQ.X2) THEN 
YLOCAL* ( Y2-Y1 )/2 . 

GOTO 100 
ENDIF 

SLOPE - (Y2-Y1 )/( X2-X1 ) 

YLOCAL - SLOPE* (XLOCAL-X2) + Y2 

CONTINUE 

RETURN 

END 




SOURCE PROCRAM 

LHF.f 


7/14/94 '’ACE# 


time 5 : 00:1 1 pm 



SOURCE TEXT 


SUBROUTINE MARCH ( NMAX, NP , Y, YSTP , F) 

DIMENSION F(NP),Y(NP) 

DIMENSION YSTP (NMAX) 

DIMENSION INDEX (40) , FS { 4 0 ) 

COMHON/SBOCK/ INSCT 

This subroutine *arches the Y direct ton am) 

chock. If the areas are balanced and then place the shock 


KOUNT - 0 
VEND - Y ( NP ) 

CONTINUE 
DO IND«1,40 
INDEX ( IND)*D, 

INSCT «= 0 
END DO 

CALL SHKPT < NMAX , NP , Y , YSTP , F , INDEX , FS , YS , 0 ) 

For tail ahocx* no seed to ©heck tha jjioaalhie poeitioi ot ahock 

IF(YS . GT 0.1* (YSTP (NMAX ) -YSTP ( 1 ) ) ) GOTO 400 


Only one possible location of shock 

if(insct.eq.3)'(^to;;4AO ;;;;;; 

More than one possible locations of shock 
IF( INSCT. GE . 5) THEN 
YSHK - YS 

CALL SHKPT < NMAX ,NP , Y , YSTP , T, INDEX , FS, YS , 3 J 
IF(YSHI.LT.YS) THEN 

'tN vino shock owbiewea '• . . 

CALL SHKPT < WKXX , NP, Y , YSTP , F , INDEX ,FS,YS,0) 

GOTO 400 
ELSE 

There 'are two separated shocks •' 

For the shock la actually locate on the turning edge of r~ function 
;; /if© heed to relocate. It 

Fix the Yl and Y2 of thia oeU region V 
I F(YS . CT . Y'< INDEX < INSCT j + 1 ) ) THEN 
Y2 - YS 
BIG - 0. 

DO ITEST-INDEX< INSCT) +1 ,NP 

IF<YS.GT.Y<ITEST) .AND. ABS(YS-Y ( ITEST) ) . GT.BIG) THEN 
Yl - Y(ITEST) 

BIG - ABS(YS-Y( ITEST)) 

ENDIF 

IF{Y ( ITEST ) .GE .YS) GOTO 300 
ENDDO 

CONTINUE 

Find Ys by biencting Yl a nd Y3 
NC - 500 
DO 320 IC«1,NC 
YS - 0.5* (Y2+Y1) 

CALL POINT (NP,Y,1\ INDEX, FS,YS, INSCT) 

DO II-INSCT,1 ,-l 

IF< INDEX( II ) -LE. ITEST) THEN 
INSCT - II 
GOTO 310 
ENDIF 
ENDDO 
CONTINUE 
IS-INDEX(l) 

IE-INDEX{ INSCT) 

CALL AREA (NP , Y , F, YS , TS , IS , IE f IFLAT0) 

IF(IFLAT0 . EQ.O ) GOTO 400 
IF< IFLATO .GT. 0) THEN 
Y2 - YS 
ELSE 

Yl - YS 
ENDIF 
CONTINUE 

WRITE(*,*) 'After 1 ,NC, ' steps of bisection' 

ENDIF 
GOTO 400 
ENDIF 
ENDIF 


Forii the shock 


IF( INSCT . LE . 1 .AND. YS.GE.YEND) RETURN 
FDIS - FS( INSCT )-F5( 1 ) 

if ( FLOAT ( INDEX { INSCT }- INDEX { 1 ) } . eq . 0 . 0) then 

WRITE ( 15,*) 1 SDT: ZERO DIVISION ABOUT TO HAPPEN In MARCH' 

<» - 1 . e32 
else 

DF - FDI S/FLOAT ( INDEX ( INSCT) -INDEX ( 1 ) ) 
end if 

F{ INDEX ( 1 ) ) - FS(1) 

Y ( INDEX { 1 ) ) - YS 

IS - INDEX ( 1 ) + 1 

DO *50 I - IS, INDEX (INSCT) 

Y { I ) * YS 
F{ I ) * F(I-l) + DF 
CONTINUE 

I F{ KOUNT . EQ. 20 ) THEN 
WRITE ( * , * ) ' KOUNT=20 1 I * 

RETURN 

KOUNT - KOUNT+1 
GOTO 100 
ENDIF 










Appendix B 


SAMGRID (Fortran Listing) 


Appendix B 


July 15, 1994 


16 





LINE# 


SOURCE TEXT 


PROGRAM SAMGRID 
Include ” sgr id. con" 

Df / Stoson Cheung 
Date: Dec. , 1993 

TMfi subroutine reads a surface grid 1 b Airfoils 
sections sad reforests it tic produce 
* huff ace grid of axisymeetrlc crOSs-section*. 

Date* Dec./ 1993 Version 3.0 

read input geoaw try 


OPEN { UNlf- 1 0 , f 1 Lfi- ' W&GRl b . IN » , STATUS- ‘ OLD ’ , FOftH- • tt>feMAfTED ' ) 
OPEN(UNiT-30, rrfcH> 'NACGRTD. IN’ STATUS- ‘OLD ' , PORK- • TORMATTED r ) 
OPEN (UNIT-40, FI LE- ' XZSI ZE . IN ' , STATUS- ' OLD ' , FORM-' FORMATTED '■ ) 


..MSIC I of sections ( strsuvl m stations) of the new grid 

pts in the circumfereecial direction (MUST be oddV : 
,»p*«#ex stxes>Md.»e station* 

FAC * the first 'grid spec leg in DfSTRI 
lt£ « X leading edge 
A RBWING > 0 , arrow wing 

KW - nuwber of points in the wraparound direction ea W Witte* 

NOnuafeer of cut in the span wise direction 

NtJ* number of points in the upper part of the wing 

Nt*nueber of points in the lovef pert of the wing 

Read surf grid dimensions (nsec x opts x 1) 

NAMELIST /WING/ NSEC , NPTS , FAC 
READ (40, WING) 

WRITE (*, "INC) 

Reed the input grid 


54 

“Jr 

56 

W ioo 


Setup distribution of cross-section* to be ' obtained ' < xd is t j 
XRLE-XLE(l) 

XRTE-XLE ( 1 ) + CHORD ( 1 ) 

WRTE-XLE (NC ) +CHORD < NC ) 

IF(XRTE.GE.WRTE) THEN 
XRT - XRTE 
AXRWING - -1. 

ELSE 

XRT - WRTE 
ARRWING - 1 . 

ENDir 

WAKE - AMIN1< XRTE, WRTE) 

DO 100 J-l , NSEC 

XDIST< J ) -XRLE+ ( XRT-XRLE ) * ( FLOAT < J-l ) /FLOAT < NSEC- 1 ) ) 

CONTINUE 

KTtP- (NPTS+1 )/2 

Tbs nose of the wing 
DO 1S7 K-l ,NPTS 

XOUT ( 1 , K )*XDIST( 1 ) 

TOUT ( 1 , K ) •= YB AS E < 1 , 1 , 1 ) 

ZiOUT ( 1 , K ) * ZBA5E (1,1,1) 

CONTINUE 

Begin main loop for each x-s ta tion 

DO 1 000 L- 2 ', NSEC :-w 


XLOCAL-XpiST(L) 

Redistribute the points free span wise cut to stresarwihe out* - 
The output ztitST sod ynbw are Utm the toot to the "tip; therefore 
then doing the lower surf see, need to rehrrthge the argument. 

The output <ZDIST,YNEN) In both surfaces have KTIP I of pta in 

the cl rcunfarenc ial direction, th*ir last point have the sane physical 

talus for both suf faces. 


91 

92 
91 

94 

95 

96 400 


I til c 

102 C 

103 

1 04 

105 
io6; 

id/ *Bagl 
108 

1 09 *Ref-' 


Do the lower *urf see ■ ' - - 

CALL REDlSf < XLOC At , KTIP , FAC , 1 ) 

DO 300 K«1,KTIP 

XOUT ( L , K ) - X LOCAL 
YOUT ( L , K ) - YNEW(R) 

ZOUT ( L , K ) - ZDIST(K) 

CONTINUE 

Po ' the ' upper' surface 

CALL REDIST(XLOCAl,KTIP,FAC,2) 

DO 400 K-KTIP+1 ,NPT5 
XOUT { L , K ) - X LOCAL 
YOUT(L,K) - YNEW(NPTS-K+1 ) 

ZOUT(L,K) - ZDI$T(NPTS-K+1) 

CONTINUE 

for the computational grid of UPS 3D code 
the wake has to hate two different pta in same ‘ ■’ 
physical location, such that ( Y1 , Z1 )- f Y3 , Z2 ) , Her* ' '' 
the calculation divided into upper and lower parts, 

for saftv sake, set Z1-Z2. 

IF (XLOCAL.LT. WAKE) GOTO 900 
DO 500 K-l , KTIP-1 
Kl-KTIP+K 
K2-KTIP-K 

and model V 1F<ABS(Y0U7(L,K1 >-VoW(L,I2 >) VlSV 3, 08-4) THEN 
I F( ABS ( YOUT< L, K1 ) -YOUT< L , K2 ) ) .LE. 1.0E-2) THEN 
H IF{ABS(YOUT(L,Kl)-YOUT(t,K2)> .LEV 1V0E-5) THEN 

ZOUT(L,Kl)-ZOUT(L,K2) 



ENDIF 

500 

CONTINUE 

900 

CONTINUE 

c 


c 

Proceed i 

c 

1000 

c 

CONTINUE 

c 

Write out i 



SOURCE PROGRAM 

samgrid.f 


date 7/07/94 

PACE# 

TIME 4:18:56 pm 

2 


SOURCE TEXT 


iw-l 

WRITE* 50)NPTS,XW, NSEC 
DO 1234 L*1 , MSEC 

WRITE* 50) ( XOUT ( L, X ) ,K=l,NFTS), 

{ YOUT ( L , X ) NPTS ) , 

( ZOUT ( L / K ) ,K*1 # NPTS ) 

1234 CONTINUE 

C . v:-'-.":*:. ■ • ■';:/!'• ■■■ ■ ■ V?'-' , : . 

C Write out original database in plotid format 


Nl-NU+NL 

WRITE* 11 )N1,NC,KW 

WRITE* 11) ( (XBASE * 1 , 1 , M) , 1*1 ,NU) , (XBASE (1 , 2 »M) 1,-1), 
M-1,NC), 

( (YBASE*I,1,MJ ,1-1, NU> , ( YBASE ( I , 2 ,K) ,I-NL,1,-1), 
M*1 ,HC ) , 

( ( Z.BASE * 1 , 1 r M) , 1*1 , NU ) , ( ZBASE ( 1 , 2 ,M ) , I*NL, 1 , “1 ) , 
M-1,NC) 


Read tie fuaalage grid and combine the lueleage with the 
iriag grid to form a whole configuration* 

CALL WBGRlt) 

ftead the nacelle* grid and combine the awellas with the : 
viog-hod y gr id . 

. . ; CALL HACCRTP v: r.V 


CLOSRilO) 

CLOSE (40) 

STOP 

END 





i JliiJ -X*® ... 

FORTRAN Program. 


UN£ * 


] 56 

IS? 

rss 

m 

T6l 


SOURCE PROGRAM 

samgrid.f 


DATE 


7 / 07/94 


PACE # 


I 03 

~m 

16$ 

166 

T67 

l?0 

(?) 

m 

m 

m 

175 

176 

177 
176 
17$ 
166 
161 
167 

163 

164 

165 

166 
167 
166 

| I 99 
191 

m 

m 

195 

196 

197 
196 
199 

291 

204 

295 

296 


10 


50 


MAX is less than NP1 ' 


..SUBROUTINE ^^MID(NPLI ^NPL2 ‘x ' Y ^ ^pi ^SEC ‘eDImJ * **■**!*& *' ^ 

and WLar.^'S3 Btxdtovise section 

PARAMETER (MAX- 460 ) ‘ 

DIMENSION YTEMP ( MAX 1 ZTEMP/MA* i 
DIMENSl <* 

IE(MAX. Le.NPIj then ' : ' • • • • 

WRITE ( * , * ) 'SUB ADDGRID 
ENDir 

IF(NPL1 .GE.NPL2 ) THEN 

^ITEC,.) No pl« De ls >dd(Kl u the streM|wlse dlr(ctl 

ENDI F 

th * “ v gzU ' •« ** it i» . twpor«ry txxtv 

*2 - Z(NPL2 ) 

S ■ °- 5 *(X(M> l 1)*X ( NPU)) 

DO 10 K-1,KDIN 
Y1 - Y(NPL1,K) 

Y2 - Y ( NP L2 , K ) 

21 - Z(NPL1,K) 

Z2 - 2 (NPL2 , K ) 

CALL LININT(X1 ,X2,Y1 ,Y2 XX Y\n 

ZTEMP(K) - ZZ 
CONTINUE 

- 

DO 30 L- NSEC,NPL2+1 ,-l 
*(1) - X(L-l) 

DO 20 K-l , XDIM 
Y(L,X) « Y ( L-l , K) 

Z(L,K) - Z(L-l,in 
CONTINUE ' 

CONTINUE 

*** '**■** ' *•*< . 

X(NPL2) - XX 
Y(NPL2,E) - YTEMP ( K 1 
Z(NPL2,Kj - ZTEMP(K) 

CONTINUE 1 1 

RETURN 
END 




WmM 

FORTRAN Prppr 



SOURCE PROCRAM 

samgrid.f 


[ LINE# 


• ,****»» • *** .*;**. ** V!*> RFIL, ARCORR) 

SUBROUTINE CIRCLE(*S,KE,RMAX,Y,Z,RFIL,a* 

DIMENSION Z(BMhYjKMWI) t *i; . .,*«**; *W* ***** J£* tbS^f il«t 

SS*^ro5t!L%SUi'Ui^i»t. « ** 

polo ts oo m«t Cir<fl« 

Loo* !«x tM *am :* ■«• uut cUcl * (YC ' IC1 

YA-Y(KS) 

ZA-Z(KS) 

YB-Y(KE) 

ZB-Z(KE) 

SY-YA-YB 

R -SY/SZ , x , . ■>,- • '•■- jY:;:/,; ■ 

n l 2*ZB*R - 2*YB - BB*H/SZ 
c - ZB*ZB+YB*YB + (BB/(2.*SZ))**2 

DET-B*B-4.«A*C 

11 »— tMn °’ 

GOTO 200 

Si* I < -B+SQRT(DET) ) / 

ZC1 - BB - 2.*YC1*SY)/(2.*SZ) 
y C 2 - ( -B-SQRt(DET> } / (2.*M 
TC2 - <BB - 2. *YC2*5Y }/ (2 .*SZ) 

IF(ZCl.GE.ZC2) THEN 
ZC-ZC1 
YC-YCX 

else 

ZC-ZC2 
YC-YC2 
endif 

Find Ux« totBl ire 

tot arc- 6. 

DO 50 K-ES , RE-1 

i0 CONTINUE 

Tina U» «c !•»**■* */• ' 

XRC - TOTARC/FLOAT(KE-RS) 

ARC - ARC* ARCORR ^-.-C 

tivA Uw coordinate for each point 

DO 100 I-ES,KE-1 

ya-y(E) 

ZA-Z(K) 

YB-YC 
ZB-ZC 

SY-YA-YB 
SZ-ZA-ZB 

R -SY/SZ 


2 . *YB - BB*R/SZ 


A - 1 .+R 4 *2 
B - 2 . *ZB*F -■ — 

- R-SR/SZ M ... ,„ n ,n - ZB*BB/SZ - RriL**2 

c : - SB ' IB/SI 

DET"B*B-4 . *a*c 

ir S“v!'*™‘»*" 1 1* >•« •• '' DET 

GOTO 200 

ENDir ( - B +SQRT(DET) } / (£;**> kfirM2n/(2 
(BB - 2 - *YC1 * SY + ( rTIL* * 2-ARC !))/{ 

( -B- SORT (DET 5 )/ l # ,. XRC **2n/l2 ■ *SZ) 

ZC2 - (BB - 2.*YC2*SY + (RFIL**2 AFC 
If(YCl -GE.Y(Y) ) Y® 13 * 

Z(E-H)-ZC1 
Y(K-H)-YC1 
ELSE 

Z<*41)-ZC2 
Y (K+l )“YC2 
ENDIF 
CONTINUE 
RETURN 
END 


YC1 

ZC1 

YC2 

ZC2 



LINE * 

~ 297 c* 

298 

299 c 

300 

301 C 

302 

303 c 

304 c 

305 c 

306 c 
302 c 

308 c 

309 c 
3 T 6 c 
371 c 
312 c. 
3(3 c 
3(4 c « 


SOURCE PROGRAM 

samgrid.f 


^7777— - - SOURCE TEXT 

SUBROUTINE CSPLINE< X Y ,N , XNEW f YNEN ,NNEW> - v " 

PARAMETER (NKAX-506) 

REAL X(N), Y(N), XNEW(NNEW), YNEV(NNEN) 

' * *--*:*** *-* * * ** **.* ** * * “*>***.* ******+**** ? 4*4 **4***4** *.V»*V A 

TKS NEW PQI*is M ; ’] * 



PACE# 


5 


, WRITTEN BY JEFF <^D0VA 10/26/6 6 


ll; C REAL D( NMAX), DELINMAX), H(NMAX) 

HI * «;::^*^i , * 4 "‘‘ , ******t-**”**--*. ........... 

c ■ bpline cofcf’WCiJstfr calculations ■■■••'' 

|!f c ** *''****-''***■**'+ :*>. '^I/**^*****/^ 

|2| SPACING AND FIRST DIVIDED DIFfERENCE 

::m: doiooi-i,n-i ^ 

m ad) - x(i+i) - id) 

3f6. 100 CONTINUE 

327... c . . : : . ? : >' ‘ X&l ; :.:.: ::;: : :v::: :-: : ,, .. ............. w ...,. >v . v 

&o 2oo i-i , n-1 "’ :: ■"•■'■ V - 

la? DEL<I) - ( Y( I+l ) - Y(I)J / H(I) 

330. 200 CONTINUE 

331 c 

III G .SttINt COEFFICIENTS 

33.3. C 

q *** INTERPOLATION FOR N-O CASE ••• 

M IF < N eq. 2) WEN * 

Ill D(l) - DEL(l) 

||S D(N) - DEL(I) 

|39. GO TO 399 

340 ENDIF 

341 c *> 

J4| c •«> nodowit tntiffi eotmaswi n< k >. 3 cmi 

5*45 c SSiMS <ME ™* P ° 1 “ T tORWUl M.TEHBD lO BE 

346 C . ' • . ' • ' . • ■> : f j >: ;>/ '■ • :*'-*’y'- :; :/ : Y:/ : :' : Y.ZZx- • '■ ? Y/ v 

* HI ' - 8<1) : 4 : H < i) ' ' x : Y0 // y/Y &$$$:. :,Y 

fig W1 " (H<1) + 1SUK) / HSUM 

fiS "2 - -H(l) / HSUM 

D <1) " N1*DEL(1) 4 N2*DEL|2) 

||1 IF <PCBST(D<1) / DEL{1)) . LE 0.) THEN 

352 D( 1 ) * 0 . 

Ill- EL ?5" (PCHST(DEL( 1 ) .DEL ( 2 ) ) .LT. 0.) THEN 

|| J UMAX - 3 . * DEL ( 1 ) 

Ill ENDIF {ABS(D(1>> GT ' ABS < DMAX >) D(1 ) - DMAX 

"ill INTERIOR POINT8 (BROOLIE MODIFICATION OF RUTLAND FORMULA) 

If? CONST - 1 . / 3 ; ' ; ?. V 

|f I DO 3 00 1-2, N-l 

|-f| TOP - DEL<I-1) * DEL<I> 

- |il T0P " TOP * .5 * (1. 4 SIGN ( 1 . , TOP) ) 

S-xJ ALPHA = CONST * (H<I-I) 4 2 . *H( I J ) / 4 

ill: + <l *^» * 

II 300 CONTINUE 

|fo C (0 ” ™* EB POtMT OUSTED TO BE 

lit RSUH - H(N-2) 4 H(N-l) ' 

|X| R1 - -H(N-l) / HSUM 

3/J N2 - (H(N-l) 4 HSUM) / ISUM 

375 D<N) - Ml*DEL<N-2) 4 H2»DEL(N-I) 

37.0 IF (PCHST<D(N) ,DEL<N-1)) .LE. 0.) THEN 

D(N) “ 0. 

iff: il ££ i r3S53r a, ™' iM iT ^ ^ ■» "v 

" 111" rwnfr fABS<I><MJ * GT D(N) - DHAX 

:;:m: c . F 

| 3M'' c** COMINUE 

' »M s 

Ilf S...; M>UNE EVALOMtON 

.586 c . , ; ; 

Ill C •*«) -Ml. I(N) 

:zm: IENb l ' ■■ ' :v Z-: :• ••;:•• y/''' 

392 . DO 4 00 J-1,N-1 

— SIH! EE “ <J< J ) + D(J4l) - 2 . *DEL{ J ) ) / 

S« . iend *DEL{ j ) - 2 . *D( J ) - D<j4 ln V h{j, < )} 

l|f CRAY " VISRCiFGEt i ) ron CRAY 

IEND - ISRCHCE<NNEW,XNEN,l,X(J+l)> fOnWH 
-Iff DO 500 I-IBEG,IEND-1 

T- 1NEH(I) - X(J) 

40? 500 CONTINUE ^ ' Y<J> + T * (D<J ' 4 T*(CTHO . T-CTHREE)) - 

4§| 400 CONTINUE 

4W c ... . 

C..., r XNEW<I) .GT. I<M) .. 


DO 600 I=IEND,NNEW 

T - XNEW(I) - X(N-l) 

CONTINUE } " Y(N ‘ 1J ^ t *(D(N-1) 4 T*(CTWO + T*CTHREE ) ) 


RETURN 

END 



sI^ i ,Tr <H ' > ' iNcx ' FTARGET > 

IF(N.LE.O) then 
ISRCHGE - 0 

return 

ELSE 

IT * 1 + (N-lj * i NCX 
ISRCHGE - 1 
00 10 I»1 , IT, INCX 

FTARGEt ) goto II 
__^ RCHGE " ISRCHGE + 1 
CONTINUE 
CONTINUE 

endif 

return 

END 


UAi t. 





SUBROUTINE CUSTER 

di«^on S ?«MHP?" , Z™ (HPI ) 
L2 “ I*-l 

NFUS - KDIM - NPTS 
HBOT - NFUS/2 + 1 
NTOP ■ NFUS/ 2 4 I 


DO 900 LL-M1 , L2 And the v«Jte alOhe 

ww 

RETURN 

ENDIF 


po tb* fcottoi fixati 

DO 10 K-l.NBOT 

YINT(F) “ Z ( LE , F ) 

ZINT(F) - Y(LL,F) 

C ^”!! U SQKT <( 1 j LL . NBOT ) -Z < H*£ J » J *. * * 

DO 80 F-l , NBOT 

Y(LE, F) - iWf(F) 

2( LL/F) - YWF(F) 

CONTINUE ..... 

Do the top ***“ 


100 


180 
1 c 

900 


UK(!) - 1 ( LlI K + ( KDIH-HTOP ) ) 

CONTINUE 

W 1 - < KOI M- NTOP) 

C.CL 

DO 180 F“ 1 / NTOP 

Y<LL,F+(FDIM-NTOP) ) - ***<?> 

Z(LL,F4(FUIM-NTOP>) - ENF(F) 
CONTINUE . .m*:**/*^^^ 


CONTINUE 

RETURN 

END 



SOURCE PROGRAM 

samgrid.f 


7/07/94 PACE * 


4:18:56 pm 



SOURCE TEXT 


***** **♦ ♦*****#♦* ***** *:•*•**■*'*;*;** *;** * **•* * ** *;***;*'*•* « * *** »*** ***»**«* * * * * 
SUBROUTINE DISTARC < X , Y , N , XNEW ,YNEW , NNEW , FGS, I FLAT ) 

DIMENSION X<N) ,Y<N) ,XNEW(NNEW) ,YNEW(NNEW) 

this program redistribute the pointe {X,Y) by subroutine DtSTRI : 
ha sod on the axe length, res la the first grid spacing Note that 
the *w3 points of the two se t* are the same . 

I FLAT* 0, grid points will custer near the first point, »L Bear the hind. 

input array is <1(1 j i*l, ,N, 

otttput 4X*ay is <XNE»< i J «YNEH( i ) ), i-I , , MNIBW 


PARAMETER (MAX-400) 

DIMENSION S{ MAX) , TOTARC (MAX) , XN ( MAX ) , YN ( MAX ) 

Maximus number of points allowed is MAX 
IF(MAX.LE.N .OR. MAX. LE.NNEW) THEN 

WRITE{ * , * ) 'SUB DISTARC : MAX is less than N or NNEW' 
STOP 
ENDIF 


LOCk for total art length 
TOTARC(l) - 6. 

DO 10 X-2,N 

ARC - 5QKT{ (X(X)-X(X-l) )**2 + ( Y( X )-Y ( X-l ) ) * O ) 

TOT ARC < X ) - TOT ARC ( K1 ) 4 ARC 
CONTINUE 

Apply subroutine ttSTRI to. obtain ; the stratohing function S 
for FGffiO , equal spacing is used 
iP(F6A.Gf .0. ) THEN 
DELT-FG S/TOT ARC (N ) 

CALL DISTRI ( DELT ,NNEW, 5 , IFLAT) 

ELSE 

S(l) - 0. 

DO 25 X-2 , NNEN 

S(K) - S(X-l) + 1. /FLOAT ( NNEW- 1) 

CONTINUE 

.ENDIF 

Redistribution, put nev array in a temporary arrays aid' VN 

TN(1)-Y(1) 

XN {NNEW )-X(N ) 

YN(NNEW)«Y(N) 

DO 66 J V J,NNEN 

ARCNEN - S ( J ) * TOT ARC ( N ) 

DO 55 X - 2,N 

IF< ABS(TOTAFC<X)- ARCNEN) . IE. l.E-7) THEN 
XN(J) - X(X) 

YN(J) - Y(X) 

GOTO 60 
ENDIF 

IF< TOTARC ( X) .GT.ARCNEW) THEN 
XI - X(X-l) 

X2 - X(X) 

Y1 - Y(X-l) 

Y2 - Y( X ) 

XX - XI + (X(K)-X(X-l))* 

6 ( ARCNEN- TOTARC(K-l) ) / ( TOT ARC ( X ) - TOTARC ( X- 1 ) ) 

CALL LININT(X1,X2,Y1,Y2,XX,YY) 

XN(J) - XX 

IF(ABS(X1-X2) .IE. l.E-7) THEN 

YN(J) - Y1 ♦ ( ARCNEN- TOTARC (K-l)) 

ELSE 

YN(J) - YY 
ENDIF 
GOTO 60 
ENDIF 
CONTINUE 
CONTINUE 


Write the temporary arrays into the output XNEW, yhrew 
DO 70 J-1,NNEN 
XNEW(J) - XN(J) 

YNEW( J) - YN( J ) 

CONTINUE 

RETURN 

END 





SOURCE PROGRAM 

samgrid.f 


SOURCE TEXT 


« V****««4>«* • »•%*+«****•**« ******* ***««* 4* **«»***«*★* •* **44>4**»>« 
SUBROUTINE DISTFI ( FANG, KFCS , S , I FINE ) 

PARAMETER (MAX-400) 

DIMENSION S { MAX ), Dt'M(MAX ) 

.Calculating t&e stretching function S when given 
the first spacing, FANG, and the number of points KFCS 
.1 1 IFINE-1 , distribution is cusltering at outer grid 

I F(MAX . LE . KFCS ) THEN 

WRI TE (*,*)' SUB DISTRI : MAX is less than KFCS* 

STOP 
END IF 


IF(KFCS.EQ. 1 ) THEN 
S(l) - 0. 

GOTO 40 
ENDIF 


DZ1 • FANG 

KFM - KFCS-1 

DZETA - 1. /FLOAT (KFM) 

RDBETA -1.5 

CALL GRBET( DZI , KFM, 0 . 0001 , 100, RDBETA) 
CALL FZ1( KFCS, RDBETA, DZETA, S) 

tr ariNE.EO.l) THEN 
DO 37 K-l , I PCS 

DUM( KFCS-K+1 ) - 5 ( X) 

CONTINUE 
DO 38 K-l, KFCS 
5(K) - l.-DUM(K) 

CONTINUE 

ENDIF 

CONTINUE 

RETURN 

END 







SOURCE PROCRAM 

samgrid.f 


7/07/94 ™ce# 


T, ME 4:18:56 pm 



SOURCE TE XT 

* 4 ***'*+*:* ******* ♦>'* **■’«■*:* *iU*+**-+i4 . * ** t 

SUBROUTINE EDGE { NC ,NU , NL , XL , XBX , XBASE , YBASE , ZBASE , NP X , LS ) 

DIMENSION ZBASE ( NPX , 2 , LS ) ,XBASE(NPX , 2 , LS ) , YBASE { NPX , 2 , LS ] 

ZLE - ZBASE (1,1,1) 

DO 200 X «1,NC 

I F{ XBASE ( 1 , 1 , X ) . GT . XBX ) THEN 
XI - XBASE (1,1#X-1) 

X2 - XBASE ( 1 , L , X ) 

Z1 - ZBASE (1,1, X-l ) 

Z2 - ZBASE ( 1,1,1 ) 

CALL LININT<X1,X2, Zl, Z2 ,XBr,ZBX) 

GOTO 210 
ENDIF 
) CONTINUE 
) CONTINUE 


DO 500 X-l/NC 

IF(ZBASE(1,1 , X) . GT.ZBX) GOTO 700 

CALL LININT( ZLE, ZBX, XL, XBX, ZBASE{ 1 , 1 , X ) ,XLE) 

XLEOLD - XBASE ( 1 , 1 , X } 

XTL - XBASE (NL,1,X) 

DO 280 I-1,NL 

F - XBASE ( 1 , 1 , K) -XLEOLD 
E « XTL- XBASE ( I , I , K ) 

XBASE ( 1 , 1 , X } - ( F*XTL + E*XLE)/(F4E) 
CONTINUE 

XTL - XBASE (NU, I, K) 

DO 300 1-1, NU 

F - XBASE ( 1 , 2 , X) -XLEOLD 
E - XTL- XBASE (I, 2, K) 

XBASE ( 1 , 2 , X) - ( F*XTL + E*XLE)/(F+E) 
CONTINUE 
CONTINUE 
CONTINUE 
RETURN 
END 





^*****mm*«*»m + *h*m****mm*u**m»****mmm 

SUBROUTINE EQSPACE 
include "sgrid.coa" 

L2 - L-l 
DO 130 LL-1,L2 
XIN ( LL) - X(LL) 

DO 120 X-1,KDIM 

ZIN{ LL, K ) « Z(LL,X> 

YIN{ LL, K) * Y ( LL , K ) 

CONTINUE 

CONTINUE 

XTOT - X(L2)-X<1) 

DX - XTOT/FIiOXT{L2~l } 

DO 160 JL-2,L2 

X ( JL) - X{JL-1)+DX 
DO 150 KL-1,L2 

IF<ABS(XIN(KL)-X(JL) ) . LE . l.E-7) THEN 

DO 140 I-l,n>IM 

Z(KL, K] - ZIN(JL,K) 

Y (XL, X) - YIN<JL,K) 

CONTINUE 
GOTO 160 
ENDIF 

IF(XIN{KL).GT.X(JL)) THEN 
DO 145 I-1,XDIM 
XI - XIN(XL-l) 

X2 - XIN ( XL) 

Y1 - YIN ( IL-1 , X ) 

Y2 - YIN<IX,X) 

Z1 - ZXN<KL-1,K) 

Z2 - ZIN ( XL , K ) 

XX - X(JL) 

CALL LINXNT(X1 ,X2 , Y1 ,Y2 , XX ,YY} 

CALL LININT(X1,X2,Z1,Z2,XX,ZZ) 

Y { JL, X) - YY 
Z<JL,X) - ZZ 
CONTINUE 
GOTO 160 
ENDIF 
CONTINUE 
CONTINUE 
RETURN 
END 




SOURCE PROCRAM 

dati 7/07/94 

samgrid.f 

T,ME 4:18:56 pm 



SOURCE TEXT 


687 **** 

688 

689 

690 

6sr 

692 

695 c 

694 c 

695 c 

696 c 
m c 
696 c 

699 c 

700 C 

701 c 

702 c 

703 c 

704 

705 

706 

707 
706 

709 

710 

711 

712“ c 

713 

714 
7TS 

716 

717 
716 

719 c 

720 c 

721 c 

722 c 

723 

724 

725 

727 

726 220 

729 

730 

731 

732 

733 

734 330 

735 

736 

737 130 


******************************* ft'*'******* **'*•*'** * ft**'*** * ft ** * * *;#■* ***** 
SUBROUTINE FILET{ Y , Z, XDIM , MB1 ,HB2 ,HT1 ,MT2 ,RFIL) 

PARAMETER (MAX- 4 00) 

DIMENSION Z(XDIM) ,Y(XDIM) 

DIMENSION D1 ( MAX ) , D2 ( MAX ) 

COMMON /REF/ ZROOT, XTIP , ARCORR 

This subroutine takes a vingr-fusal ago station y <Y(k) , t(k.) > fc-l/KDIM,- 
and find the two ( top And bottom) in teraectionj of the wing nod the 
fuselage. 

And theD, for example, at the bottom intersection (Yflt') ,6(xn , it 
extends to a Segment of points* <Y(h),Z(K)) k.«IFl to XF2, where 
KF1-X-MB1 , XF2-K+MB2 . : : '-f ; - /: 

Similar procedure for the top part* 

And then, call subroutine CIRCXE to replace the eegment by 4 segment 
of a c ircle with radius RFXL . 


I F ( RFI L . EQ . 0 . ) GOTO 735 
IF( MAX . LE . XDIM) TEEN 

WRITE< * , *) ’SUB riLET : MAX is leas than XDIM' 
STOP 
ENDIF 


Bottom part of the Aircraft 
IF(MB1 .EQ. 0 . ANt>. HB2.EQ.0) GOTO 13S 
DO 130 X»1,XDIM 

IF (X.LE.XTIP .AND. Z ( K ). GT . ZROOT) THEN 
IF1-X-MB1 
XF2-X+KB2 

CALL CIRCLE ( KF1 ,KF2 ,KDIM,Y, Z, RFI L, ARCORR) 


Me hawe N-XF2-XF1+1 pts in filet area, employ two more points 
from the original grid end redistribute tW, the grid spacing looks 

■ smoother ; ‘ v ••.••••. - *VJ l V : 

KFl-XFl+i 
XF2-XF241 
DO 220 XD-Kn,KF2 

D1 ( ED- XF1 + 1 ) - Y( XD) 

D2(XD-XF141) - Z ( KD) 

CONTINUE 
N-XF2-KF141 

CALL DISTARC ( D2 ,D1 ,N ,D2 ,D1 , N, -10 . ,0) 

DO 330 ID-KF1,KF2 
Y(XD)-D1< XD-XF1+1 ) 

Z( XD)-D2 ( XD-XF1+1 ) 

CONTINUE 
GOTO 135 
ENDIF 
CONTINUE 
CONTINUE 

Top' part of the aircraft : •: . £v : 

IF(MT1 .EQ 6 .AND. HT2.EQ.0) GOTO 735 
DO 700 X- XTIP, XDIM 

IF( X . GT. XTIP .AND. Z( X) . LE . ZROOT) THEN 
XF1-X-MT2 
KF2-X+MT1 

CALL CIRCLE ( KF1 f XF2 r XDIM, Y, Z, RFI L y ARCORR ) 

We hate ^|CF2-XF14i pta'' in filet ' ‘ 'imra‘''iiiioirm ''poiaiEufr 

from the original grid and redistribute them, the grid spacing 

Smoother . ' :>#*?■*$■* 

XF1-KF1+1 
XF2-XF2+1 
DO 420 XD-XF1 , KF2 

D1 ( XD-XF1+1 ) - Y ( XD) 

D2 ( XD-XF1+1 ) - Z(XD) 

CONTINUE 
N-XF2-XF1+1 

CALL DISTARC (D1,D2 ,N, D1 ,D2 ,N ,-10 . r 0) 

DO 530 KD-KF1 , KF2 
Y( XD )-Dl ( KD-XF1+1 ) 

Z(XD)-D2 (XD-XF1+1 ) 

CONTINUE 
GOTO 735 
ENDIF 
CONTINUE 
CONTINUE 







SOURCE PROGRAM 

samgrid.f 


OATE 7/07/94 I pace# 
time 4:18:56 pm 14 


SOURCE TEXT 


4 4 * * * ♦ * • * 4 * 4 * * « <4 4 4 * * * * 4 * * * * * * * * * 4 ;**' '* * • * * ♦ * * * * * * « 4 > * * » * • • ;* * * * * * • * 4 4 4 
SUBROUTINE f Z1 ( LI , TBETA , DET, 1 ) 

COMPUTES NORMALIZED NORMAL DISTANCE, Z<L) 


PARAMETER (MAX-400) 

DIMENSION Z{ MAX ) 

IF(KAX.LE. LI ) THEN 

WRITE <*,*) 'SUB FZI : MAX is less than LI' 


IF( TBETA. EQ. I . ) THEN 
DO 10 l«l,Ll 
Z( L) -0. 

10 CONTINUE 
ELSE 

DO 20 L* 1 , LI 
ETA-( L-l )*DET 
RR-(TBETA+1. )/(TBETA-l. ) 

EEE-1 . -ETA 
RBETA-RR * * EEE 

Z(L)- { TBETA- 1 . ) * (RR-RBETA)/ (RBETA+1 . ) 
20 CONTINUE 
END IF 
RETURN 
END 





SOURCE PROCRAM 

samgrid.f 


SOURCE TEXT 

SUBROUTINE GRBET(DFM*NPT\ ^CC *ICC*bETA ) ^ ^ ***** *'* ******************* 

STRETCHING PARAMETER , H.TA, 

PARAMETER (MAX- 4 00) 

DIMENSION Z(MAX) 

IF(MAX.LE.NPT) THEN 

WRITE ( * , * ) ' SUB GREET : MAX is less than NPT' 

STOP 

END!F 

ICCL-ICC 

FPCC L- FPCC *DFM 

BETA1-BETA 

Zl-DFM 

DET-1 ./NPT 

BR-1. 

FF— Z1 
IICC-ICC/10 
DO 10 I-l,IICC 
BF-BETA1 

BETA-0. 5* (BET Al+1. ) 

CALL FZ1 ( 2 , BF , DET , Z ) 
rr-Z(2)-Zl 

IFfFF.GT. 0. ) GO TO 15 
BETA1-2 . *SETA1-1 . 

0 CONTINUE 
5 CONTINUE 

DO 5 NIT-1, ICCL 
CALL FZ1 ( 2 , BETA , DET , Z ) 

F-Z( 2 )-Zl 

I F ( F . GT . 0 . ) THEN 

rr-F 

BF-BETA 

ELSE 

FF-F 

BR-BETA 

END IF 

BETA- 0 . 5 * ( BF+BR J 

IF( ABS ( F) . LT. FPCCL ) GO TO 4 

1 CONTINUE 

" " ) Mr*, f ; 

cSmNUE 80 ' 368 EICEEDED MAX N0 - 0F ITS BETA, F ,3G13.6) 

. F-E< 2)-Zi ■ . kiiim ■ 

BM1-BETA-1. 

BFM1-BF-1 . 

BRM1-BR-1 , 

RETURN 

END 











SOURCE TEXT 


CONTINUE 

IF ( INTRAI L . EQ . 2 J THEN 
DO 200 E«l, FOUNT 

IF(ZOUT<L,NPAIR<l,K)) . LT . ZNAC (NTEMP > ) THEN 
EE - E-l 

IF (YOUT ( L, NPAIR ( 1 , EE ) ) . LT . YNAC ( NTEMP ) ) THEN 
NTEKP - NTEMP - I 
GOTO 190 
ENDIF 
GOTO 215 
ENDIF 
CONTINUE 
CONTINUE 
NE - NTEMP 

I F ( NTEMP . NE . NNACP ) THEN 

CALL LININT(YNAC(NE-1) ,YNAC(NE) , ZNAC (NE-l) , ZNAC (NE) r 
YOUT ( L , NPAIR { 1 , EE ) ) , ZZE ) 

YYE ■ YOUT { L , NPAIR ( I , EE ) ) 

ELSE 

CALL LININT( ZOUT { L , NPAIR ( 1 , EE ) ) , ZOUT ( L ,NPAI R ( 1 , EE ) + 1 ) , 

YOUT < L, NPAIR ( 1 , EE ) ) f YOUT (L , NPAIR ( 1 , EE ) + l ) , 
ZNAC(NE) ,TYE ) 

ZZE - ZNAC(NE) 

ENDIF 

ELSE 

DO 230 E-1,NPH 

IF( ZOUT( L, K) , GT. ZNAC( NTEMP ) ) THEN 
EE - I 

IF < YOUT ( L, EE ) .LT. YNAC (NTEKP)) THEN 
NTEMP - NTEMP - l 
GOTO 190 
ENDIF 
GOTO 235 
ENDIF 
CONTINUE 
CONTINUE 
NE - NTEMP 

CALL LININT( ZOUT ( L , EE ) , ZOUT ( L , ICE- 1 ) f TOUT { L , EE ) , YOUT( L, EE- 1 ) , 
ZNAC(ME) , TYE) 

ZZE - ZNAC(NE) 

ENDIF 

ttws lover sur faca ; 

Store vingroot grid 
DO I»1 , NPAIR ( 1 , XS ) 

YWG(E) - YOUT(L,K) 

ZWiG(E) - ZOUT { L , K } 

ENDDO 


Store vlngtip grid to * bpp irity ; 

IF( INTRAlL.EQ.2) T&EN 
DO K-NPAIR(1,EE) ,NPH 

YINT(E-NPAIR(l,KE)+l) - YOUT(L,X) 
ZINT(E-NPAIR(l,RE)+l) - ZOUT(L,X) 

ENDDO 

ELSE 

DO K-KE,NPH 

YINT(E-IE+1 ) - YOUT(L,I) 

ZINT(E-EE+1) - ZOUT ( L , E ) 

ENDDO 

ENDIF 

Get t be nacelle grid under the wake line ready 

NADD - NE-NS+1 +2 

YWE(l) - YYS 

Z WE ( 1 ) « ZZS 

YWK(NADD) “ YYE 

ZWX(NADD) - ZZE 

DO E*NS ,NE 

YWE( E— NS+2 ) - YNAC (I) 

ZWE(E-NS+2) - ZNAC(I) 

ENDDO 

IF( INTRAIL.EQ.2) THEN 

NNACC - NPWNH-NPAIR(1,KS)-(NPH-NPAIR(1,EE)+1) 
ELSE 

NNACC - NPWNHNPAIR ( 1 , ES ) - ( NPB- J3E+1 ) 

ENDIF 

CALL DISTARC< YNE , ZWE ,NADD, YWX , ZNX , NNACC ,-10,0) 
Stick t be nacelle grid under the vieg 
bo E«l, NNACC 

ZHNG (NPAIR (1,XS)+K) - ZWK(X) 

YWG (NPAIR (1,ES)+K) - T*K(E) 

ENDDO 

Put back tke wlogtip grid Uto the wing 

IP( INTRAIL.EQ.2) THEN 

EREP - NPB-NPAIR ( 1 , EE )-f 1 
DO K-1,EREP 

ZWNG (NPAIR ( 1 , KS } +NNACC+X ) - ZINT(E) 

YWNG (NPAIR { 1 , K5 ) +NNACC+E ) - YINT(E) 

ENDDO 

ELSE 

EREP - NPH-EE+1 
DO K* 1 , EREP 

ZWNG ( NPAI R ( 1 , E5 ) +NNACC+ E ) » ZINT(E) 

YWNG (NPAIR ( 1 , ES) +NNACC+E ) - YINT(X) 

ENDDO 

ENDIF 

NHNACH - NPKNH 


Fore the upper surface 


Store the Wingtlp grid 
IF( INTRAIL . EQ. 2) THEN 
DO K*NPH , NPAIR ( 2 , EE ) 

YWNG(NWNACH+E~NPH) - TOUT(L,E) 
ZWNG (NWNACH+E-NPH ) - ZOUT(L,E) 
ENDDO 

N1 - NHNACH+NPAIR( 2 , EEJ-NPH 
ELSE 

DO E-NPH,NPTS 

IF(ZOUT(L,E) .GT. ZZE) THEN 
XN1 - E 
GOTO 260 
ENDIF 
ENDDO 
CONTINUE 
DO K«NPH,ENI 

YWNG(NWNACH+K-NPH) = YOUT(L,K) 
ZWNG(NWNACH+K-NPH) - ZOUT(L,K) 
ENDDO 

Nl = NWNACH+KN1-NPH 






SOURCE PROCRAM 

DATE 7/07/94 

samgrid.f 

TIME 4:18:56 pm 



SOURCE TEXT 


; Store tk* vingrpot arid to « tew> array 
DO X-NPAIR < 2 ,KS) , NPTS 

YINT(K-NPAIR(2,KS)4l) - YOUT(L,K) 

ZINT( K-NPAIR ( 2 , KS ) 4 1 ) - ZOUT(L,K) 

ENT) DO 

Oat the oAcelle grid above the Wake line ready 

NADD - NNACP-NE 4 NS-1 +2 

YNK<1) - YYE 

ZWK(l) - ZZE 

DO K«NE41,NNACP 

YWK(K-NE+1 ) - YNAC(K) 

ZNK(K-NE4l) - ZNAC(K) 

ENDDO 

DO K-1,NS-1 

YWK<NNACP-NE414K) - YNAC(K) 

ZWK(NNACP-NE4l4K> - ZNAC(K) 

ENDDO 

YNK { N ADD ) - YYS 
ZWK (NADD) - ZZS 

CALL DT STARC ( YWK , ZWK , NADD , YWK , ZWK , NNACC ,-10,0) 
Stick the nacelle grid above the ving 
DO K- 1 , NNACC 

ZKNG<N14K) - ZWK(K) 

YWNG(Nl4X) - YWX(K) 

ENDDO 

Put back the vlogroot grid into the wing 
KREP « NPTS-NPAI R ( 2 , KS ) 4 1 
DO K*1,KREP 

ZWNC(N14NNACC4I) - ZINT(K) 

YWNG(N1+NNACC4K) - YINT(K) 

ENDDO 
GOTO 700 
ENDIF 


«• axe before trailing v MWWl redistribution 

CONTINUE 

; ; Rottoa ■part.’ ' " ' ' ;; v -: 

nbd out tie points that being replaced i>y the nacelle 
NTEMP « 1 
CONTINUE 
DO 440 K-1,NPB 

ir(ZOUT(L,K) .GT.ZNAC(NTEMP)) THEN 
KS - K-l 

rl - abs{ zout( L, IS )-zout( L , KS-1 ) ) 
r2 - abs{ zout( L, KS )- ZNAC (NTEMP ) ) 
if (r2. le.rl/8. ) then 
XS-IS-1 
end If 

IF<ABS<Y0UT< L,KS)-YN AC (NTEMP) ) . LE . RADNAC/800. ) THEN 
NTEMP - NTEMP 4 1 
GOTO 430 
ENDIF 

CALL LININT ( ZOUT( L, KS ) , ZOUT( L , KS41 ) , YOUT ( L , KS ) , 
YOUT(L,KS41) ,ZNAC( NTEMP) , YYS ) 

ZZS - ZNAC( NTEMP) 

GOTO 445 
ENDIF 
CONTINUE 
CONTINUE 
NFRIS - NTEMP 

NTEMP - NNACP 
CONTINUE 

DO 460 K-KS4l,NPH 

IF<ZOUT<L,K} .GT. ZNAC< NTEMP ) ) THEN 
KE » K 

IF(ABS(YOUT( L,KE)-YNAC< NTEMP) ) . LE . RADNAC/800 . ) THEN 
NTEMP « NTEMP - 1 
GOTO 450 
ENDIF 

CALL LININT ( ZOUT ( L , KE- 1 ) , ZOUT { L , KE) , YOUT< L, KE-1 ) , 
YOUT { L , KE ) , ZN AC ( NTEMP ) , YYE ) 

ZZE - ZNAC (NTEMP) 

GOTO 465 
ENDIF 
CONTINUE 
CONTINUE 
NFRIE - NTEMP 


Store vingroot grid 

DO 500 K-1,KS 

YWNG(K) - YOUT(L,K) 

ZNNG(K) - ZOUT ( L , K ) 

CONTINUE 

Get the nacelle grid under the vake line ready 

N2 - NFFIE-NFRIS+1 +2 

YINT(1 ) - YYS 

ZINT( 1 ) - ZZS 

YINT(N2 ) » YYE 

ZINT (N2 ) - ZZE 

DO 520 K-NFRIS, NFRIE 

YINT(I-NFRIS42) - YNAC(K) 

ZINT( K-NFRIS42 ) - ZNAC(K) 

CONTINUE 

NNACC - NPNNB-KS-(NPH-KE4l) 

CALL DISTARC ( YINT , ZINT, N2 , YWK , ZWK , NNACC ,-10,0) 
Stick the nacelle grid under the Wing 
DO K=l, NNACC 

ZNNG(KS4K) - ZWK(K) 

YWNG(KS4K) - YWK(K) 

ENDDO 

Put the vingtlp grid 

DO 540 K-KE,NPH 

YWNG(KS4NNACC+K-KE+1) - YOUT ( L , K ) 
ZWNG(KS4NNACC4K-KE+1) - ZOUT(L,K) 

CONTINUE 


Upper part 

DO 5*0 K«NPH,NPTS 

IF ( ZOUT ( t , K ) . LE . ZNAC (NFRIE ) ) THEN 
NN1 - K-l 

CALL LININT( ZOUT(L,NNl4l ) ,ZOUT( L,NN1 ) ,YOUT( L,NN1+1 ) , 
YOUT ( L , NN1 ) , ZNAC ( NFRI E ) , YY1 ) 

ZZ1 - ZNAC (NFRIE ) 

GOTO 595 
ENDIF 
CONTINUE 
CONTINUE 



SOURCE PROGRAM 

samgrid.f 


7/07/94 PACE# 


time 4:18:56 pm 


SOURCE TEXT 


DO 600 K”NPH ,NPT5 

IF<ZOUT(L,K) .LE.ZNAC(NFRI£} ) THEN 
NN2 - X 

CALI LININT< ZOUT( L,NN2 ) , ZOUT< L,NN2-1 > , YOUT( L,NN2 ) 
YOUT( L,NN2-1 ) , ZNAC(NFRIS) , YY2 ) 

ZZ2 « ZNAC (NFRI5) 

GOTO 605 
ENDIF 
CONTINUE 
CONTINUE 

Store the wing tip grid 

DO 620 X«NPH,NN1 

YWNG(NPWNH+K-NPI) - YOUT(L,X) 

ZWNG ( NPWNH+ K-NPH ) - ZOUT(L,X) 

CONTINUE 

R«dl»tlbute the point* above the nacelle 

N2 - NN2-NN1+1 

YWX{ 1 ) - YY1 

ZWK { 1 > * ZZ1 

YWX(N2) - YY2 

ZWX(N2) - ZZ2 

DO r-NNl+l f NN2-l 

YWK(K-NN1+I) - YOUT<L,I> 

ZWX(X-NN1>1) - ZOUT ( L , I ) 

END DO 

NNA - NPWNH-(NN1-NPH'H)-<NPTS-NN241) 

CALL DISTARC( YWX , ZWX,N2 , YWX , ZWK ,NNA , -10,0) 

NNT - NPWNH+< NN1-NPH) 

DO X»1,NNA 

YWNG ( NNT+ X ) - YWX(X) 

ZWNG(NNT+K) - ZWX(X) 

ENDDO • ■ v r 

£ tor* ' th« viegxbdt grid.' . 

DO X-NN2,NPTS 

YNNG<NNT4NNA+K-NN24l> - YOUT<L,X) 
ZHNG{NNT+NNA+K-NN24lJ - ZOUT(L,X) 

END DO 
CONTINUE 

IF(IPRNT.NE.O)THEN 

WRITE ( IPRNT} ( X0UT( L, 1 } , X“1 , NPWN ) , 
<YWNG<X),J=1,NPWN), 

(ZWO(K) ,X«1,NPWN) 

CALL FLUSH (IPRNTJ 
ENDIF 

DO 750 X-1,NPWN 

XOUT( L, X) - XOUT{L,I) 

YOUT<L,X) - YWNG(E) 

ZOUT<L,X) - ZWNG<X) 


CONTINUE 

IF(UJUM(1) .EQ.LNUM(I) -AND. LNUK< 2 > EQ. LNUH< 4 ) ) THEN 
NPTS - NPTS - MC 
ENDIF 
RETURN 






SOURCE TEXT 


*:#'« » ** 'if**'* ;*'** **'■*>** ** * •* ***** * ******* * * *'*;* *** * **'*'* * > * * j*t **'**'*! ♦ « M' 

SUBROUTINE NACGHID 
include "sgrld. cow” 

C This subroutine read the nacelle grid and combine it with the 
C wing grid. 

COMMON /ENG/ XENG < 2 , NPt ) , YENG* 2 , WpI ,HPt) , fcfcNG (2 ,NPI , NPI ) 

S ,MNAC<2) ,MNACP(2) 

COMMON /NACG/ XNAC(NPI ) ,YNAC(NPI ,NPI ) ,ZNAC(NPI ,NPI ) 

COMMON /NDIM/ NNAC , NNACP , LNUM ( 4 > 


(XENG^YENG, SENG ) 
WAC{*> 


Coordinate b£ engine 
Nuaobetpf stations in nacelle 
Number of points in each station 


MNACP{*) Number of points In each 

(1 , *> * ) in«*r nacelle; (2,* , * y outer nacelle 


Read the 'napellea geometry t 
Inset nacelle geometry 

CALL NACIN 
XIN1 - XNAC(l) 

XIN2 - IN AC (NNAC) 

MNAC(l) - NNAC 
MNACP(L) - NNACP 
DO 100 L*1 ,NNAC 

XENG ( 1 r L ) - XNAC(L) 

DO 80 K-l, NNACP 

YENG ( 1 , L , K ) - YNAC ( L , K ) 
ZENG ( 1 , L i X ) - ZNAC(L,K) 
CONTINUE 
CONTINUE 

{mter nacelle geometry 
CALt NACIN 
XOUT1 - XNAC(l) 

XOUT2 - XN AC (NNAC) 

KNAC ( 2 ) - NNAC 
MNACP ( 2 ) - NNACP 
DO 200 L-l , NNAC 

XENG ( 2 / L ) • XNAC(L) 

DO 180 K-l, NNACP 

YENG ( 2 , L, X ) - YNAC(L,K) 
ZENG ( 2 , L , K ) - ZNAC<L,K) 
CONTINUE 
CONTINUE 


In i cm of 2 necelles , three eon* will be made. 

NEwtNG , station < a) will be «Miad io the wing at inlet and outlet 

\ <►£ the nacelle . 

MOUNT : Mount the nacelle under . the wing and/or wake line 


inner nacelle 
Outer nacelle 


The lirat tone, oniy oae nacelle appear* 

WRITE (*,*)' zone 1' 

JN - I 

CALL NEWING(LNUM,XIN1 ,XOUTl ) 

CALL MOUNT (JN,LNUM, 21) 

INUM ( 3 ) - LNUM<1) 

LNUM ( 4 ) - LNUM( 2) 

The second none consists two needles appear 
NRlTE(* / *)’zone 2‘ 

JN - 1 

CALL NENING ( LNUN , X0UT1 , XIN 2 ) 

LNUM(1 ) - LNUM ( 1 ) + 1 
CALL MOUNT (JN f LNCM,0) 

LNUM ( 3 ) - LNUM(l) 

LNUM<4)- LNUM(2) 

WRITE ( * , * ) ' zone 2' 

CALL NEWING<LNUM,XOUTl,XIN2) 

LNUM ( 1 ) - INUM(l) + 1 
JN - 2 

CALL MOUNT ( JN , LNUM ,22) 

INUM(3) - LNUM ( 1 ) 

LNUM (4 ) - LNUM( 2 ) 

The third itohe, only CSS nacelle appear* ' - ^ 
WRITE( * , * ) ’ zone 3* 

JN - 2 

CALL NEWING (LNUM, XIN2 , XOUT2 ) 

INUM(l) » LNUK(l) + 1 
CALL MOUNT ( JN , LNUM ,23) 





SOURCE PROGRAM 

samgrid.f 
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™ E 4:18:56 pm 



SOURCE TEXT 


4 4 4 * 4 4*4 * r * *4 4 4 *4 *4 4 4 4 4;* 4 4 4 * 44 4 4'*4 * • * 4 * 4 4 * * 44444*4***44**#** *444** 4444*44 4 
SUBROUTINE NEWI NG( LN UM,XX 1,XX2 ) 

Include "agrid.coi" 

DIMENSION XWNG { NPI ,NPI ) , YWNG ( NPI , NPI ) , ZWNG(NPI , NPI ) 

DIMENSION LNUM( 4 ) 


Rewrite the coordinate of the wing 
DO 10 L-1,NSEC 
DO 10 X-l f NPTS 

XWNG(L,K) * XOUT(L,K) 

YNNG(L,K) - YOUT( L, K ) 

ZWNG ( L , X ) ■= Z0UT(L,X> 

CONTINUE 


find out where the x- location of start and end of nacelle 
XNSTRT - XXI 
XNEND - XX 2 

Add tare stations in the wing r these two stations lie exactly oa 
XNSTRT and' RHUMB.' ‘ ^ ;■ 

DO 55 NN-1,2 

IF(NN.EQ.l) XX-XNSTRT 
IF(NN.EQ.2) XX-XNEND 
DO 50 LW-1,NSEC 

The wing section is Terr close to nacelle* s station (XNSTRT or XNEND) 
I P ( AfiS ( XWNG < t*, 1 ) -XX ) . IE . 1 . E-? ) THEN 
DO X-1,NPTS 

XWNG( LN f X) - XX 
ENDDO 

LNUM(NN) - LW 
GOTO 55 
END IF 

Create an extjra ' station In the : wing' ■ 

IF( XWNG ( LW , 1 ) . 6T . XX ) TEEN 
X(LW) - XX 
LNUH(NN) - LW 
DO K«1,NPTS 

XI - XWNG { LN— 1 , K ) 

Y1 - YWNG(LW-1,I) 

Z1 - ZWNG(LW-1,I) 

X2 - XWNG ( LW , K) 

Y2 - YWNG(LW ,K) 

Z2 - ZWNG ( LW , K ) 

CALL LININT( XI f X2 , Y1 , Y2 , XX , YY ) 

CALL LININT{Xl r X2,Zl, Z2,XX,ZZ) 

Y ( LW , X ) * YY 
Z(LW,X) - ZZ 
ENDDO 

If we are la wake < Bake sure top pts Intersect the bottoe pts 
DO 24 *i«2, <NPTS+l)/2 
DO 20 X2-(NPTS+l)/2+l,NPTS 

IF ( ABS{Y ( LW, XI ) -Y ( LW, 12 ) ) . LT.l.E-6 .AND. 

L ABS(Z(LK,X1)-Z(LW,K2) ) .LT.l.E-A ) THEN 

Y( LW , XI ) - Y ( LW , 12 ) 

Z( LW r XI ) - Z ( LW , K 2 ) 

GOTO 24 
END IF 
CONTINUE 
CONTINUE 

Put the rest of the station into (X,Y,Z) 

DO 35 L-LW,NSEC 
DO 30 X-1,NPT5 

X<L4l) - XWNG ( L , X ) 

Y ( L+ 1 / X ) * YWNG ( L , X ) 

Z(L+1,X) - ZWNG(L,K) 

CONTINUE 
CONTINUE 
NSEC - NSEC + I 
DO 45 L-LW,NSEC 
DO 40 X-l r NPTS 
XWNG { L , E ) - X ( L ) 

YWNG<L,X) - Y ( L , K ) 

ZWNG(L,r> - Z ( L , K ) 

CONTINUE 
CONTINUE 
GOTO 55 
ENDIF 
CONTINUE 


► 1 55 

jr 

CONTINUE 


it: 

DO 910 L- 1 , NSEC 

4 

DO 910 1-1, 

NPTS 

T' 

XOUT ( L / X ) 

- XWNG( L, K) 

re 

YOUT(L,K) 

- YWNG ( L , K ) 

>1 

ZOUT(L,K) 

- ZWNG(L,X> 

►£ 910 

CONTINUE 



RETURN 


io 

END 






SOURCE PROCRAM 

samgrid.f 


7/07/94 PACE # 
4:18:56 pm 



■* * * **** * * sa’to;*.*'* *** * 4* *'ir'i'**V*>>* 

SUBROUTINE NOSE { FNBOT , FNTOP ) 
include “sgrid.con" 

DIMENSION Dl(NPI) ,D2(NPI) ,S(NPI) 


SOURCE TEXT 


From the nose to th« lending edge 
ie., froe station 1 to Ml 

Elf st point ' oi tie no«ie :: ' :: 

DO 12 K«1,KDXH 
X(l) - XIN(l> 

Y(1,K> - YIN(1,1> 

Z(1,X} - ZIN{ 1,1) 

CONTINUE 


loop for nil stations, from station 2 to station m 
DO 500 M«2,M1 
L - M 

X LOCAL - XIN(M) 

Store the Input to dmany arry 
DO 13 ' 1C- 1 , N'FP 

YINT<K) • YIN<M,K) 

ZINT<K) - ZIN(H, K) 

CONTINUE 

YREF is yiltte of the first point of the wing? - 
XIU6f is the corresponding index of each station: 
YREF*Y0UT( 1 ,NPTS) 


C OPTIONS 

c 


Boeing Basel toe conf Igurstion 

DO 43 K-l r NB> 

IF(YINT(K) .GE.YREF) THEN 
KREF-K 
GOTO 44 
ENDIF 
CONTINUE 
CONTINUE 

KREF-NFP/2 +KREFADD 


: C -™ 

C 




Lov Boom Configuration 


lower part of the noie l-Y to Y«mry 
XS-1 
JTE-KREF 
KN-KE-KS+l 
DO 74 K- IS,KE 
KK - MS+1 
DlfKX) - YINT(K) 

D2(KI) - ZINT( K) 

CONTINUE 

CALI DISTARC (D1,D2,XN,YNEW, ZDI ST, XTI P, FNBOT ,L) 

DO 85 K«1,KTIP 
Y<L,X) - YNEW(K) 

Z(L,K) - ZDIST( K) 

X(L) - Z LOCAL 
CONTINUE 


From Y*YREF to pos Y 
KS-KREF 
KE-NFP 
KN-KE-rS + 1 
DO 100 r- KS,KE 
IK » K-IS41 
Dl(ZZ) - YINT(X) 

D2(KK) - ZINT{ K) 

CONTINUE 

CALL DISTARC(D1 ,D2 , KN , YNEW, ZDI ST, XTIP , FNTOP , 0 ) 
DO 400 K*KTIP , KDIM 

Y ( L , X) - YNEN(X-KTIP+1) 

Z(L,X)» ZDI ST ( K-XTIP+1) 

Z<L) - X LOCAL 
CONTINUE 


500 CONTINUE 
RETURN 
END 



SOURCE PROGRAM 

samgrid.f 


OATt 7/07/94 ' 

time 4:18:56 pm \ 


*»««&..■•■■• *" ^ ■•■«•■•■■- v " ■ ' "" ■■■■ "— ■ TCV _ 

SOURCE TEXT 

SUBROUTINE REDI ST< XLOCAL , FTIP , FAC , I FLAT ) 

include "sgrid.con" 

This is the stain subroutine to redlatributethe points fro* 
span vise cut to streaavise cut. : • , 4J . ' ■, 

Hben there is *cheuog~SIKP'* , it neaas the output is tQt SIMP <*>de 

PARAMETER (INT-4 00} 

DIMENSION S(INT) 

NU « aushw of point is tha upplr 

jlL - rn.h.r nf jpoiht in the lover • surface :••••• \ > ; 

HC - nusbei of spanuise sections 

^21 NT# TINT) e^poiot of atreawvise cut for ** Jt?* *?5**°* 

KT * I or pte in circu». direction extracted from the old grid. 

FTIP * I of pts in the circus, dijfeotign in one surface 

I FLAT - 3 i uppd* surface 

' ' I FLAT « 1 j lover surface 

IF(INT.iE wmvt ., 

WRITE(*,*)' STO REDIST t INT is less than KTIP 

STOP 
END IF 

IF( IFLAT. EQ. 2 ) THEN 
NUL - NU 
ELSE 

NUL - NL 

ENDI F .... ,... .....v • ■»»: : : ; v ; .v,. : „ 


the leading edge 
' I F< XLOCAL . GT.XLE <NC JJTHEJf . . . 

:; ; X-athtlo* is at ’ the wing ''tip ' ''yS^li S ! '• 

I F( XLOCAL. LE . (Xt£ (NC ) +CHORD(NC) ) )TffiN. 

Should cell NtNOtorAXt if wake is contained in this station 
I F( XLOCAL - GT . XLE ( 1 ) +CHORD( 1 ) } THEN 
CALL WINCNNAXE( XLOCAL, FT, NUL, IFLAT} 

GOTO 200 

ENDI F ... . . . . . . ............. 

distance 1$ between theleadlng edge and the trailing edge 
DO 60 M-1,NC 
DO 50 1-2, NUL 

I F< XBASE (I , IFLAT,NC-H+1) .GE. XLOCAL 
i ABS (XBASE ( I , IFLAT,NC-M+I ) -XLOCAL} . LE . 1 . E-7 ) THEN 

XI - XBA5E( 1-1,1 FLAT ,NC-M+1 ) 

X2 - XBASE (I , I FLAT ,NC-M+ 1 ) 

Yl - YBASE(I-1,IFLAT,NC-M+1) 

Y2 - YBASE ( I , IFLAT ,NC-M+1 ) 

21 - 2BASE (1-1,1 TLAT , NC-H+ 1 } 

22 - 2 BASE (I , IFLAT ,NC-N+1 ) 

CALL LININT(X1 ,X2 , Yl ,Y2 ,XLOCAL, YT ) 

CALL LININT(X1 ,X2 , Z1 ,Z2 , XLOCAL, 21) 

YINT(M) * YY 

2INT(M) - ZZ 
GOTO 60 
ENDIF 
CONTINUE 
CONTINUE 
KT-NC 
GOTO2O0 


The x- station passes the wing tip# should have vafce 

CALL WINGNWAXE( XLOCAL, XT, NUL, IFLAT) 

GOTO 200 
ENDIF 
ENDIF 


Streanvise distance is i« the leading edge 

iF(xi^al: u.ite (He) j then 

Should call NTNGNWAKt if wake is contained in this station 
IF(XLOCAL.GT.XLE(l)^CHORD(l)) THEN 
CALL WINGNWAKE(XLOCAL,KT,NUL,IFLAT} 

GOTO 200 
ENDIF 

FT - 0 

DO 160 I - 1,NUL 

XT - XT + 1 

Create a point at the root 
lF(XBASE(i, IFLAT, 1 ) .GT. XLOOAt) THEN 
XI - XBASE (1-1,1 FLAT, 1) 

X2 - XBASE (I , IFLAT, 1) 

Yl - YBASE ( I -1 , I FLAT , 1 ) 

Y2 - YBASE (1,1 FLAT , I ) 

21 - ZBASE <1-1, IFLAT, 1} 

22 - ZBASE ( 1 , IFLAT , 1 ) 

CALL LININT(X1 ,X2,Y1,Y2, XLOCAL, YY) 

CALL LININT(Xl ,X2 , Z1 , Z2 , XLOCAL, 22) 

YINT(KT) - YY 
ZINT(KT) - ZZ 
GOTO 200 
ENDIF 

DO 150 M - 2,NC 
XL - XLOCAL- XBASE ( I , IFLAT ,M-1 ) 

XR - XLOCAL-XBASE ( I , IFLAT ,M) 

I F ( XL * XR . LT . 0 . ) THEN 

XI - XBASE(I , IFLAT, M-l ) 

X2 = XBASE ( I , IFLAT ,M) 

Z1 - ZBASE ( I , IFLAT ,M-1 } 

22 = ZBASE (I, IFLAT, M) 

Yl =* YBASE { I , IFLAT , M- 1 ) 

Y2 = YBASE (I , IFLAT, M) 

CALL LININT(X1,X2,Y1,Y2, XLOCAL, YY) 

CALL LININT(X1,X2,Z1 ,Z2, XLOCAL, ZZ) 

YINT(KT) - YY 
ZINT(KT) - ZZ 
GOTO 160 

ElSEIF { ABS ( X L* XR ) . LF l.E-7) THEN 

I F( ABS ( XL) . LE .1 . E-7 ) THEN - 



DATE 


1 66B 

vm 

rm 


rm 

vm 


rm 




rm 

wr 

150 

rm 

riPfijjr 

160 

\ P if 

1674 

C 

1675 

200 

1 676 

c 

(677 

c 

1676 

c 

167$ 

c 

I PoO 
1661 
1662 

220 

vm 

(664 

1665 

230 

*666 

1667 

1666 

240 

166$ 
16$0 
16$ 1 

250 

I6$2 

C 

16$ 3 

C OFT 

V6$4 

C 

C 

1695 

■mi 

C 

■ml 

t 

16$$ 

1760 

1761 

c 

1762 

(763 



YINT(KT) 
ZINT(XT) 
ELSE 

YINT(KT) 
ZINT(KT) 
ENDIF 
GOTO 160 
ENDIF 
CONTINUE 
CONTINUE 
ENDIF 


YBASE (1,1 FLAT , M } 
ZBASE ( I , I FLAT , M) 


CONTINUE 

R ight H6%f tINT apd ZINT is fro* wing tip to tte root, i» order to 
tbe cdbic spline ptograa, need from root to tip. 

DO 220 XX- 1, XT 
S<XK) - YINT(XX) 

CONTINUE 
DO 230 XX-1,XT 

YINT(XX) - S<XT-XX+1) 

CONTINUE 

do 240 xx«i,xt 

S<XK) - ZINT(XX) 

CONTINUE 
DO 250 XX«1,XT 

ZINT(XX) - S(XT-XX+1) 

CONTINUE 


Distribute t coordinates < cap# around direction) froa ttk>t to 
leading edge and back (xdist) . 

I M , grid points Vill custer near the wing tip, -0 iveaxtb* root . 
For HESS, IT*0 i for Boeing, IT-1 
iT-'i 

CALL DISTARC ( ZINT, YINT, XT, ZDIST, YNEW, XTIP , FAC, IT) 

I F( L.NE.NSEC ) CALL CSPLINE ( ZINT, YINT, XT , ZDIST , YNEN , XTIP ) 

RETURN 

END 







SOURCE PROGRAM 

samgrid.f 


7/07/94 I ?ACE# 


time 4:18:56 pm 



I LINE # 


1 704 

1705 

1706 C 

170 7 C 
1/08 C 

1 709 c 

1710 
171 I C 
1712 C 
1715 

I 7 i 4 

1715 

1716 

1717 

1 718 

1719 2 

1 720 3 

1721 

1722 
1725 


SOURCE TEXT 


hommmomm * * 4* * * ;*• * ***** ♦ *4 * *:*;** * v * * **** 7. 

SUBROUTINE SUBTRACGRI D ( NPL1 , X,Y, Z,NPI /NSEC , XDIN) 

This subroutine allows us to substiact a grid Hue NP11 in 
stxtumriie section, and the new dimension la NSEC again. 

DIMENSION X(NPI),Y(NPI,NPt),Z<NPI,NPt) 

Ranuaiber the late stations 
NSEC * NSEC-1 
DO 30 1- NPL1,N5EC 
X(L) - X(L+1) 

DO 20 K-1,KDIM 
Y<L,I) - Y<L41,K) 

Z(L,K) - Z ( L+ 1 , K } 

CONTINUE 

CONTINUE 





SOURCE PROGRAM 

samgrid.f 


DATE 7/07/94 PACE # 

"™* 4:18:56 pm | 27 


SOURCE TEXT 


* f * ****** *** ** * **** * *****>**♦** *.*;**;* ** * ;* f ** *. *****;** ;+**;* + *'*+'<• ********* * + + 

SUBROUTINE TAI L ( FAC1 , FAC 2 ) 
include " sgrid . coi" 

DIMENSION S (NPI ) 

COMMON /REF/ ZROOT , KTIP , ARCORR 

ITIP = ( KDIH+1 )/2 

LE = L-l .. . 

I C- The tali of the configuration 
C HAXEPT i« the Y value that the wake i.A ( ie T value of ZROOT) 


PT1- < YOUT (NSEC , i j +Y0DT ( NSEC , NPTS ) ) /2 . 

PT2- <YIN(NF,NFP)+YIN(NF,1) )/2 . 

XI - X<LE) 

X2 - XIN<NF) 

LI - LE* 1 
L2 - L1+<NF-M2)-1 
DO 300 L-L1,L2 
M-K2+L-L1+1 
X { L ) - XIN<K) 

X LOCAL * XIN<M) 

C OPTIONS 

c NAK*PT - (riN{NrVNW ) mN<N^,l > )/2 . 

CXtL LlNINf ( XI , X2 , PT1 , PT2 , XtOCAL , WAlEPT ) _ 

C fiealgn the let point (IF) will be the off fuselage aide 

q U H *lao th* hnaljier of pt* in the uppe* o* lower fuselage, therefore 

C it depends bu the pravioua station . 

C {no. of pt In thin atation )/ ( no . of pta In' previous) 

C (radlua of tfcU station) /(radius of preVibua) 

ixf icb XBS (TtN < M , I j -YIN ( M , NFP ) V/XBS < YIN -YIN {M- 1 , NfP > ) 

IF- ( IDIM-NPTS )/2+l 
IF- ( IDIM-NPTS )/2 

♦ ir<*ATIO.LE.O;9 .OR. RATIO. GE.I.l) THEN 

* If-IFlt<RXflO*nOXT<lF)) 

* . . ENDIF "• ■ • : : - ! i* L^:.= .. 


find Xfl the point at old grid where YiNfM/mj-NxitePt 
and calculate the pointa fro* neg Y to Y at XT1 
DO 222 l-l , NFP 

IF(ABS(YIN(M,I)-HAIEPT) .LE.l.E-7) THEN 
IT1-K 
GOTO 223 
ENDIF 

IF(YIN(M,I) .GT.WAIEPT) THEN 
IT1-K 

YI-YIN { N , I- 1 ) 

Y2-YIN(M,I) 

Z1-ZIN{M,I-1) 

Z2«ZIN(H,I) 

YY-WAIEPT 

CALL LININT(Y1 ,Y2, Zl, Z1 , YY, ZZ) 

YIN(M,IT1 )*YY 
ZIN(M,IT1)-ZZ 
GOTO 223 
ENDIF 
CONTINUE 
CONTINUE 

DO 230 I— 1 , IT1 
YINT ( I ) -YIN ( M , I ) 

ZINT(I)-ZIN{M,I) 

CONTINUE 

CALL DISTARC( YINT, ZINT , KT1 ,YNEW, ZDIST, KF , FAC1 ,1 ) 

DO 240 I-1,IF 
Y ( L , I ) -YNEW ( I ) 

Z{L,I)-ZDIST(I) 

CONTINUE 


find KT3 the point At old grid where YIN<M,rra)-NXfEPT 
note t IT2*NFP-ITl+i 

and calculate the points fro® neg Y to Y at IT1 

iT2«NFP-m+l 

DO 250 1-1 , IT2 

YINT { I ) -YIN < M , IT1+ 1- 1 ) 

ZINT< I ) -ZIN(H, IT1+K-1 ) 

CONTINUE 

CALL DISTARC* YINT, ZINT, IT2 , YNEW, ZDIST , IF , FAC1 , 0 ) 

DO 260 1-1 , IF 

Y < L, IDIM-IF+I ) -YNEW< K ) 

Z< L, IDIM-IF+I )-ZDIST( I) 

CONTINUE 


Theae are pdinta in off fuselage aide 
KN-No. of pts in the off- fuselage side 
IN-IDIM-2 * IF 
INH- ( KN+1 )/2 

RSPAN- AB5( Z ( LE , KTI P ) - Z ( L , IF ) ) 

IF( FAC2 . LT . 0 . )FAC2-1.£+15 
DELT-FAC2 * ( RSPAN J/FLOAT ( INH ) 

CALL DI5TRI (DELT, KNH+1 , S, 0) 

CALL LININT{X1, X2,YOUT{ NSEC, NPTS/2 ) , PT2 , XLOCAL , YOB) 
DO 270 I- 1 , INH 

Z< L,IF+I)-Z<L,KF}+S( K + l )*RSPAN 

CALL LININT< ZROOT, RSPAN , WAIEPT , YOB , Z { L , IF+K ) , YY ) 
Y{ L, IF+I) -YY 
CONTINUE 
DO 2 B0 I— 1 , INH 

Z ( L, IDIM-IF-I+ 1 )-Z( L, IF+I) 

Y ( L, IDIK-IF-I+ 1 ) -Y ( L , IF+I ) 

CONTINUE 


Smooth the fuselage-wate part 
DO 285 IK=1 , KDIM 
YINT { IK )=Y ( L, II ) 

ZINT (KK)-Z(L,IK) 

CONTINUE 

ZROOT = ZIN (M , IT2 ) 

RFILT * RFIL 

IF(X(L) .GE. XSTART .AND X ( L ) . LE . XOFF } 

CALL FI LET (YINT, ZINT , KDIM , MD1 , MB2 , MT1 , MT2 , RFI LT) 














SOURCE PROCRAM 

samgrid.f 


4 : 18:56 pm 


SOURCE TEXT 


I F < TOTLEN . EQ . 0 . ) GOTO 510 
DELT“ (X { L2 )-X( L2-1 ) )/TOTLEN 
CALI DISTRI(DELT,NNAKE4l,S,0) 

DO 500 LD-1,NWA*E 
• L - L24LD 

X (L) « X ( L2 ) +S < LD+ 1 ) * TOTLEN 
DO 480 X-1,XDIM 
Z(L, I) - Z(L-1,K) 

Y(L, K) - Y(L-1,K) 

4 80 CONTINUE 

500 CONTINUE 
510 CONTINUE 

ccx:cccctxi«:cccc«^cccccccx:ccccaic<xcccccxccc<x^cccc<xccccccccccccccc 

c ; • . • i •'/ <: i :• ; ; x. . • 

c 

c ' ' x ••• •" v . ; ' ‘ ‘ : 

c He aight like fcp icaie the whole Udog by « reference length 

c - •. . 1 : • ■■■ : • 

NSECT-L2+NWAXE 
DO 550 L-1,NSECT 
X (L )*X( LJ/5CALE 
DO 550 K«1,XDIM 

Y ( L/T)*Y ( L, K J/SCALE 
Z<L,I)-Z(L,R)/SCALE 
550 CONTINUE 


C OPTIONS . 

c toeing Bhbpiifte CoefiguratiPtt . 

C Before writ* out the grid, tbie boei**9 **•*&*& to be fixed 

C ' ' fOr edee stupid rwahoa*, : 

o DO 570 

o *<ia;*>-mi,x) 

C* 

c570 COHTIHOR ' = 

0 Cfsii SU*T*AC6BID( 31 ,1,*, %<MPt , N«CT , &tH ) 

c Mt r *0 , X, V , Z, »«,#***, 





SOURCE PROCRAM 

samgrid.f 


7/07/94 PACE * 


time 4:18:56 pm 



SOURCE TEXT 


• * ** ** 444 **** *44 4 #** * 44 ' 4'**;*4 * * 4 ‘ 44 *** 4’*4 4'*‘*;*4 * 4 * 44 *** « 4*4 * * #4 *> * 4.4 *>♦.<*; <*‘*4 **44 
SUBROUTINE WFHATCH 
Include "sgrld.ccm" 

DIMENSION Dl(NPI) ,D2<NPI ) 

COMMON /REF/ ZROOT , KTI P , ARCOFR 

C ■ 7 \ • •' . . v.;! 1 -••• " ' • 

;• q' ' This subroutine •ovee the whole inward or outward such that the 
C fuselage and the wing^iOOt iaetcjh. 


The wing “body section, itself lias NSEC sections 
LW is the Section Index fox Twr^ZIN 
DZMAX - 0. 

DZNIN - l.E+20 

LB - Ml+1 

LE - K1+ (NSEC-1 ) 

DO 300 L« LB, LE 
LW-L-K1+1 

X LOCAL - XOUT ( LW , 1 ) 


Calculate the points (YI NT, ZXNT) on the fuselage at XLOCAL 
Note: assumed that In this section, each station has save ntmber 
-x : of points ia rai^axdttnd direction. 

DO 10 M-Ml/NF 

IF{XIN{M) -GE.XLOCAL .OR. 

ABS(XIN(M)-XLOCAL) .LE.l .E-7) THEN 
MW - M 
GOTO 15 
ENDIF 
CONTINUE 
CONTINUE 
DO 30 K«1,NFP 
Xl-XIN(MW-l) 

X2«XIN(MN) 

Y1-YIN(KW-1,I) 

Y2-YIN<KN,K) 

Zl-ZIN<KW-l,I) 

Z2-ZIN<MW,X) 

CALL LININT( XI ,X2 , Y1 ,T2 , XLOCAL, YY ) 

CALL LININT< XI, X2,Z1,Z2, XLOCAL, ZZ) 

ZINT(X) - ZZ 
YINT(K) - YY 
CONTINUE 
CONTINUE 


Hot# the wing in the redirection < spanwise j to make sure 
the wing is not inside of betside the fuselage 


the wing is not inside ot cetside the fu 

ZROOT - ZINT( (NFP+1 )/2) 

DO KZ1 - 1,NFP 

IF(ZINT(EZ1) .GT. ZROOT} ZROOT- ZINT< IZ1 ) 
ENDDO 

DZ - ZOUT(LW,l) -ZROOT 
IF{ ABS(DZ) . GT. ABS(DZMAX) ) DXHAX-DZ 
IF{ABS{DZ) . LT . ABS ( DZMIN ) ) DZMIN-DZ 
CONTINUE 

IFfDZMAX . GT. 0 . ) DZ=DZMIN 
IF ( DZMAX . LT . 0 . ) DZ-DZHAX 
DO 400 K«1,NPTS 
DO 400 L-l , NSEC 

ZOUT ( L , K ) - ZOUT<L,X) - DZ 
CONTINUE 

wrlte( * , * ) 'DZ - ' , DZ 
RETURN 
END 






SOURCE PROGRAM 

samgrid.f 


SUBROUTINE WINg/boDY 
include "sgrid. coi" 

DIMENSION Dl(NPI) ,D2(NPI) 
COMMON /REF/ ZROOT , XTIP , ARCORF 


SOURCE TEXT 


, * * ***** ******* * ** A********;*** * 


The wing-body section/ Itself be a NSEC sect lone 
LW is the section Index for YIN , ZIN 

LB - ki+i 

LE - Ml + NSEC -1 
DO 200 L-LB ; LE 

LW-L-M1+1 (start with second wing station 
X LOCAL - XOUT(LW,l) 


Calculate the points <YINT,ZINT) on the fuselage at XLOCAL 
Note i assueed that ia this section/ each station has sane nvaaber 
of point* in rape- around direction. 

DO 10 H-Ml ,NP 

IF(XIN(M) . GE . XLOCAL .OF. 

ABS ( XIN(M) -XLOCAL) . LE . 1 . E-7 ) THEN 
MW - M 
GOTO 15 
ENDIF 
CONTINUE 
CONTINUE 
DO 30 X«1,NFP 
Xl-XIN(MW-l) 

X2-XIN(MN) 

Y1-YIN{MW-1,I) 

Y2-YIN(HW,X} 

Z1-ZIN(MH-1,X) 

Z2-ZIN<MW,K) 

CALL LININT< XI ,X2,Y1,Y2,X LOCAL, YY> 

CALL LININT( XI ,X2,Z1,Z2, XLOCAL, ZZ) 

ZINT{X) - ZZ 
YINT(I) - YY 
CONTINUE 
CONTINUE 

There are NPTS points in the wing ares, hei<$ to check out 
how many point la in the fuselage 

XI is the inters pt ./. the fuselage a wing at hot tew in old grid 

*2 ia the Inters pt ./. the fuselage a wing at top in old grid 

Hit) is the inter point ./. the fuselage a wing at bottom in ft** grid 

There are IT points iron 1 to |^1 

XT- (NFP+1 )/2 

MID- < KDIM-NPTS+2 >/2 

MID- ( IDIM-NPTS )/2 

rind ' Xl':i ' : i 

DO € 0 K-XT,1,-1 

I F ( YINT ( X } . LT . YOUT ( LW , 1 ) .AND. YINT< I ) .NE . YINT( 1-1 ) ) THEN 

11 - I 
GOTO 62 

ENDIF 

CONTINUE 

CONTINUE 

rind n 

DO ‘JO I-l,NtP 

IF(YINT<I} . GT . YOUT( LW , NPTS } . AND.YINT(X) . NE . YINT { 1+1 ) ) THEN 

12 - I 
GOTO 72 

ENDIF 

CONTINUE 

CONTINUE 

Fgr arrow-wing type, f? seeded to be relocated 
IF(ABS(YOUT<LH,1)-YOOT(LW,NPTS) ) .LE. l.E-7) THEN 
XI - II + XI ADD 
12 - II 

YINT <X2) -YINT (XI) 

ZINT(X2)-ZINT(K1 ) 

ENDIF 


Calculate the points at the bottow ( from a eg. Y to Y at MID) 
X$»l 
XE- II 
XN-IE-IS+1 
DO 82 I- KS , XE 
XX - X-XS+1 
D1 ( XX) - YINT{ X ) 

D2{ XX) - ZINT< X } 

CONTINUE 

CALL DI5TARC ( D1 , D2 , XN , YNEW, ZDIST , MID , - 10 . , 1 ) 

DO 100 X-l , MID 
Y(L,XJ - YNEN{X) 

Z(L,I) - ZDIST(I) 

X (L ) - XLOCAL 
CONTINUE 

The points of the wing section 

DO 110 X-1,NPTS 

Y<L,HID+X) - YOUT(LW f X) 

Z<L,MID+K) - ZOUT<LW,X) 

X<L) - X0UT( LW, X ) 

CONTINUE 

Calculate the points at the top 
KS-X2 
XE- NFP 
KN-XE XSU 
DO 115 K> IS , XE 
XX = X-KS+1 
Dl(XX) - YINT ( X } 

D2(KK) * ZINT{ X | 

CONTINUE 

CALL DI5TARC{ D1 , D2 , KN , YNEW, ZDI ST , MID ,-10 . ,0) 

DO 120 K-l , MID 

Y { L , K+MID+NPT5 ) = TNEW(K) 

Z{ L , X+MIIHNPTS ) - ZDIST ( X ) 

X ( L ) - XLOCAL 
CONTINUE 


For arrow-wing type. Mate sure wafce points oh 
IF( ABS (YOUT{LW,1)-YOUT(LH, NPTS) ).LE.l.E-7) THEN 
Y(L,MID+NPTS+1) = YOUT ( LW , 1 ) 

Y { L f MID ) = Y ( L , MI D*NPTS+ 1 ) 






SOURCE PROCRAM 

DATE 7/07/94 

samgrid.f 

time 4:18:56 pm 



SOURCE TEXT 


Z ( L , KID} 
ENDIF 


Z(L,MID+NPTS+1> 


Fill the uvnootb part by FILET 

First of all, find the set of points needed to be rearrange 
DO 125 Kk-l,KtHM 
YINT(KK)“Y(L,KIC) 

ZINT(KK)-Z(L,riC) 

CONTINUE 
RFILT - RFIL 

IF ( X( L) . GE . XSTART . AND. X(L) . LE.XOFF) 
a CALL FILET(YINT, ZINT, KDIM,MB1 ,MB2 , KT1 ,MT2 , RFILT) 

DO 140 K-l/KDIM 
Y(L,K)-YINT(K) 

Z(L,X)«ZINT(K) 

CONTINUE 

CONTINUE 

RETURN 

END 






SOURCE TEXT 


«44*»***«*4*****4* **> * * •** '**;#• '* +>**:'• *** *!*.* **4*;$*** +.'*. ♦t.*!*!* j^****!*:*** 

'SUBROUTINE wiNGNMAlii^ FLAT ) 

Include "sgrid.coa* 

This subroutine creates the data when the station is in the place 
where some points are on the wing, soae are the wake >•• 

At that x- station, we put 10 point* in the wake part. 


XT - 0. 

I F ( ARRWING . GT . 0 . ) THEN 

The wine is sveeped backwards 
DO 75 J-2,NC 
Jl-J-1 
J2-J 

IF(XLOCAL.GT.XLE(J)+CHORD(J ) ) GOTO 75 
XIl-XLE(Jl)+CHORD{Jl) 

Zl -ZBASE (NUL, I FLAT ,J1) 

XI2«XLE(J2)4CHORD(J2) 

Z2-ZBASE(NUL,IFLAT,J2) 

ZTIP-Z1+ ( (XLOCAL-XI1 )/(XI2-XIl) )*(Z2-Z1) 
ZFLAT-ZTIP 


NLEDG - NC 
DO KRUN-2,NC 

IF < XLOCAL. LT. XBASE (1,1 FLAT, KRUN)) THEN 

If X LOCAL < leading edge, we ahould do the following 

NLEDG - MR UN-1 

XI - XBASE (1,1 FLAT , HP UN- 1 ) 

X2 - XBASE (1,1 FLAT, MRUN ) 

Y1 - YBASE{ 1,1 FLAT, MRUN-1) 

Y2 - YBASE (1,1 FLAT , MR UN ) 

Z1 - ZBASE (1,1 FLAT, HRUN-1) 

Z2 - ZBASE (1,1 FLAT, MRCN ) 

CALL LININTf XI, X2,Y1,Y2, XLOCAL, YY) 

CALL LININT(X1 , X2 , Z1 r Z2 , X LOCAL , ZZ ) 

YINT(l) - YY 
ZINT{1) - ZZ 
GOTO 40 
ENDIF 
ENDDO 
CONTINUE 


DO 65 M-NLFDG , J2 , -1 
DO 64 1-2, NUL 

IF( XLOCAL . LE . XBASE ( I , IFLAT,M ) ) THEN 
XI - XBASE (I-l,irLAT,M) 

X2 - XBASE ( I ,IFLAT,M) 

Y1 - YBASE(I-1,IFLAT,M) 

Y2 “ YBASE ( I , I FLAT , M ) 

Z1 - ZBA5E(I-1 ,IFLAT,M) 

Z2 - ZBASE ( I , I FLAT , M ) 

CALL LININT( XI, X2,Y1,Y2, XLOCAL, YY) 
CALL LININT(X1,X2,Z1,Z2, XLOCAL, ZZ) 
ir< NLEDG. EQ.NC} THEN 
MM-NLEDG-M+l 
ELSE 

MM-NLEDG-M+2 

ENDIF 

YINT(MM) - YY 
ZINT(HH) - ZZ 
GOTO 65 
ENDIF 
CONTINUE 
CONTINUE 


Add 10 points for the wake. 

ZROOT - ZBASE ( NUL , I FLAT , L ) 

Y1 -YBASE < NUL , I FLAT , 3 1 ) 

Y2 » YBASE (NUL, I FLAT, J 2 ) 

Zl- ZBASE (NUL , I FLAT , J 1 ) 

Z2 - ZBASE ( NUL , I FLAT, J 2 ) 

CALL LININT(XI1,XI2,Y1,Y2, XLOCAL, YY) 

CALL LININT(XI1 ,XI2 ,Z1 , 22 , XLOCAL, ZZ ) 

DZ - ( ZFLAT— ZROOT )/9 . 

DO 70 M-1,10 

IF (NLEDG. EQ.NC) THEN 
MH - NLEDCJ2+1+M 
ELSE 

KK - NLEDG-J2+14M +1 
ENDIF 

ZINT(MM) - ZROOT 4 FLOAT ( 10-M) *DZ 
CALL LININT( ZROOT , ZZ, YOUT( L-l , 1 ) , YY , ZINT(HM) , YYY } 
YINT(MM) - YYY 
CONTINUE 

IF(NLEDG. EQ.NC) THEN 
XT - ( NLEDG- J 2 + 1)410 
ELSE 

XT - (NLEDC-J2+1 +1)+10 
ENDIF 
GOTO 200 
CONTINUE 


The wing is sweeped forward* 

DO 105 J-2,NC 
Tl-J-1 
J2-J 

IF (XLOCAL . LT . XLE ( J ) +CBORD( J ) ) GOTO 105 
XI 1-XLE ( J1 ) 4CH0RD( J1 ) 

Zl- ZBASE ( NUL , I FLAT , J1 ) 

XI 2*XLE( J2 ) +CHORD( J2 ) 

Z 2- ZBASE (NUL, I FLAT, J 2 ) 

ZTIP- Z1+ ( (XL0CAL-XI1 )/ ( XI2-XI1 > ) * (Z2-Z1 ) 
ZFLAT-ZTIP 

-SIMP c 2YTP-2BXSE(NUL,IFLAT,NC) 

ZTIP= ZB ASE ( NUL , I FLAT , NC ) 


Add 10 points for the wake. 

Y1-YBASE(NUL,IFLAT,J1) 

Y2= c YBASE(NUL, IFLAT, J2 ) 

Zl = ZB ASE (NUL, I FLAT, Jl) 

Z2 - ZBASE ( NUL , I FLAT , J2 ) 

CALL LININT(XI 1 , XI 2 , Y1 , Y2 , XLOCAL , YY ) 

CALL LININT<XI1,XI2,ZI,Z2, XLOCAL, ZZ) 

DZ - (ZTIP-ZFLAT)/?. 

DO 80 M-1,10 

ZINT(M) - ZTIP - FLOAT(M-l )‘DZ 

CALL L1NINT(ZTIP,ZZ,YQUT(L l,NPTS/2+l) ,YY,ZINT(M) , YYY) 








IIISP1 

" ' ? j 


SOURCE PROCRAM 

samgrid.f 


s* 




DATE 

7/07/94 

TIME 

wmdm 

4:18:56 pm 

'Si‘5St 


PACE # 


35 


2392 ] 

m\ 

2T9+ 

2395 

25?: , 
2598 

zm\ 


24 n 


YYY 


YINT(M) 

CONTINUE 
DO 100 M-Jl ,1,-1 
DO 90 1-2, NUL 

I F ( X LOCAL . LE . XBASE (1,1 FLAT ,M ) > THEN 
XI - XBASE ( 1-1,1 FLAT ,M> 

X2 - XBASE (I , IFLAT ,M) 

Y1 - YBASE (1-1,1 FLAT, M) 

Y2 - YBASE (I , IFLAT, H) 

Z1 - ZBASE (1-1,1 FLAT , H J 
Z2 - ZBASE { I , IFLAT, H) 

CALL LININT<X1,X2,Y1,Y2,XL0CAL,YY) 
CALL LININT (X1,X2,Z1,Z2, XIOCAL, ZZ > 
YINT<J1-M+1 +10) - YY 
ZINT( Jl-M+1 +10) - ZZ 
GOTO 100 
ENDIF 
CONTINUE 
CONTINUE 
XT - 10 + J1 
GOT0200 
CONTINUE 
ENDIF 


SOURCE TEXT 


IS 


CONTINUE 

RETURN 

END 




SOURCE PROCRAM 

samgrid.f 


DATE - 


7/0//94 


time 4:18:56 pm 


36 




LINE# 


SOURCE TEXT 


c* v****** * « * * m******** * * **'**?&*?*'* **&M*'** * ******?■? ■■■ ■■ 

SUBROUTINE hingin' 
include "sgrid co»" 

C NC “ t of sect ion* in the sprawls* direction 

C 4BASE - Z value of Itb section . . . : 

c xtt(K) ^ Xending edge X value of eectidn ;. ; ;• 

C ytE(K) « leedlng edge T of Ktl> eectidft 

C GfiORD(K) - Chord length of the fth nection B _ 

C NO « | of points in the upper * lower (Kth) section 

t Note ; the end points of upper end X $m seetiooe are s«*e physical P ts 

c;'-; " ’_• 

C MACH2 configuration 
READ( 10,900) 

READ( 10,910) 

READ< 10 , 920 )NC,NU 
NL - NU 

K-l 

111 CONTINUE 

REXD{ 10,930) 

READ( 10,950) 

RE AD < 10 1 * ) XBASE ( 1 ,2,1), YBASE (1,2,1), IBASE ( 1 , 2 , E ) 

READ( 10,940) 

FEAD( 10 ,950) 

READ(10,MXBASE(1 ,1 ,K) , YBASE (1, 1,*) , XBASE< I , 1 , K ) 

XLE(K) - XBASE (1,1,1) 

YLEIX) - YBASE ( 1 , 1 , X ) 

CHORD (I) - ABS(XBASE(1 ,1 ,K) -XBASE (NU,1,K)) 

F-IC+1 

IF(X. LENC)GOT0111 


F6rMAT(1X) 

FORMAT (IX) 

FORMAT ( 25X, 1 5 , 9X , 1 5/ ) 
FORMAT (IX) 

FORMAT (IX) 

FORKAT(IX) 

FORMAT ( 3F16 . 7 ) 

RETURN 

END 




SOURCE PROCRAM 

samgrid.f 


7/07/94 I PACE# 
4:18:56 pm | 



SOURCE TEXT 


SUBROUTINE FUSEIN 
include "sgrid. com" 
DIMENSION YT(NPI) , ZT(NPI) 


| of section In the fuselage 
| of points in ith section 


Read ' the fuselage geometry 

MACK? configuration 
READ{ 10, 810} 

READ< 10,820) 

READ j 10, 830)NF,NFP 

M-0 

CONTINUE 

READ (10 ,840) 

M « M + 1 

DO 100 X-1,NFP 

READ(10, *) XIN(M) , YIN(M, K) , ZIN(M, X ) 

YT(K)«YIN(M,K) 

ZT(K)-ZIN(M,X} 

CONTINUE 

CAtli DIBTA*C(VT, ET,NFT , YT, ZT ,NIT>10 . , 0 ) 

DO 110 X-1,NFP 
YIN<N,K)-YT(K) 

ZIN(M,X)-ZT<K) 

CONTINUE 

IF(M.LT.NF) GOTO 40 

write this ti»ael4gie g<k^try into Rt6rit> PI nnhr foomat 

Kw-i 

KK-NFP 

WRITE ( 51 )KK, XW,M 
DO 800 L»1 ,M 

NRITE{ 51) (XIN(L) , K*1,KK) , 

(yin(l,x) ,K- i,rr) , 

(ZIN(L,X) ,1C«1,KK) 

CONTINUE _ v- 

FORMAT ( IX ) 

FORNAT(IX) 

FORMAT( 31X ,15,81,15/} 

FORMAT (IX) 

FORMAT (2X, 3F1€ . 7) 

RETURN 

END 








SOURCE PROGRAM 

samgrid.f 


SOURCE TEXT 


SUBROUTINE NACIN 
include "sgrid.con" 

DIMENSION YT(NPI ),ZT<NPI) 

COMMON /NACG/ XNAC(NPI ) , YNAC(NPI ,NPI J , ZNAC(NPI ,NPI > 
COMMON /NDIM/ NNAC ,NNACP , LNUM( 4 ) 


« # or section io the nacelle 
« f of points in nth section 


Read the lUceile geometry 
READ< 30,810) 

READ< 30,820) 

RE AD j 3 0, 8 3 0 ) NN AC , NN ACP 

k»o 

CONTINUE 


READ ( 30,840) 

M * M + 1 / 

do 106 k-i,nnacp 

READ( 30,8 50) 1NAC<M) ,YNAC (M , X ) , ZNAC<H, X) 
YT(X)-YNAC<H,K) 

ZT(X)-ZNAC(M,X) 

CONTINUE 

CXll DISTAHCtrr, XT,tWACP,YT , Zf , NNACP ,-10. ,0) 
DO 110 E*1 • . 

W<K, E)-Ef(X) 

CONTINUE : : - - . ; ' . . 

IF(M . LT.NNAC) GOTO 40 


NJfltis the nacelle geometry Into PfcoflD Planar fomat 

. ■ kw-j - i;. 

KX»NNACP . 

DO *00 D-l 

u i»xc< tj ; ; o j , 

^YMXGtL, K) V*-! , EE), 

• - , ••• t uiac<i , it) , r»i, i ty - 

■ continue ' 

call flush j 

FORMAT < IX) 

FORMAT < IX) 

FORMATS 3IX, I 5,81,15/) 

FORMAT (IX) 

FORMAT( 2X, 3FIS .7 ) 






SOURCE PROCRAM 

samgrid.f 


LINE # 


2568 
2W c 
2576 
25! 71 c 
2572 c 

2575 c 
2574” c 
257$ c 

2576 c 

2 578 c 
2575 c 

2586 c 

2581 c 

2582 c 

2583 c~~ 

2584 c 

2585 C 
2566 

2587 

2588 
2585 
2596 c 

2591 c 

2592 c 

if 

2596 


SOURCE TEXT 


SUBROUTINE KINGMAKER 
include "sgr id. cob" 

This prdgraa generates a • flipped 1 delta wing with no iwiit 
on airfoil coordintea read in fro*! tort. 90. ¥o t»e as 
■part: of saegrid, RlNGlN 'is not necessary < oor is ytaiSwRfrj.''’ 

Nritten by; Donovan L. Mathias 
July 1992 

whenever possible* the same variables are used as lo samgridf 

fort. JD V airfoil coordinatres 

fort. 77 deacrltion of the wing 

: Declarations 

real xte< ii) , sidpei , si 6 PE 2 , scale ,xaf< i so j 

REAL YU<150) ,YL( 150) 
real SPAN 
INTEGER L, I , J , X 

Initialization 

' read ( 77 , * )NC 
read (77 , * )XLE( 1 ) , XLE (NC ) 
read(77,*)XTE(l),XTE(NC) 
retd (77 , * ) ZBASE (1,1,1) , ZBASEfl, 1 ,NC) 

Read in ai foil coordinates (t la I of X coozds.) 

READ( 90 , * ) :: 

READ(90,*) 

READ(90,*)NU 
READ< 90 , * ) 

NL-NU 
DO I-1,NU 

READ (90,19 )XAF( I ) , YU ( I ) ,YL(I ) 

END DO 

format (lx , f 9 . 7 , 3x, f 9 . 7, Jx, f 9 . 7 ) 

Establish z distance (Sparod.**) 

SPAN • ZBASE ( 1 , 1 ,NC )- ZBASE (1,1/1) 

bo E-6,^C-i ’ : : : 'C yy •: ?:-:W 

DO I-1,NU 

ZBASE(I,1,K+1) - ( K* (SPAN/ (NC-1 ) ) ) 

ZBASE(I,2,X41) - ( K * ( SPAN/ ( NC-1 ) ) ) 

ENDDO 

ENDDO 

Establish $we«p (1 PGB LE, 3 FOR TE) 

SlbPEi - ( XLE(NC ) -XLE ( 1 ) )>( ZBASE ( 1 , 1 , NC )- ZBASE (1 X ll ) ' 

SL0PE2 ", .<.™ < NC ) " XTE ( 1 ) ) /( ZBASE (1,1, NC ) - ZBASE ( 1,1,1)) 

Generate leading' and trailing edge a 

" DO X- 1 ,NC X ‘ : ‘ X ' 

XLE(X) - XLE(l) 4 SLOPE1 * ZBASE (1,1, K) 

XTE ( K ) - XTE ( 1 ) +SLOPE2 * ( ZBASE ( 1 , 1 , X ) -ZBASE { 1,1,1)) 

ENDDO 

..Distribute grid point* 

DO K«1,NC 

SCALE - XTE ( X)-XLE ( K ) 

DO I“1 ,NU 

XBASE ( I , 1 , X ) - XLE(X) 4 SCALE*XAF(I) 

YBASE(I / 1,K) - YU ( I ) ‘SCALE 

XBASE ( I , 2 , X ) - XLE(X) 4 SCAL£*XAF(I) 

YBASE ( I , 2 , K ) - YL( I ) ‘SCALE 
ENDDO 
ENDDO 

Return valueis ’ to 'original' name* 

DO F-l ,NC 

CHORD ( X ) - ABS (XLE (X) -XTE ( K) ) 

YLE(X) - YBASE ( 1 , 1 , X ) 

ENDDO 

RETURN 

END 
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Preface 


This manual serves as a supplementary document for the official refer- 
ence manual of a relatively new research software, PVM, which has been 
developed at Oak Ridge National Laboratory. A beginner, who has no 
previous experience with PVM, would find this manual useful. 

We would like to thank you in advance that if you find any problems in 
PVM or this manual, please contact one of us. 

Mr. Merritt Smith 

NASA Ames Research Center, MS 258-1 
Moffett Field, CA 94035 
e-mail: mhsmith@nas.nasa.gov 
phone: (415) 604-4462 


Dr. Samson Cheung 
MCAT Institute 

NASA Ames Research Center, MS 258-1 
Moffett Field, CA 94035 
e-mail: cheung@nas.nasa.gov 
phone: (415) 604-4462 



INTRODUCTION 


7 INTRODUCTION ■ 

This manual provides you with an introduction to PVM and provides the 
fundamentals necessary to write FORTRAN programs in the PVM envi- 
ronment through a tutorial sample. This manual is designed for those 
who have no previous experience with PVM. However, you should know 
basic FORTRAN programming and UNIX. If you are ready for an 
advanced PVM application, please consult the official PVM Reference 
Manual. 


Software 

Package 


PVM stands for Parallel Virtual Machine. It is a software package that 
allows a heterogeneous network of parallel and serial computers to 
appear as a single concurrent computational resource. PVM allows you 
to link up all or some of the computational systems on which you have 
accounts, to work as a single distributed-memory (parallel) machine. We 
call this a Virtual Machine. 


PVM is useful for the following reasons. Unlike large mainframe com- 
puters or vector supercomputers, workstations spend most of the time 
idle. The idle time on a workstation represents a significant computa- 
tional resource. PVM links these workstations up to become a powerful 
multi-processor computational machine. With PVM, the lack of super- 
computer resources should not be an obstacle to number crunching com- 
putational programs. Furthermore, the annual maintenance costs of a 
vector supercomputer is often sufficient to purchase the equivalent com- 
puting resource in the form of workstation CPU’s. 



Definitions Here are some terms we use throughout this document: 

Virtual Machine PVM links different user-defined computers together 
to perform as one large distributed- memory computer. 
We call this computer the Virtual Machine. 

Host Individual computer (member) in the virtual machine. 

Process Individual program operating on different computers 

or hosts. 

Processor The processing unit in computers. A virtual machine 

can be viewed as a multi-processor computer. 
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INTRODUCTION 


Task The unit of computation handled by the virtual 

machine. You may want to think of one processor han- 
dling one task. 

Tid Task identification number which is a unique number 

used by the daemon and other tasks. 

Console A program from which you can directly interact with 

the virtual machine. (Add hosts, kill processes,...) 


Structure of 
PVM 


The PVM software is composed of two parts. The first pail is a daemon. 
We call it pvmd. This is the control center of the virtual machine. It is 
responsible for starting processes, establishing links between processes, 
passing messages, and many other activities in PVM. Since the daemon 
runs in the background, you have to use PVM console to directly interact 
with the virtual machine. 

The second part of the system is a library of PVM interface routines 
located in nbpvm3 . a. This library contains user callable routines for 
message passing, spawning processes, coordinating tasks, and modifying 
your machine. In writing your application, you will need to call the rou- 
tines in this library. 


Directory 

Setup 


This setup is for NAS system. Before you use PVM, you need to set up 
the following directories on all the machines that you want PVM to link: 

• Make a directory $HOME/pvm3/bin/ARCH in all the hosts of the virtual 
machine. 



Note arch is used throughout this manual to represent the archi- 
tecture name that PVM uses for a given computer. The table 
in the Appendix lists all the ARCH names that PVM sup- 
ports. For example, for Silicon Graphic IRIS workstations, 
you should make a directory $HOME/pvm3/bin/sGi. 

Make a directory $H0ME/pvm3/inciude, and copy the file f pvm3 . h 
from /u sr/n as/pkg/pvn 3 . 2 /inciude . (If you are on different sys- 
tem from NAS, please consult your system consultant.) 

Make a directory $HOME/pvm 3/codes, and write your application 
programs in this directory. You can actually put your programs any- 
where you like as long as the correct “include” files are includes. 
The current setup is for clarity. 
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Programming Concept 
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Master and 
Slaves 


Programming Concept 



Unlike graphical software or a word-processor, you cannot see PVM 
working by clicking your mouse buttons. In fact, a virtual machine is 
quite an abstract concept because you don’t physically have a multi- 
processor machine! In this chapter, you will learn a simple concept, 
which will help you to visualize how PVM works. 


A common way to work with PVM is a Master/ Slave relationship. A 
Master process starts Slave routines and distributes work. However, a 
Master does not actively participate in the computation. A Master 
process most often resides on the originating host (user’s computer), 
while the Slave programs are distributed to the hosts of the virtual 
machine. 


You need to distribute executables of Slave programs to the directory 
$HOME/pvm3/bin/ARCH on every host. You can locate this Master pro- 
gram anywhere you like. 

Since the Master program spawns Slave programs on each of the 
hosts to do jobs, it is important to understand the communication 
(message passing) amoung the hosts in PVM. 

Topically, a Master and a Slave have the following logic: 

Master Slave 

1 Enroll itself to PVM 

2 Spawn slave processes 

3 Initialize buffer, pack, and send 
message to all slaves. 

4 ...wait for slaves to finish... 

5 Receive message from slave(s) 

6 Exit PVM 


1 Enroll itself to PVM 

2 Receive message from master 

3 ...do something useful... 

4 Initialize buffer, pack, and send 
message to master 

5 Exit PVM 


The figure on the opposite page graphically describes a Master/Slave 
relationship and shows the exchange of information. 




FIGURE 1 


'.LAVE 


in Master/Slave programs. 



5 







Programming Concept 


SPMD 


Another common way to work with PVM is the SPMD, Single Program 
Multiple Data model. There is only a single program, and there is no 
Master program directing the computation. The user starts the first copy 
of the program and using the routine pvmfparent( ), this copy can deter- 
mine that it was not spawned by PVM, and thus must be the first copy 
(parent). It then spawns multiple copies (children) of itself and passes 
them the array of tids. At this point each copy is equal and can work on 
its partition of the data in collaboration with the other processes. 

Typically, a SPMD program has the following logic: 

1. Enroll in pvm 

2. If I am the first copy (parent) 

a) Spawn child processes 

b) Initialize buffer, pack, and send message out 

3. If I am a secondary copy (child) 

Receive messages 

4. Work!... Work! ...Work! 

5. Exit PVM 

The program on the opposite page describes a SPMD logic and shows the 
exchange of information. Please spend some time to study the program. 

In the next chapter we will introduce the PVM daemon and the funda- 
mentals of message passing. 



Programming Concept 


SPMD 

Program 


c 

c SPMD Fortran example using PVM 3.0 

c 

program spmd 

include ' . . /include/f pvm3 . h 1 
PARAMETER ( NPROC-4 ) 
integer mytid, me, i 
integer tids ( 0 : NPROC) 


1 



Enroll in pvm 

call pvmfmytid( mytid ) 


c 

c Find out if I am parent or child 

c 

— — — — call pvmf parent (tids(O) ) 

if( tids(0) .It. 0 ) then 
tids(0) = mytid 
me - 0 

c 

c start up copies of myself 

c 

call pvmf spawn ( ’ spmd ’ , PVMDEFAULT, ' * 1 , 

* NPROC- 1 , tids ( 1 ) , info) 

c 

c multicast tids array to children 

c 

call pvmf in itsend ( PVMDEFAULT, info ) 

call pvmf pack ( INTEGER4 , tids, NPROC, 1, info ) 

call pvmfmcast{ NPROC-1, tids ( 1 ) , 0, info ) 

else 


c 

c receive the tids array and set me 

c 

_____ call pvmfrecv( tids(0), 0, info ) 

call pvmfunpack( INTEGER4 , tids, NPROC, 1 , in fo ) 
do 30 i-1, NPROC-1 

if( mytid . eq . tids(i) ) me - i 

30 continue 

endif 


4 

5 


c 

c all NPROC tasks are equal now 

c and can address each other by tids(O) thru tids (NPROC- 1 ) 
c for each process me => process number [ 0- (NPROC-1 ) ] 
c 



print*, T me -',me, 1 mytid mytid 

call dowork('me, tids, NPROC ) 


c 

c program finished exit pvm 

c 


call 

stop 

end 


pvmf exit ( info) 






Programming Concept 


subroutine dowork ( me, tids, nproc ) 
include ' . . /include/f pvm3 . h f 

c 

c Simple subroutine to pass a token around a ring 



integer me, nproc 
integer tids( 0: nproc) 

integer token, dest, count, stride, msgtag 

count - 1 
stride - 1 
msgtag ■ 4 

if( me .eq. 0 ) then 
token * tids(0) 

call pvmf initsend ( PVMDEFAULT , info ) 
call pvmfpack( INTEGER4 , token , count , stride , in fo ) 
call pvmfsend( tids(me+l), msgtag, info ) 
call pvmfrecv( tids (nproc- 1 ) , msgtag, info ) 
print*, ’token ring done 1 
else 

call pvmfrecv( tids(me-l), msgtag, info ) 

call pvmf unpack ( INTEGER4 , token , count , stride , in fo ) 

call pvmf initsend ( PVMDEFAULT, info ) 

call pvmfpack( INTEGER4 , token , count, stride, in fo ) 

dest - tids(me+l) 

if( me .eg. nproc-1 ) dest - tids(O) 
call pvmfsend( dest, msgtag, info ) 
endif 


return 

end 
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Console 

Commands 


PVM Daemon 



The PVM daemon is the control center of the virtual machine. You can 
activate the PVM daemon by starting the PVM console or by invoking 
the daemon directly with a list of hosts. The latter will be discussed in 
chapter 6. To start the console, enter pvm at UNIX prompt on your local 
machine. The PVM console prints the prompt 

pvm> 

and accepts commands from standard input. The console allows interac- 
tive adding and deleting of hosts to the virtual machine as well as interac- 
tive starting and killing of PVM processes. Even if the daemon is started 
directly, the console can be used to modify the virtual machine. 


Here are the commands available in the PVM console: 


ADD 

ALIAS 

CONF 

DELETE 

ECHO 

HALT 

HELP 

ID 

JOBS 

KILL 

- MSTAT 
PS 

PSTAT 

QUIT 

RESET 

SETENV 

SIG 

SPAWN 
UN ALIAS 
VERSION 


add other computers (hosts) to PVM 
define and list command aliases you set 
show members in virtual machine 
remove hosts from pvm 
echo arguments 

stop all pvm processes and exit deamon 

print this information 

print console task identity 

display list of running jobs 

terminate tasks 

show status of hosts 

list tasks 

show status of tasks 

exit PVM console, but PVM daemon is still activated 
kill all tasks 

display or set UNIX environment variables 
send signal to task . 
spawn task 

remove alias commands you previous set 
show PVM version 
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Console 

Usage 


Suppose the console is running on workstation win210. This computer 
will automatically be a host in your virtual machine. Here are some 
examples of using PVM console: 


1. Activate PVM console 

win210> pvm 

2. Add amelia and fred to your virtual machine 

pvm> add amelia 

1 successful 

HOST DTID 

amelia cOOOO 

pvm> add fred 

1 successful 

HOST DTID 

fred 100000 

3. Check the configuration of your virtual machine 

pvm> conf 


host, 1 data 

format 



HOST 

DTID 

ARCH 

SPEED 

win210 

40000 

SGI 

1000 

amelia 

cOOOO 

SGI 

1000 

fred 

100000 

SGI 

1000 


4. Delete amelia 

pvm> delete amelia 

1 successful 

HOST STATUS 

amelia deleted 

5. Exit PVM console, but PVM daemon is still running 

pvm> quit 

pvmd still running 
win210> 


u 
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PVM Library 


This chapter introduces the PVM library. In writing your application pro- 
grams, you need to call the subroutines in the library to instruct PVM to 
control processes, send information, pack/unpack data, and send/receive 
messages. Many subroutines have pre-defined option values for some 
arguments. These are defined in the include file f pvm3 h and are listed in 
the Appendix. 


call pvmfmytid( tid ) 

This routine enrolls this process with the PVM daemon on its first call, and generates a 
unique ua. You call this routine at the beginning of your program. 


call pvmfexit( info ) 

This routine tells the local PVM daemon that this process is leaving PVM. You call this 
routine at the end of your program. Values of info less than zero indicate an error. 


call pvmfkill( tid, info) 

This routine kills a PVM task identified by tid. Values of into less than zero indicate an 


call pvmfspawnf pname,flag,where,ntask,tids,numt ) 

This routine starts up nt»sk instances of a single process named pnamc on the virtual 
machine. Here are the definition of the other arguments: 


Option Value Meaning 

P VMDEFAULT (0) PVM can choose any machine to start task 

P VMHOST ( 1 ) where specifies a particular host 

PVMARCH (2) where specifies a type of architecture 

PVMDEBUG (4) start up processes under debugger 

PVMTRACE (8) processes will generate PVM trace data 

ere is where you want to start the PVM process. If flag is 0, where is ignored, 

i contains identification numbers of PVM processes started by this routine, 

mt indicates how many processors started; negative values indicate an error. 

Note You should always check tide and mum to make sure all pro- 
cesses started correctly. 
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call pvmrparent ( tid ) 

This routine returns the tid of the process that spawned this task. If the calling process 
was not created with pvirtpawn, then ud-PvmNoParent. 


Dynamic 

Configuration 


call pvmfaddhost( host. Info ) 
call pvmfde!host( host, info ) 

These routines add and delete hosts to the virtual machine respectively. Values of info 
less than zero indicate an error. 



Note Both routines are expensive operations that require the syn- 
chronization of the virtual machine- 


call pvmfinitsend( encoding, bufid ) 

This routine clears the send buffer, and creates a new one for packing a new message. 


encoding 

Encoding Value 

Meaning 


PVMDEFAULT (0) 

XDR encoding if virtual machine config- 
uration is heterogeneous 

This is not implemented 

PVMRAW (1) 

no encoding is done. Messages are sent in 

in PVM v3.2. 


their original format. 

~~~~ ■ 

PVMINPLACE (2) 

data left in place. Buffer only contains 
sizes and pointers to the sent items. 



bufid contains the message buffer identifier. Values less than zero indicate an error. 


Message 

Buffers 


call pvmffreebuf( bufid, info) 

This routine disposes the buffer with identifier bufid. You use it after a message has been 
sent, and is no longer needed. Values of info less than zero indicate an error. 


Packing and 
Unpacking 


call pvmfpackf what, xp, nitem, stride, info ) 
call pvmfunpackf what, xp, nitem, stride, info ) 

These routines pack/unpack your message xp, which can be a number or a string. You 
can call these routines multiple times to pack/unpack a single message. Thus a message 
can contain several arrays, each with a different data type. 
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Note There is no limit to the complexity of the packed messages, 
but you must unpack them exactly as they were packed, 
what indicates what type of data xp is 


nitem 

stride 

info 


STRING (0) 
BYTE! (1) 
INTEGER2 (2) 
INTEGER 4 (3) 


REAL (4) 
COMPLEX8 (5) 
REAL8 (6) 
COMPLEX 16 (7) 


is number of items in the pack/unpack. If xp is a vector of 5, idtem is 5. 


is the stride to use when packing. 

is status code returned by this routine. Values less than zero indicate an error. 


Sending and 
Receiving 


call pvmfsend( tid, msgtag, info ) 

This routine labels the message with an integer identifier msgtag, and sends it immedi- 
ately to the process tid. Values of info less than zero indicate an error. 


call pvmfmcast( ntask, tids, msgtag, info ) 

This routine labels the message with an integer identifier msgtag, and broadcasts the mes- 
sage to all nt«k number of tasks specified in the integer array ads. Values of info less than 
zero indicate an error. 


call pvmfrecv( tid, msgtag, build ) 

This routine blocks the flow of your program until a message with label msgtag has 
arrived from tid. A value of -1 in msgtag or tid matches anything (wildcard). This routine 
creates a new active receive buffer, and puts the message in it. Values of bufld identify 
the newly created buffer; values less than zero indicate an error. 


call pvmfnrecv( tid, msgtag, build ) 

This routine performs in the same way as pvmfrecv, except that it does not block the flow 
of your program. If the requested message has not arrived, this routine returns bufw-0. 
This routine can be called multiple times for the same message to check if it has arrived, 
while performing useful work between calls. When no more useful work can be per- 
formed, the blocking receive pvmftrcv can be used for the same message. 
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call pvmfprobef tid, msgtag, build ) 

This routine checks if a message has arrived; however, it does not receive the message. 
If the requested message has not arrived, this routine returns bufld-0. This routine can 



PVM Library 


be called multiple times for the same message to check if it has arrived, while perform- 
ing useful work between calls. 


call pvmfbufinfo (bufid, bytes, msgtag, tid, info) 

This routine returns information about the message in the buffer identified by bufld. The 
information returned is the actual ragtag, source Ud, and message length in bytes. Values 
of tnfo less than zero indicate an error. 
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Tutorial B 


This chapter shows you how PVM may be applied to your application 
programs through a simple example. The example chosen is the Golden 
Section rule for finding the maximum of a function. You may remember 
it from Math class in high school. Let us review the method and the algo- 
rithm. 


Golden 

Section 


Suppose we want to find the maximum of a curve y-=f(x); where x is 
between the interval aj and a 2 . The points a 3 and a 4 are symmetrically 
placed in this interval, so that 

a 3 -(l-a) a, + aa 2 (EQ l) 

a 4 - a a, + (1-a) a 2 (EQ 2) 

See Figure 1 at left. Golden Section rule requires a to be 0.382. 



The algorithm of finding the maximum is as follow: 



If f(a 4 ) < f(a 3 ) 

If f(a 4 ) > f(a 3 ) 


1 Consider new interval (a t ,a 4 ) 

1 Consider new interval (a 3> a 2 ) 


2 Apply EQ.(l) and (2) again 

2 Apply EQ. (1) and (2) again 

X 

3 Until maximum is reached 

3 Until maximum is reached 


If f(a 3 )-f (a 4 ), the maximum is found 


Figure 1. Interval division for 
Golden Section 



The FORTRAN program (Serial Program) on the opposite page is the 
Golden Section rule that a programmer would write on a normal serial 
computer. Please spend a few minutes to study the flow of the program. 
This simple program consists of two parts, the main (calling) program 
and the function subroutine. The latter has only four lines. 

Note Notice that for each interval (a,,a 2 ), we need to call the 
function evaluation four times to find f(aj), f(a 2 ), f(a 3 ), and 
f(a 4 ). 
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Linear optimization: 


Search for maximum of a x-y curve. 
DIMENSION A(4),FN<4) 


Golden Section rule 


Initial interval 
L « 0 

TOL - 1 . E“3 
A ( 1 ) - 0.4 
A (2 ) - 1.6 
ALPHA - 0.382 


10 CONTINUE 
C 

C Loop begins : 

L - L + 1 


Equations (1) and (2): 


A ( 3 ) - ( 1 . -ALPHA) *A( 1 ) +ALPHA*A(2) 
A ( 4 ) - ALPHA*A(1) + ( 1 . - ALPHA ) *A( 2 ) 
FN ( I ) - F(A(1)J 
FN(2) - F ( A( 2 ) ) I 

FN ( 3 ) - F(A<3)) 

FN { 4 ) - F(A(4) ) J 

WRITE ( 10 , * ) 'A \A(1), 

WRITE ( 10 , * ) * F \FN(1) 

WRITE ( 10 , * ) ' ’ 

ERR - ABS ( FN ( 2 ) - FN ( 3 ) ) 

IF (ERR . LE . TOL) GOTO 999 


Four function evaluations 


/A(l) ,A(2),A(3),A(4) 

/ FN ( 1 ) , FN ( 2 ) , FN ( 3 ) , FN ( 4 ) 


IF ( FN( 4 ) .GT. FN ( 3 ) ) THEN 


A ( 1 ) * Bl 
A ( 2 ) - B2 
GOTO 10 
ELSEIF (FN ( 4 ) 
Bl - A ( 1 ) 
B2 - A( 4 ) 
A( 1 ) - Bl 
A( 2 ) - B2 
GOTO 10 
END IF 
CONTINUE 
STOP 
END 


Iff(a 4 )>f(a 3 ) 


FN ( 3 ) ) THEN 


If f(a 4 ) < f(a 3 ) 


Function evaluation 


FUNCTION F (X ) 

F - TANH (X )/ ( 1 . +X*X ) 

RETURN 

END 
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Recall that in the procedure of finding a new interval, the program calls 
the function evaluation four times serially to get ffaj), f(a 2 ), f(a 3 ), and 
f(a 4 ). We would like to assign four processors to perform the four func- 
tion evaluations simultaneously on the virtual machine. Therefore, we 
modify the Serial Program by writing the main (calling) program as a 
Master program, and the function subroutine as a Slave program. 

The following steps are general guidelines to writing a Master program. 
Please study the steps, and compare them with the program on the oppo- 
site page. Also compare it with the Serial Program. 


1. Include fpvm3.h 

Include this file in your program, you are able to use the PVM 
preset variables; such as pvmdefault, real4, and more, men- 
tioned in Chapter 4 and the Appendix. 

2. Enroll Master to PVM 

Use pvmfniytid(mytid) tO enroll. 

3. Assign virtual processors 

Use the following call to spawn nproc function processes. 

pvm f spawn ( pname, PVMDFAULT, where , nproc, tids , numt ) 

Also tell PVM the name of the Slave program (p name). PVM 
returns tids, the identifier of the nproc processors. 

4. Initialize buffer and pack data 

Use pvmf initsend to clear buffer. 

Use the following routine to pack a real array a of dimension m. 

pvmf pack (REAL4 ,A,m,l,info) 

5. Send message 

Use the following call to send the packed message to the Slave 
process identified by tids. 

pvmfmcast{ nproc, tids,msgtag, info) 


PVM Master 
Guideline 
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Master 

Program 


Linear optimization: 

Search for maximum of a x-y curve. 
PROGRAM MASTER 


include 1 . . /include/fpvm3 . h ' 
DIMENSION A( 4 ) ,FN(4) 
integer tids (0 : 32 ) , who 
character*8 where 
character* 12 pname 


Enroll this program in PVM 
■call pvmfmytid (mytid) 

Start up the four processors 

nproc - 4 Assign four processors 

where * 1 * 1 

— Slave program s name 

pname - function - 

call pvmf spawn ( pname, PVMDFAULT, where, nproc, tids, numt) 
do 20 i-0,nproc-l 

write(*,*) 'tid', i, tids(i) 
continue 


Initial interval 
L - 0 

A ( 1 ) - 0.4 
A ( 2 ) - 1.6 
ALPHA - 0 . 382 
TOL - l.E-3 
ERR - 1. 

CONTINUE 


Loop begins: 
L - L + 1 


Equations (1) and (2): 


(I. -ALPHA) *A(1) + ALPHA*A( 2 ) 
ALPHA*A ( 1 ) + (1 . -ALPHA) *A(2) 


Pack nproc, tids, A, 
and ERR 


Broadcast data to all node programs 
first pack them, then send them 
call pvmf in itsend ( PVMDEFAULT, info) 
call pvmf pack ( INTEGER4 , nproc, 1, 1, info) 
call pvmf pack ( INTEGER4 , tids, nproc, 1, info ) 
call pvm f pack ( REAL4 , A, 4 , 1 , info) 
call pvm f pack ( REAL4 , ERR, 1, 1, info) 


msgtype value matches the one 
received in Slave program ° 


msgtype = 1 

call pvm fm cast (nproc, tids, msgtype, info) 
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6. Wait until messages come from Slaves 

Use pvmfrecv( ) to block until Slaves return function values. 
Make sure value of msgtype matches values coming from Slaves. 


7. Receive and Unpack data 

The sequence of unpacking is the same as the packing in the 
Slave. 


8 . 


Exit PVM 

Use pvmfexit(info) tO exit PVM. 
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6 


c Wait for results from processors 


msgtype value matches the 
one sent from Slave program 

7 ” 


Receive /unpack FN and ‘who’ 
from the 4 processors one by one 100 


msgtype * 2 
do 100 i-l,nproc 

call pvmfrecv( -1 , msgtype, info) 
call pvmf unpack ( INTEGER4 , who, 1 , 1 , inf o) 
call pvm f unpack ( REAL4 , FN(who) , 1, 1, info) 
continue 


WRITE( 10 , * ) f A * ,A(1),A(2),A(3) , A ( 4 ) 
WRITE ( 10 , * ) ' F ',FN(1),FN(2),FN(3),FN(4) 
WRITE( 10 , * ) 1 ' 

ERR - ABS { FN ( 2 ) - FN ( 3 ) ) 

IF ( ERR . LE . TOL) GOTO 999 
C 

c 

IF ( FN( 4 ) .GT. FN ( 3 ) ) THEN 
B1 - A ( 3 ) 

B2 - A ( 2 ) 

A( 1 ) - B1 
A( 2 ) - B2 
GOTO 10 

ELSEIF (FN( 4 ) . LT . FN(3)) THEN 
B1 - A ( 1 ) 

B2 - A(4 ) 

A ( 1 ) - B1 
A( 2 ) - B2 
GOTO 10 
ENDIF 
c 

c Program finished leave PVM before exiting 

999 continue 

... — call pvmfexit( info) 

STOP 

END 
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PVM Slave 
Guideline 


The Slave program is basically the function evaluation program. In order 
to do the function evaluation, it needs information from Master. For 
example, it needs the identity numbers (tids(i), . . . ,tids<4 )) that PVM 
assigns, and the values of a],...,a 4 . 

The following steps are general guidelines to writing a Slave program. 
Please study the steps, and compare them with the program on the oppo- 
site page. Also try to find the connection with the Master Program. You 
may find Figure 1 helpful. 

1. Include fpvm3.h 

Include this file in your program, you are able to use the PVM 
preset variable names; such as pvmdefault, real 4, and more, 
mentioned in all tables in Chapter 4 and the Appendix. 

2. Enroll Slave with PVM 

Use pvmfmytid(mytid) tO enroll. 

3. Identify the parent of this process 

Use the following call to obtain the task identifyer (mtid) of par- 
ent process. This is useful for returning solutions to the Master. 

pvmf parent (mtid) 

4. Receive and Unpack data 

Make sure the value of msgtype matches the one from Master. The 
sequence of unpacking is the same as the packing in Master. 

5. Perform function evaluation 

6. Initialize buffer and pack data 

Use pvmf lnitsend to clear buffer. 

Use the following call to pack a real array f of dimension n. 

pvmf pack (REAL4 , F,n , 1 , info) 

7. Send data 

Use the following call to send the packed message to Master 

pvmfsend (mtid,msgtag, info) 

8. Exit PVM 

Use pvmfexit( inf o) tO exit PVM. 
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program function 


Slave 

Program 


include 1 . . /include/f pvm3 . h f 


integer tids ( 0 : 32 ) , who 
real a(32) 
tor - l.e-3 

Enroll this program in PVM 
call pvmfmyt id (mytid ) 

Get the parent's task id 
.call pvmfparent(mtid ) 

continue 


Unpack the same way as 
Master sends 


Receive data from host 
msgtype » 1 

call pvmfrecv(mt id, msgtype, info) 
call pvmf unpack { INTEGER4 , nproc, 1 , 1 , inf o) 
call pvmf unpack ( 1NTEGER4 , tids , nproc, 1, info) 
call pvmfunpack ( REAL4 , A, 4 , 1, info) 
call pvmfunpack (REAL4 , ERR, 1, 1, info) 


Function 

evaluation 


i f ( err . le . tor ) go to 99 

Determine which processor I am 
do 5 i-0,nproc*l 

if (tids( i ) .eq .mytid) me * i 
continue 
who “ me + 1 

Calculate the function 
X - A (who) 

f - TANH(X)/ ( 1 . +X*X) 


Pack f and 
processor ‘who’ c 


Go to 3Hhd wait for 
another call from master 


"Send the result to Master 
call pvmf initsend (PVMDEFAULT, info) 
call pvmfpack{ INTEGER4 ,who, 1,1, info) 
call pvmf pack ( REAL4 , f , 1, 1 , info ) 
msgtype * 2 

call pvmf send(mtid , msgtype , info) 
go to 3 

Program finished. Leave PVM before exiting 
continue 

call pvmfexit( info ) 

stop 

end 
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Compilation 
and Running 



After you finish your program, it is time to compile and run. Follow the 
steps below to compile your programs. 

1. Make sure you have the correct directory setup 

Follow the advice from Directory Setup in Chapter 1. 

2. Compile the program 

Use the sample Makefile on the opposite page to compile your 
programs. 

Note The Makefile links the PVM library, iibfpvm 3 .a. 


3. Copy executables to all the hosts 

Follow the advice from Directory Setup in Chapter 1, and distrib- 
ute the exeCUtableS tO $HOME/pvm3/bin /ARCH. 

4. Activate PVM 

Activate PVM by entering P vm at UNIX prompt. 

5. Decide the configuration of the virtual machine 

Add or delete hosts to the virtual machine. (Chapter 3) 

6. Quit PVM console 

Leave PVM console (don’t halt daemon) by entering quit at the 
pvm prompt. 
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Makefile # Custom section 

# Set PVM__ ARCH to your architecture type (SUN4, HP9K , RS6K, 
etc. ) 

# if PVM_ARCH - BSD386 then set ARCHLIB - -lrpc 

# if PVM_ARCH - SGI then set ARCHLIB - -lsun 

# if PVM_ARCH “ I860 then set ARCHLIB * -lrpc -lsocket 

8 if PVM_ARCH = IPSC2 then set ARCHLIB * -lrpc -lsocket 

# otherwise leave ARCHLIB blank 

# 

# PVM_ARCH and ARCHLIB are set for you if you use 'aimk'. 


# 


PVM Library 


Make appropriate changes/^ 
for your own path 


PVM_ARCH - SGI 

ARCHLIB - -lsun 

# END of custom section - leave this line here 

# 


PVMDIR - 
PVMLIB - 
SDIR 
BDIR 
XDIR 


/amd/f s02/pub/Iris4d_irix4/nas/pkg/pvm3 
$ ( PVMDIR ) /I ib/$ ( PVM_ARCH ) /I i bpvm 3 . a 

/u /wk /cheun g/p vm 3 /b i n 
$(BDIR)/$(PVM_ARCH) 


CFLAGS - -g - 1 . . /include 

LIBS * $ ( PVMLIB ) $ (ARCHLIB ) 


2 


« SGI, 


F77 * f 77 

FFLAGS - -g 

FLIBS - $ (PVMDIR)/lib/$ ( PVM_ARCH )/l ibf pvm3 .a $ (LIBS) 


default : 


master function 


$ (XDIR ) : 

- mkdir $(BDIR) 

- mkdir $(XDIR) 


clean : 

rm -f *.o bfgs quadfunct 

master: $( SDIR ) /master . f $(XDIR) 

$(F77) $ ( FFLAGS ) -o master $( SDIR ) /master . f $(FLIBS) 
mv master $(XDIR) 

function: $( SDIR) /function . f $(XDIR) 

$ ( F77 ) $ ( FFLAGS ) -o function $ (SDIR)/f unction . f $(FLIBS) 
mv function $(XDIR) 
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Problems and Tips 



PVM is a relatively new piece of software. It is not advanced enough to 
warn you ahead of time before problems come. Here are a couple of 
cases that you may encounter as a beginner. 


Problems 



Can’t activate PVM 



If the message you get, after entering P vm at UNIX prompt, is 

libpvm [pid-1]: Console: Can't start p^md, 

it is possible that the last time you halted PVM daemon, the daemon 
created a residual file /tm P / P vmd.xxxx; where xxxx is an unique 
number for you. Delete this file, and start PVM again. 


If the daemon is running but the PVM console will not start, 
it is possible that you have too many processes running. You have 
to kill all the processes before you re-activate PVM console. 


Note Use P s -ef i username at UNIX prompt to locate your run- 
ning processes. 


Can’t add hosts 



It is possible that there are no links between your local computer and the 
other hosts. Check the following two things: 


Make sure each of your hosts has a . rhosts file in the $home direc- 
tory, and this file points to your local computer. 


Make sure the rhosts file is “read” and “write” protected from oth- 
ers users. 
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Host 

File 


You can create the following file to build the virtual machine without 
activating the PVM console. The addresses must be recognizable by your 
system. 

computerl . address ~"1 

computer2 . address 1 — host file 

computer3 . address 
computer4 . address J 



Note The first machine listed must be the initiating host. 

Note If tasks are to be spawned on specfic systems, the system 
name contained in where ( routine pvm&pawn) must match the 
name in the host file exactly. 

Note If spawning tasks are on the initiating host, use the truncated 
host name. For example, if the full address is 

win210.nas. nasa . gov ] 

use win 2 io instead. This is a bug in PVM v3.2. 


Having the host file ready, enter the following at UNIX prompt, 
Win210> pvmd3 host 
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Notes 



Place to jot down problems. 

If encounter problems, please contact: 
Merritt Smith: mhsmith@nas.nasa.gov 
or 

Samson Cheung: cheung@nas.nasa.gov 
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TABLE 1. 

ARCH names used in PVM. 


ARCH 

Machine 

Note 

AFX8 

AlliantFX 8 


ALPHA 

DEC Alpha 

DEC OSF-1 

BAL 

Sequent Balance 

DYNIX 

BFLY 

BBN Butterfly TC2000 


BSD386 

80386/486 Unix box 

BSDI 

CM2 

Thanking Machines CM2 

Sun front-end 

CMS 

Thanking Machines CMS 


CNVX 

Convex C-series 


CNVXN 

Convex C-series 

native mode 

CRAY 

C-90, YMP, Cray-2 

UNICOS 

CRAYSMP 

Cray S-MP 


DGAV 

Data General Ayiion 


HP 300 

HP-9000 model 300 

HPUX 

HPPA 

HP -9000 PA-RISC 


I860 

Intel iPSC/860 

link-lprc 

IPSC2 

Intel iPSC/860 host 

SysV 

KSR1 

Kendall Square KSR-1 

OSF-1 

NEXT 

NeXT 


PGON 

Intel Paragon 

link -lprc 

PMAX 

DECstation 3100,5100 

Ultrix 

RS6K 

IBM/RS6000 

AIX 

RT 

IBMRT 


SGI 

Silicon Graphics IRIS 

link -Isun 

SUN3 

Sun 3 

SunOS 

SUN4 

Sun 4, SPARCstalicn 


SYMM 

Sequent Symmetry 


TUN 

Staedent Titan 


UVAX 

DEC Micro VAX 
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